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PATTERN  SYNTHESIS  OF  ARRAY 
ANTENNAS  WITH  ARBITRARY 
ELEMENTS  BY  SIMULATED  ANNEALING 
AND  ADAPTIVE  ARRAY  THEORY 


allows  an  exhaustive  and  simultaneous  control  of  lateral 
lobes  and  other  significant  array  design  parameters. 

2.  DESCRIPTION  OF  THE  METHOD 

Let  us  consider  an  array  of  N  elements,  where  F„  is  the 
position  of  the  nth  element  (n  =  1,2, r  is  the  unit 
vector  in  the  far-fieid  observation  direction.  /„(r)  is  the 
radiation  pattern  of  the  nth  element,  and  W=  [>v,,iv,, — 
Wy)Y  are  the  e.xcitations  of  the  array  elements.  The  radiation 
pattern  can  be  expressed  (in  decibels)  as 

/7(r)  =  20  logins ’■■(/(r)l  (1) 

where  i/(r)  is  given  by 

t/(r)  =  (2) 
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ABSTRACT:  A  method  to  cany  out  the  pattern  synthesis  of  aiititrapr 
array  antennas  is  presented.  This  method  uses  the  simulated  annealing 
technique  to  calculate  the  power  interferences  described  in  the  Olen- 
Compton  method.  Its  application  to  a  linear  array  of  16  unequally 
spaced  nonidentical  elements  and  to  a  circular  arc  array  of  25  axial 
dipoles  is  shown  in  order  to  illustrate  its  efficiency  and  accuracy.  ©  1999 
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The  elements  of  the  array  are  going  to  be  used  as  those  of 
an  adaptive  array.  So,  we  assume  that  all  of  the  elements  are 
receiving  a  desired  signal  from  the  direction  r^,  M  interfer¬ 
ence  signals  impinging  on  the  antenna  from  directions  r„ 
(m  =  1, 2, ... ,  M),  and  uniform  thermal  noise.  The  amplitude 
level  of  the  interference  direction  r„  is  a  positive  magnitude 
denoted  by  all  of  them  conforming  to  the  so-called 
interference  spectrum. 

The  theory  of  adaptive  array  establishes  that  the  signal- 
to-interference-plus-noise  power  ratio  (SINR)  can  be  mini¬ 
mized  by  choosing  optimal  excitations  w„.  In  this  synthesis 
method,  an  adequate  selection  of  the  interference  levels  will 
force  the  adaptive  array  technique  to  adjust  a  certain  desired 
sidelobe  level  map  which  allows  us  to  compensate  the  inter¬ 
ference  spectrum  and  to  achieve  the  minimum  SINR. 

According  to  the  theory  of  adaptive  arrays  [3],  the  opti¬ 
mum  excitation  coefficients  obtained  to  maximize  the  ratio 
between  the  desired  and  the  interference-plus-noise  signals 
(SINR)  are  given  by 


Key  words:  adaptive  theory;  antenna  array  pattern  synthesis 
1.  INTRODUCTION  • 

In  recent  years,  most  of  the  work  in  the  pattern  synthesis  of 
arrav  antennas  has  been  focused  on  uniformly  spaced  arrays 
with  isotropic  elements  (1].  In  [2],  Olen  and  Compton  pro¬ 
posed  a  synthesis  algorithm  for  arbitrary  arrays  (with  differ¬ 
ent  elements  located  in  arbitrary  positions)  which  is  based  on 
the  adaptive  array  theory  described  in  (3).  The  method  has 
some  convergence  problems  since  the  adjustment  for  the 
interference  power  (namely,  a  parameter  that  controls 
the  interference  spectrum  intentionally  imposed  to  adjust 
the  secondary  lobe  map)  is  more  likely  qualitative  than  quan¬ 
titative.  An  improvement  in  the  convergence  of  this  method 
by  means  of  the  stability  criterion  of  Routh-Hurwitz  also  has 
been  presented  in  (4].  And  a  new  extension  of  the  method, 
based  on  a  linear  expansion  of  the  gain  near  the  interference 
levels,  has  been  developed  in  (5).  The  main  drawback  of  the 
method  is  the  lack  of  exhaustive  control  over  every  lateral 
lobe  because  those  whose  sidelobe  level  is  lower  than  the 
desired  one  are  not  considered.  As  a  consequence  of  that,  a 
beam  widening  and  gain  reduction  is  caused.  Furthermore, 
there  is  no  control  over  other  array  parameters  such  as  the 
excitation  dynamic  range  and  adjacent  element  on  excitation 
ratio.  In  this  paper,  we  present  a  method  based  on  [5],  but 
using  a  simulated  annealing  technique  for  the  computation  of 
interference  levels.  The  main  advantage  is  that  this  method 


JF  =<!>-’•  £7*  (r^)  (3) 

where  <l>  is  a  covariance  matrix  defined  as 

4>  =  /+  L  (4) 

m»  I 

where  I  denotes  the  identity  matrix,  superscript  (*)  denotes 
the  conjugate,  and  (T)  the  transpose. 

The  novelty  of  the  proposed  method  consists  in  the  selec¬ 
tion  of  the  interference  levels.  A  cost  function,  which  takes 
into  account  the  sidelobe  level  of  the  patterns  as  well  as 
other  significant  parameters  of  the  array  such  as  the  dynamic 
range  of  excitations,  is  defined  and  minimized  by  a  simulated 
annealing  technique  (6].  The  cost  function  we  have  defined 
can  be  written  as 


C  =  ci-  t  [SLL.,  '-SLL,,Y  +  c. 


p- 1 


w_ 


(5) 


where  SLLp  „  is  the  obtained  level  of  the  pih  lobe,  SLLp  j  is 
the  desired  one.  \w^/w^\  is  the  excitation  dynamic  range, 
and  C;  are  weights  to  control  the  importance  given  to  each 
term  of  (5)  in  the  optimization  procedure.  The  potentiality  of 
the  method  is  that  any  parameter  of  interest  (directivity, 
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'  'i  beam  width,  etc.)  can  be  easily  considered  by  including  it  in 
the  cost  function. 

1 

,  .  --i  3.  EXAMPLE  DESIGNS 

'  ^  3  1  Linear  Array  with  Noniiniformly  Spaced  Elements.  In  order 

to  illustrate  the  method,  we  have  synthesized  ar«y  oj 
•  dipoles,  placed  along  the  z-axis  at  distances  d  with  different 

lengths  l„  and  different  orientations  (denoted  by  the  ang*e  ■„ 
'  beween  the  dipole  and  the  array  axis).  The  U  vector  for  the 

:  nth  element  is 

- '  ^  cos(7t(.  sin(0  +  T„)]  -  cosCtt/,) 


•  exp  j2iT  52  ® 

1 1-0 


The  desired  pattern  has  a  main  beam  and  two  adjacent 
lobes  with  a  level  of  -35  dB.  A  level  of  -20  dB  is  requjed 
for  the  rest  of  the  lobes.  A  discrete  spectrum  of  58  interfer¬ 
ence  directions  is  distributed  across  the  sidelobe  region.  By 
minimizing  the  cost  function,  a  set  of  interference  levels  ^ 
adjusted;  then  the  excitations  are  computed  by  (3)  and  (4). 
Table  1  shows  the  geometric  parameters  of  the  array  ele¬ 
ments  and  their  corresponding  optimal  excitations.  ^*8“" 
shows  the  far-field  pattern  computed  by  using  Eqs.  (1)  and 
(2).  demonstrating  the  efficiency  of  the  proposed  technique. 

3.2  Circular  Arc  Array.  Let  us  consider  an  anay  antenna  of 
S  dipoles  in  the  z-direction,  uniformly  distributed  on  a 
perfectly  conducting  cylinder,  conforming  an  arc  of  radius  R 
between  angles  -4>o  <^o  a*  shown  in  Figure  2-  In  our 

case,  the  selected  values  are  Af  ”  25,  /?  -  5.7jA,  and  “ 
60*  to  produce  a  separation  between  adjacent  elements  of 
0  5A  We  have  introduced  72  interference  directions  uni¬ 
formly  distributed  for  5*  l«^l  ^  180*.  In  this  ca^.  we  want 
to  synthesize  an  asymmetric  pattern  with  -40  dB  sidelobes 
on  one  side  of  the  main  lobe  and  a  tapercd-sidelobe  structure 
from  -  25  to  -  35  dB  on  the  other  side.  The  factor 
has  to  be  minimized. 


TABLE  1  Parameters  of  the  Unear  Array  with  1 6  Arbitrary 

Dipoles  and  Element  Excitations  Obtained  After 

the  Optimization  _ _ _ 

n  _ 

1  0.25  0.0  OiO 

1  ms  0  J  0  JO  0.031  •  e‘' 

-t  nod  SO  OJa  O.lOo-e" 

3  0.24  a.u  nnax  ,/(+0J*5) 

t  n-jn  —32  0J4  0.086 

s  n  -JX  —  3  2  0.60  0.184  • 

5  0.26  «  ,/(+0.097) 
c  n?7  10  0.45  0.122 -e'' 

7  0.23  1-0  0.46  0.190  e 

o  0  24  0  0  OJO  0.129  • 

02^  00  OJO  0.218. 

14  n  10  1 0  0-57  0»151  f 

-*  A  *9*5  —XI  0.6o  0.062 

15  0.22  f-r  A,„.,/{-0J)«) 


4 

0.20 

-32 

5 

0.26 

.  -3.2 

0.27 

10 

7 

0.23 

1.0 

V  8 

0.24 

0.0 

0.25 

0.0 

10 

0.21 

7.0 

11 

0.28 

6.0 

•  12 

0.30 

4.4 

13 

0.29 

0.0 

14 

0.19 

1.0 

15 

0.22 

-11 

Mm  16 

0.22 

3.0 

-90  -75  -60  -45  -30  -15  0  15  30  45  60  75  90 

theta  (degrees) 

Figure  1  Radiation  pattern  of  the  synthesized  linear  array  with  16 
arbitrary  dipoles 


(|)=-60* 
N,  ■ 


Rgure  2  Geometry  of  the  circular  arc  array 
The  power  pattern  of  this  array  can  be  expressed  as 

f(<i)  -  E  I„e>“-e^'‘‘^’^°'^*'*‘'^FE(.d>  -  <(>„)  (7) 

/!■  I 

where  1„  and  a„  are  the  amplitude  and  phase,  respectively,  of 
the  excitation  of  the  nth  element,  and  FEW  is  the  pattern 
of  the  element  in  the  .^fT-plane.  In  conformal  arrays,  this 
term  is  of  great  importance  because  each  element  is  indepen¬ 
dently  oriented.  The  element  pattern  used  in  this  paper  is 
that  proposed  by  Jiao  et  al.  in  (7],  based  upon  results  pre¬ 
sented  by  Herper,  Hessel,  and  Tomasic  in  [8,91  for  a  cylindri¬ 
cal  array  of  axial  dipoles  on  a  perfectly  conductmg  cylmder. 
The  FE  factor  for  a  separation  of  0.5  A  is  given  by 

FE{d,  d>)  -  3  sin  0[1  +  2max(cos  d>,  -1/2)1  (8) 

where  the  a-dependence  is  ipored  in  our  case  because  we 
only  consider  pattern  synthesis  in  the  ^LY-plane. 

Figure  3  shows  the  obtained  array  pattern,  again  showing 
the  good  performance  achieved.  Table  2  shows  the  excitation 
coefficients  of  this  array,  whose  has  been  mim- 

mized  to  a  value  of  8. 
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phi  (degrees) 


Rgure  3  Radiation  pattern  of  the  synthesized  circular  arc  array  of 
25  axial  dipoles 


TABLE  2  Excitations  Obtained  for  the  Synthesis  of  the 
Circular  Arc  Array  of  25  Axial  Dipoles 


rt 

n 

1 

0.264  • 

14 

0.840  • 

2 

0.180 

15 

0.766 

3 

0.220  •  e<^* 

16 

0.698 

4 

0.238  • 

17 

0.604 

5 

OJll  ■ 

18 

0.507  •  e^‘-  '•^®®’ 

6 

0.417  • 

19 

0.450 

7 

0.486  • 

20 

0o77-e^^‘'***’ 

8 

0.595  •  e^^- 

21 

0.296  • 

9 

0.629  • 

22 

0.209  •  e>‘“®®”> 

10 

0.762  •  ffA+*®®®’ 

23 

0.121 

11 

0.826  •  e'^* 

24 

0.103 

12 

13 

0.797  •  e^* 

0.844  •  e^* 

25 

0.103  •  e'f* 

4.  CONCLUSIONS 

A  novel  technique  for  synthesizing  arrays  by  using  the  adap¬ 
tive  array  theory  and  simulated  annealing  technique  has  been 
presented.  The  method  is  based  on  the  optimum  selection  of 
a  set  of  interference  levels  which  are  combined  with  an 
adaptive  array  scheme  to  minimize  the  signal-to-interference- 
plus-noise  power  ratio,  allowing  us  to  obtain  the  desired 
sidelobe  map.  The  potentiality  of  the  method  is  that  it  is  able 
to  take  into  account  any  significant  parameter  of  the  array 
such  u  the  excitation  dynamic  range,  gain,  beam  width,  etc. 
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ABSTRACT:  A  daign  criterion  for  the  optimization  of  bandpass  optical 
receivers  based  on  lossless  matching  networks  has  been  established.  A 
minimum  value  for  the  input-equivalent  noise  current  under  the  con¬ 
straint  of  fixed  transimpedance  gain  and  an  expression  for  detuning 
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1.  INTRODUCTION 

Bandpass  optical  systems  are  becoming  more  and  more  im¬ 
portant  owing  to  the  need  for  bandpass  low-loss  transmission 
channels  (i.e.,  CATV  AM-SCM  systems,  microwave  links  for 
radar  systems).  In  these  applications,  the  bandwidth  require¬ 
ments  are  quite  moderate  with  respect  to  the  capability  of  the 
optical  fiber,  but  the  system  designer  has  to  take  into  account 
signal  degradations  produced  along  the  link:  the  RIN  of  the 
transmitting  laser,  nonlinear  distorsions  in  multioctave  and 
SCM  systems,  fiber  attenuation  and  dispersion,  shot  noise, 
and  preamplifier  noise  figure  [1-3].  The  noise  influence  and 
intermodulation  products  power  are  inversely  related  by  the 
choice  of  the  optical  modulation  index  {3],  and  the  system 
design  is  strong^  affected  by  this  tradeoff.  An  improvement 
of  receiver  sensitivity  allows  us  to  achieve  in  a  simpler  way 
system  requirements  in  terms  of  one  or  more  of  the  following 
parameters:  1)  maximum  number  of  users  connected  to  a 
single  node  in  a  star  architecture,  2)  minimum  required 
transmitted  optical  power,  and  3)  maximum  length  of  the 
single  portion  of  the  link. 
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Abstract:  A  method  of  pattern  synthesis  of  arbitrary  linear  arrays  is  presented.  This  method  uses  the 
simulated  annealing  technique  to  calculate  the  power  interferences  m  the  Olen-Compton  method.  Application 
to  a  linear  array  of  16  unequally  spaced  non-identical  elements  is  shown. 


l.Introducdon 

El  problema  de  la  sintesis  de  patrones  de 
radiacion  ha  tenido  un  gran  interes  en  los  liltimos 
anos,  pero  la  mayoria  de  los  trabajos  se  centran  en 
arrays  lineales  con  elementos  isotropicos  y 
equiespaciados.  De  entre  estos,  cabe  destacar  el 
metodo  de  Orchard-Elliott  [1].  Olen  y  Compton  [2] 
propusieron  un  algoritmo  de  sintesis  que,  utilizando  la 
teoria  de  arrays  adaptativos  [3],  permitia  sintetizar 
arrays  con  elementos  distintos  y  dispucstos  en 
posiciones  arbitrarias.  El  mayor  problema  de  este 
metodo  radica  en  el  hecho  de  que  el  ajuste  del 
parametro  que  mide  los  niveles  de  interferencia 
necesarios  para  el  control  de  la  topografia  de  lobulos 
de  los  diagramas  de  radiacidn  se  realiza  mediante 
prueba  y  error,  lo  que  conlleva  grandes  dificultades  a 
la  hora  de  conseguir  que  el  algoritmo  converja. 
Posteriormente,  se  introdujo  una  extension  de  dicho 
algoritmo  que,  basada  en  el  criterio  de  estabilidad  de 
Routh-Hurwitz  mejoraba  la  convergencia  [4].  En  otro 
trabajo,  se  consiguio  mejorar  los  rcsultados  obtenidos 
mediante  un  desarrollo  lineal  de  la  ganancia  en  tomo 
a  los  niveles  de  interferencia.  En  dicho  metodo,  dichos 
niveles  se  calculan  resolviendo  un  sistema  de 
ecuaciones  lineales,  lo  que  le  confiere  una  gran 
velocidad  de  convergencia  [5].  Sin  embargo,  dicho 
metodo  no  permite  tener  un  control  exhaustive  sobre 
cada  uno  de  los  Idbulos  laterales,  ya  que,  con  el  fin  de 
asegurar  la  convergencia,  se  descartan  aquellos 
lobulos  cuyo  nivel  es  inferior  al  deseado.  Esto 
provoca  un  aumento  del  ancho  de  haz  que  conlleva 
una  perdida  en  la  ganancia.  Ademas,  dicho  metodo  no 
tiene  control  sobre  otros  parametros  del  array,  como 
son  el  rango  dinamico  o  la  variabilidad  de  las 
excitaciones  entre  elementos  colindantes.  En  este 
trabajo  se  presenta  un  metodo,  basado  en  [5]  pero  que 
utiliza  la  tecnica  de  “simulated  annealing”  para  el 
calculo  de  los  niveles  de  interferencia.  El  nietodo 
permite  tener  un  control  exhaustive  sobre  los  lobulos 
del  diagrama  de  radiacidn  asi  como  otros  parametros 
de  interes  en  el  diseno. 


2.  Descripcion  del  metodo 

Considerese  un  array  de  N  elementos.  Sea  f„el 
vector  de  posicion  del  elemento  n-esimo  (n=l,  2,..,  N), 
f  el  vector  unitario  en  la  direccion  de  un  punto  en 
campo  lejano,  /„(f)  el  diagrama  del  elemento  n-esimo 
y  W=[m'„vVj,...,  el  vector  con  las  excitaciones  de 
cada  elemento  radiante.  Entonces,  el  vector  director 
del  array  viene  dado  por; 


U(f)  =  [/;(r)  /iCr)  //?)  V  (1) 


mientras  que  el  diagrama  de  radiacion,  expresado  en 
dB's,  viene  dado  por  la  siguiente  expresibn: 

p(f)=201og|W^U(f)l  (2) 

Supongase  que  el  array  es  adaptativo  y  que  las 
senales  recibidas  por  los  elementos  incluyen  una  senal 
deseada  procedente  de  la  direccion  -f  ^  ,  M  senales  de 
interferencia,  en  ciertas  direcciones  denotadas  por 
-r  (772=1, 2, ....  M),  y  ruido  termico  de  igual 
magnitud  para  todos  los  elementos.  De  acuerdo  con  la 
teoria  de  antenas  adaptativas  [3],  para  maximizar  la 
razon  de  la  senal  deseada  frente  a  las  senales  de 
interferencia  y  ruido  (SINR),  las  excitaciones  optimas 
deben  verificar  la  siguiente  expresion: 

W=3)'‘U*(fj)  (3) 

en  donde  la  matriz  de  covarianza  se  calcula  mediante 
la  siguiente  expresion: 


en  donde  es  el  nivel  de  la  interferencia  771-esinia, 
que  debe  ser  positiva.  Se  observa  que,  ajustando 
dichos  niveles,  se  modifica  "W  y  con  ello  el  diagrama 
de  radiacion  adaptado.  Estos  niveles  se  calculan 
minimizando  una  fiincion  de  costo  que  tiene  en  cuenta 
el  nivel  de  lobulos  laterales  de  dicho  diagrama,  asi 
como  otros  parametros  de  interes  en  la  optimizacion. 
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Esta  funcion  viene  dada  por  la  siguiente  expresion 

C=c^-'t  (5) 

/  =  1 

en  donde  5Z,Z,,  „  y  SLL,-,,  denotan  el  nivel  en  porencia 
del  lobulo  i-esimo  obtenido  y  deseado 
respectiyamente,  es  el  rango  dinamico  de 

las  excitaciones,  mientras  que  c,.  son  los  factores  de 
peso  que  permiten  darle  mas  importancia  a  ciertos 
parametros  en  la  optimizacion.  Esta  funcion  es 
minimizada  mediante  la  tecnica  de  “simulated 
annealing”,  la  cual  permite  optimizar  cualquier 
paiimetro  de  interes  (directividad,  ancho  de  haz,  etc.,) 
sin  mas  que  incluirlo  en  la  funcion  de  costo. 


3.  Ejcmplo  de  diseno 

Para  ilustrar  el  metodo  se  ha  utilizado  un  array 
lineal  de  16  elementos  distintos  no  equiespaciados.  El 
vector  director  del  elemento  n-esimo  viene  dado  por 


cos[TC  sin(6  +  rj]  -cos(7: /J 
cos(e  +  r„) 

^  n-l  ^ 

:  exp  sm(0) 

V  '=0 


(6) 


en  donde  0  es  el  angulo  medido  desde  el  “boresight” 
del  an'ay,  /„  es  el  tamano  del  dipolo  7i-esimo 
expresado  en  longitudes  de  onda,  es  el  angulo  que 
forma  este  con  el  eje  del  array  (en  grados)  y  es  su 
espaciado  tambien  en  longitudes  de  onda.  Los  valores 
que  se  ban  utilizado  se  muestran  en  la  tabla  I  (rfo=0). 

Se  desea  sintetizar  un  diagrama  de  tipo  suma  con 
los  dos  lobulos  interiores  a  cada  lado  del  haz  principal 
a  un  nivel  de  -35  dB  y  el  resto  a  -20  dB.  Para  este 
proposito,  se  ha  utilizado  un  total  de  58  interferencias 
distribuidas  sobre  la  region  de  lobulos  laterales  del 
diagrama  de  radiacion.  Tras  la  minimizacion  de  la 
funcion  de  costo  se  obtienen  los  niveles  optimos  de 
cada  una  de  ellas  que  mas  se  ajustan  a  las 
especificaciones  del  diseno.  Las  excitaciones 
asociadas  se  calculan  utilizando  las  expresiones  (3-4) 
y  se  muestran  en  la  tabla  I  (con  el  argumento 
expresado  en  radianes).  Finalmente,  el  diagrama  de 
radiacion  se  muestra  en  la  fig.  1 . 


n 

h 

Si 

1, 

0.25 

0.0 

0.50 

0.093-ej'’°°’''’ 

2 

0.25 

0.5 

0.50 

0.031-ej‘'®’‘’' 

3 

0.24 

5.0 

0.55 

0.106-ej'‘'°”®> 

4 

0.20 

-32 

0.54 

0.086-ej‘*‘’-’®*’ 

5 

0.26 

-3.2 

0.60 

0.184-ej‘-®°“’ 

6 

0.27 

10 

0.45 

0.122-e''*°”’' 

7 

0.23 

1.0 

0.46 

0.190-e^‘"®'^^ 

8 

0.24 

0.0 

0.50 

0.129-e^<-““’-> 

9 

0.25 

0.0 

0.50 

0.218-ej<*°^-'’ 

10 

0.21 

7.0 

0.51 

0.155-ej'-‘’’^*’ 

11 

0.28 

6.0 

0.47 

0.177-e^‘''°'*” 

12 

0.30 

4.4 

0.48 

0.07  l-e^'-'’”'’ 

13 

0.29 

0.0 

0.61 

14 

0.19 

1.0 

0.57 

0.151-ej''®'*®’ 

15 

0.22 

-2.1 

0.65 

0.062-e^‘"®’”> 

16 

0.22 

3.0 

- 

0.153-e^''°°’=’ 

_ f 

Tabla  I.  Parametros  del  array  y  excitaciones 
obtenidas  en  la  sintesis. 
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Abstract-A  method  is  described  which  optimizes  the  performance  of  an¬ 
tenna  arrays  with  defective  elements,  by  numerically  Snding  a  new  excitation 
for  the  remaining  ones.  Two  approaches  for  recovering  the  patterns  have 
been  investigated.  In  the  first,  new  excitations  are  obtained  by  changing 
their  amplitude  distribution,  whereas  in  the  second,  only  the  phase  distri¬ 
bution  is  modified  while  keeping  the  original  amplitudes.  These  optimum 
values  for  new  excitations  are  calculated  by  the  application  of  the  simulated 
annealing  technique.  Examples  for  both  sum  and  shaped  beam  patterns  axe 
presented  to  ulustrate  the  optimal  results  obtained. 

1.  INTRODUCTION 

In  antenna  arrays,  element  failures  destroy  symmeciy  and  cause  sharp  vari¬ 
ations  in  the  field  intensity  across  the  array  apertme,  increasing  the  sidelobe 
level  of  the  power  pattern.  Recently  there  has  been  an  increased  interest 
in  developing  methods  to  improve  the  patterns  of  arrays  in  the  presence  of 
failed  or  partially  failed  elements.  A  recent  work  replaces  the  signals  from 
failed  elements  in  a  digitally/  beamformed  receive  array  [ll  .  The  method  is 
restricted  for  receive  and  is  only  applied  to  linear  arravs.  Another  approach 
uses  neural  networks  to  produce  a  new  array  transfer  function  and  oners 
advantages  in  direction  Snding  in  a  single  signal  em-ironment  [2].  The  array 
failure  correction  has  also  been  performed  by  means  of  a  quadratic  optimiza¬ 
tion  to  achieve  a  suboptimum  beam  pattern  closely  matching  the  original 
desired  specification  [31.  In  addition,  there  are  several  resynthesis  methods 
based  on  calculating  a  new  set  of  excitations  for  the  remaining  elements  of 
the  array.  These  excitations  axe  obtained  by  using  a  genetic  algorithm  to 
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minimize  a  given  cost  function  [4]  or  via  a  conjugate  gradient  method  [5]. 
Although  chese  resynthesis  methods  aire  the  only  ones  appropriate  for  trans¬ 
mitting,  they  require  to  modify  both  the  amplitude  and  phase  distribution 
of  the  aperrure,  which  may  increase  the  hardware  complexity.  The  advent  of 
low-power  solid-state  microwave  amplifiers  has  stimulated  interest  in  phase- 
only  antenna  pattern  symthesis  in  recent  years,  because  it  is  easier  and  less 
expensive  to  modify  dynamically  the  phase  of  an  iLF -signal,  rather  than 
its  power  level  [6,  7].  Frey  and  EUiott,  have  shown  that  it  is  possible,  for  a 
fixed  amplitude  distribution,  to  obtain  patterns  with  asymmetric  sidelobes 
(with  respect  to  the  main  beam)  by  phase-only  changes  [8]. 

In  this  paper,  we  have  used  .the  simulated  annealing  technique  [9],  to 
minimize  by  amplitude-only  synthesis  or  phase-only  changes  a  cost  function 
that  is  the  square  of  the  difference  between  the  obtained  and  desired  sidelobe 
levels  in  a  set  of  o  -cuts  for  the  patterns  of  planar  arrays  with  failed  elements. 
In  this  way  the  technique  performs  an  improvement  of  the  degraded  pattern 
in  terms  of  the  sidelobe  level.  Throughout  this  paper,  the  elements  are 
assTimed  to  be  completely  failed,  which  is  the  most  useful  and  interesting 
case  since  known  errors  can  be  quite  easily  corrected  [Ij.  The  main  advantage 
of  using  the  simulated  annealing  technique,  is  to  eliminate  the  possibility  of 
getting  the  cost  function  stuck  in  local  minima,  one  of  the  main  problems 
of  the  conventional  gradient-based  methods.  Examples  using  a  10  x  10 
rectangular  grid,  circular  boundary  planar  array  with  three  elements  failures 
are  presented  for  sum  patterns.  We  also  show  results  relating  to  shaped- 
beam  patterns  using  a  16  x  16  rectangular  grid,  circular  boundary  planar 
array  with  two  defective  elements,  which  are  specially  difficult  to  reoptimize 
since  such  distributions  exhibit  more  variability  chan  those  associated  with 
sum  and  difference  patterns.  In  fact,  in  the  previous  literature,  there  are  no 
works  of  failure  correction  which  show  results  of  shaped-beam  patterns. 

2.  DESCRIPTION  OF  THE  METHOD 

The  pattern  in  the  far  field  for  a  planar  array  of  W  isotropic  elements  laying 
on  the  XY  -plane,  with  an  arbitrary^  grid  and  boundary  is  given  by: 


F(5,C))  =  y^  In  sxpijk sin. d[xn.  cos 0 -^  yn sin. (pi)  (1) 

where  {xn,yn)  a--e  the  positions  of  the  radiacing  elements  and  I- 
=  An  exp(yan)  their  complex  excitations. 

For  computing  convenience  it  is  useful  to  express  (1)  as  follows:  F{j}) 
represents  F  {9^3,  (ps)  where  drs  is  a  particular  value  of  5  in  the  cut  ©3 .  In 
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a  similar  wav,  reoresents  the  exponential  factor  exp^jifcsin^far-riCOSffl  . 
y„sinpj) .  Assigning  the  index  y  to  a  specific  value  of  n  and  p  to  the 
(^rsj  <Ps)  direction  we  obtain: 

fo>)  =  Ew(p)  (-) 

s 

G(p) ,  the  power  in  dB,  is  defiined  using  the  relation: 

a(p)  =  10  logio(BQ’)]'’  (3) 

The  cost  function  C  is  defined  as  the  sum  of  the  square  of  the  differences 
between  the  actual  pattern  G(p)  and  the  desired  pattern  D(v) : 

C^=E[G(p)-B(p)P  =  Ei£{?)f  w 

P=1  p=l 


where  po  is  the  number  of  points  considered.  The  cost  function  is  always 
positive  and  depends  on  a  finite  number  of  variables  (planar  array  excita¬ 
tions).  The  minimization  process  applied  to  the  cost  function  by  chang¬ 
ing  their  variables  is  equivalent  to  improve  the  pattern.  For  the  case  of 
amplitude-only  synthesis,  the  amplitudes  of  non-defective  elements  axe  used 
as  variables,  and  given  by  Ag  =  .4°  -r  SAg  .  Phase-only  changes  are  accom¬ 
plished  changing  the  phase  variables  of  non-defective  elements  as  given  by 
=  0:°  -i-  Sag  ,  keeping  the  amplitudes  constant.  In' other  words,  the  prox¬ 
imity  of  the  actual  pattern  to  the  desired  one  is  measured  by  the  numerical 
value  of  a  cost  function.  The  minimization  of  this  function  is  performed  by 
means  of  the  simulated  annealing  method  [9].  Using  tins  method,  it  is  pos¬ 
sible  to  include  in  the  cost  function  some  other  parameters  of  interest,  such 
35  directivity,  dynamic  range  |-4may/Amini  for  amplitude-only  synthesis  or 
IcJtnax  —  aminl  for  ohase-oniy  changes,  because  it  is  not  necessar:/  to  know  its 
dependence  with  the  cost  function  as  it  happens  in  gradient  based  methods. 

3.  SliVIUL-ATED  ANNEALING:  AN  INTRODUCTION 

The  simulated  annealing  technique  is  based  on  an  analogy  with  thermody¬ 
namics,  specifically  on  how  the  slow  cooling  of  a  system  allows  it  to  find  a 
minimum  energy  state.  The  method  is  suitable  for  optimization  problems  of 
large  scale,  especially  ones  where  a  desired  global  extremum  is  hidden  among 
many,  poorer,  local  e.xtrema.  One  of  the  main  ad’/antages  of  the  algorithm 
is  its  easy  implementation  and  use.  Unlike  gradient- based  methods,  this 
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technique  only  requires  to  know  the  function  -values  and  not  its  derivatives. 
This  means  that  any  design  parameter  can  be  optimized  by  only  including 
it  in  the  cost  function.  The  algorithm  always  accepts  a  do-wnhill  step,  but 
somerimes  accepts  an  uphill  one.  In  this  way,  the  algorithm  can  avoid  get 
stuck  in  local  minima,  one  of  the  main  problems  of  the  gradient-based  meth¬ 
ods.  The  probabilit>'  of  an  uphill  excursion  depends  on  a  control  parameter 
T  (analog  of  temperature)  and  the  values  of  the  cost  function  E  (analog 
of  energy). 

The  flowchart  of  the  simulated  annealing  technique  is  sho-wn  in  Figure  1. 
The  algorithm  begins  with  an  encoding  of  the  -variables  of  the  problem. 
Then,  the  Lniticd  temperat-ure  is  fixed  as  well  as  its  decrement  in  each  itera¬ 
tion  and  the  initial  solution  for  the  variables.  After  that,  this  configuration 
is  e-valuated.  beins:  .  the  cost  function  ’/alue  obtained. 

In  the  next  step,  a  random  perturbation  in  any  of  the  variables  is  per¬ 
formed,  obtaining  a  new  configuration  with  a  cost  value  equal  to  E^  ■  The 
probabihty  of  an  energy  state  E  is  given  by  the  so-called  Boltzmann  distri¬ 
bution; 

Prob(£)  ~  exp{-E/kT)  (5) 

where  k  is  the  Bolzmann’s  constant.  The  probability  p  of  choosing  the 
new  configuration  is  then  given  by; 

p~exp[-{E2-Ei)/kT]  (6) 

This  implies  that  if  the  cost  fmction  decreases,  Eo  -will  be  lower  that  Ei , 
and  therefore  p  -will  be  greater  than  1 ,  which  means  that  the  transition  is 
accepted.  It  is  obvious  that  if  the  system  only  accepts  those  transitions  ver¬ 
ifying  this  condition  (that  is,  accepting  only  do-wnhiil  steps),  the  algorithm 
could  easy  get  stuck  in  a  local  minima.  The  strategy  employed  to  avoid  this 
consist  on  accepting  those  transitions  whose  probabilit]/-  p  is  greater  than  a 
random  ^mlue  x  ,  calculated  in  the  range  (0,  1).  In  this  -way,  the  algorithm 
accepts  not  only  all  the  do-wnhill  steps  but  some  uphill  ones  too,  trying  to 
avoid  local  minima. 

If  the  new  configuration  is  rejected,  a  new  random  perturbation  is  intro¬ 
duced  and  the  process  above  described  is  repeated.  If  not,  the  temperature 
is  slightly  decreased  (in  a  small  amount  which  guarantees  a  slow  cooling  of 
the  system)  and  the  process  is  iterated,  generating  a  new  random  configu¬ 
ration.  The  process  ends  when  the  temperature  reaches  a  pre-'dousiy  fixed 
value. 

4.  APPLICATION 

We  have  synthesized  a  suin  pattern  ’using  a  10  x  10  element  planar  array 
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Figure  1.  Flowchart  of  the  simuiated  annealing  algorithm. 
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in  a  rectangTilar  grid  and  circular  boundary  with  A/2  interelement  spacing. 
The  corresponding  aperture  distribution  was  obtained  by  sampling  a  circular 
Taylor  distribution  [10]  with  a  sidelobe  level  of  —35  dB  and  n  =  5  .  The  sum 
power  pattern  without  element  failures  has  a  peak  sidelobe  level  of  —34.5 
dB  and  is  plotted  in  Figure  2,  where  u  =  sin  &  •  cos  ©  and  v  =  sin  0  -  sin  ©  . 

Without  loss  of  generality,  let  us  assume  that  the  array  has  three  defective 
elements  located  in  the  first  quadrant,  as  is  shown  in  Figure  3.  AVhen  these 
elements  fail,  the  peak  sidelobe  level  rises  up  to  —29.3  dB.  This  degraded 
pattern  is  plotted  in  Figure  4. 

In  the  first  case,  the  optimization  of  the  pattern  has  been  performed 
calculating  a  new  set  of  amplitudes  for  the  non-defective  elements.  Hence, 
the  number  of  unknowns  to  be  solved  by  the  simulated  annealing  technioue 
is  73.  The  obtained  pattern  has  a  peak  sidelobe  level  of  — 34.S  dB  and 
is  shown  in  Figure  5.  The  results  indicate  that  it  was  possible  to  reduce 
the  sidelobe  level  of  the  array  with  failed  elements  down  to  the  design  level. 
However,  as  it  happened  in  [5],  this  reduction  in  the  sidelobe  level  produces  a 
corresponding  increase  in  the  beamwidth  of  the  pattern,  that  has  an  elliptical 
cross  section,  yielding  a  drop  in  gain.  This  result  could  be  improved  by 
including  the  directivity  in  the  cost  function,  at  the  expense  of  substantially' 
increasing  of  the  computational  time. 

In  a  second  approach,  a  excitation  phase  for  each  non-defective  element  is 
calculated  while  keeping  the  original  amplitude  distribution.  The  obtained 
pattern  using  this  technique  has  a  peak  sidelobe  level  of  —30.5  dB  and  is 
shown  in  Figure  6.  This  yields  an  improvement  of  only  1.2  dB  over  the 
degraded  pattern,  which  may  nor  justify  the  use  of  the  phase  distribution 
required.  However,  it  has  been  found  that  these  results  can  be  substantially 
improved  if  we  fix  asymmetric  lobes  in  the  optimization.  This  allows  to 
recover  almost  the  original  sidelobe  level  in  one  half  of  the  oattem  at  the 
expense  of  increasing  the  sideiobes  in  the  other  half.  Using  this  asymmetry 
in  our  example,  it  was  possible  to  reduce  the  peak  sidelobe  level  to  —33.5 
dB  in  the  ©-range  (0°-180°) ,  since  the  desired  level  was  —35  dB,  and 
—27.3  dB  in  the  ©-range  (180°-360°),  in  order  to  reach  —28  dB.  The 
obtained  pattern  is  shown  in  Figure  7.  Note  that  the  oattem  does  not 
show  a  significant  increase  in  the  beamwidth.  unlike  the  amolitude  only 
synthesis. 

The  technique  can  be  also  applied  to  shaped  beam  patterns  at  only  minor 
changes.  In  this  case,  the  cost  function  (4)  must  include  a  new  term  to  take 
into  account  the  ripple  level  of  the  shaped  region.  To  obtain  ootimal  results, 
each  ripple  peak  must  be  individually  controlled  as  well  as  the  peak  sidelobe 
level.  In  our  example,  we  have  used  a  16  x  16  element  oianar  array,  with 
rectangular  grid,  circular  boundary,  and  a  A/2  interelement  spacing.  Tne 
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Figure  6.  Svni;hesized  sum  power  pattern  with  three  element  failures  using 
phase-only  control.  Sidelobes  lowered  in  the  whole  ©-range  (0°-360°) . 


Figure  7.  S^'mthesized  sum  power  pattern  with  three 
phase-oniy  control.  Sidelobes  lowered  in  the  ©-range 
increased  in  the  ©-range  (lS0°-360°) . 


element  failures  using 
(0°-iS0°)  but  sHghtiy 
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Figure  8.  Power  pactera  obtained  by  sampling  a  ©-symmetric  pure  real 
aperture  distribution  corresponding  to  a  ©-symmetric  dat-topped  beam 
(SLL=  —25  dB,  ripple  ±0.5  dB). 
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initial  pattern  was  obtained  by  sampling  a  d>  -symmetric  pure  real  aperture 
distribution  corresponding  to  a  p -symmetric  fiat-topped  beam  with  a  peak 
sidelobe  level  of  —25  dB.  ripple  r:  0.5  dB  [11],  and  is  shown  in  Figure  8.  We 
have  simulated  two  nuU  elemenrs  located  in  the  first  quadrant  as  Figure  9 
show's.  The  pattern  obtained  after  nulling  these  two  elements  has  a  peak  a 
sidelobe  level  of  —18.7  dB,  ripple  il.l  dB  and  is  plotted  in  Figure  10. 
Using  the  amplitude-only  synthesis  we  could  reduce  the  sideiobes  to  —24.5 
dB  and  the  ripple  level  to  ±0.5  dB.  The  synthesized  pattern  is  shown  in 
Figure  11.  Again,  we  have  found  that  the  amplitude  only  synthesis  allows 
to  reduce  the  sidelobe  level  as  well  as  the  ripple  level  down  to  near  the 
design  level.  Regarding  to  phase-only  synthesis,  we  could  not  obtain  so 
good  results  as  modifying  the  amplitude  weights  because  a  reduction  of  the 
sideiobes  yielded  an  increasing  of  the  ripple  level. 

CONCLUSIONS 

A  method  of  array  failure  correction  for  planar  arrays  has  been  orooosed. 
The  optimization  performed  to  the  degraded  pattern  does  not  require  to 
modify  both  the  amplitude  and  phase  weights  of  the  non-defective  elements, 
but  only  the  amplitude  or  the  phase  aperture  distribution.  This  reduces  the 
hardware  complexity  required,  specially  in  the  case  of  the  phase  synthesis 
which  has  found  to  yield  optimal  results  for  asymmetric  sideiobes.  This 
method  is  applicable  to  planar  arrays  of  a  wide  variety  of  grid  structures,  to 
non-uniformly  spaced  arrays,  and  even  to  antennas  made  up  of  non-identical 
elements. 
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Abstract 

The  authors  report  the  application  of  the  simulated  annealing  optimization  technique  to 
the  synthesis  of  linear  array  sum  patterns  by  highly  efficient  real  aperture  distributions  with  no 
edge  brightening. 

1.  Introduction 

Villeneuve  [1]  published  a  pattern  synthesis  method  that  is  the  equivalent,  for  linear 
arrays,  of  Taylor's  method  for  synthesizing  sum  patterns  of  near-optimal  peak  directivity  for 
linear  antenna  apertures  [2].  Like  Taylor's  method,  Villeneuve  s  affords  excitation  distributions 
with  considerable  edge  brightening  (i.e.  a  pronounced  rise  in  amplitude  at  the  ends  of  the  array). 
To  alleviate  edge  brightening  (at  the  expense  of  efiiciency),  McNamara  [3]  modified  Villeneuve's 
method  by  allowing  the  outer  side  lobe  level  to  fall  off  more  rapidly.  Edge  brightening  can  also 
be  reduced  by  null  filling,  again  at  the  expense  of  efficiency  (and  with  the  added  complication 
of  making  the  excitation  distribution  complex  instead  of  real)  [4]. 

We  recently  reported  the  use  of  a  simulated  annealing  optimization  method  to  synthesize 
radiation  patterns  approximating  linear  and  circular  Taylor  sum  pattern  by  means  of  highly 
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efficient  real  aperture  distributions  with  smoothly  varying  amplitudes  and  no  edge  brightening 
[5].  We  have  now  adapted  this  procedure  to  linear  arrays  producing  patterns  approximating 
Villeneuve  sum  patterns. 


2.  Method 

Given  the  W  roots  b;  of  a  Villeneuve  pattern  [1],  small  perturbations  66,.  are  optimized  by 
the  simulated  annealing  method  [6]  using  the  cost  function 


(1) 


where  |4ax/4.ini  dynamic  range  ratio  of  the  aperture  distribution,  is  a  measure  of 

excitation  distribution  smoothness  controlling  edge  brightening  (V=  max{i2^.},  where  Rj  is  the 
ratio  between  the y-th  excitation  peak  and  the  lesser  of  the  flanking  minima),  rj  is  the  efficiency 
of  the  array  (the  ratio  between  its  peak  directivity  and  that  of  a  uniformly  excited  array  with  the 
same  number  of  elements),  the  f-  (which  penalize  side  lobe  levels  SLL,.  in  excess  of  the  desired 


levels  SLL,-  ,,)  are  defined  by 


(5ZL, -  SLL.J)^  if  SLL, i  SLL,^^ 

0  otherwise 


(2) 


and  the  are  weights  that  can  be  varied  by  the  designer  in  accordance  with  his  or  her  design 
priorities. 


3.  Results 

The  dashed  curve  in  Fig.  1  shows  a  -20  dB  Villeneuve  pattern  calculated  for  an  array  of 
30  isotropic  elements  M2  apart  (n  =  6,  the  value  optimizing  efficiency  for  W  =  30  and 
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SLL  =  -20  dB).  Fig.2  (dashed  curve)  shows  the  calculated  excitation  distribution,  which  exhibits 
marked  edge  brightening.  The  solid  curves  in  Figs.l  and  2  show  the  radiation  pattern  and 
excitation  distribution  calculated  by  the  method  described  in  this  letter.  The  new  pattern  respects 
the  SLL  specification  (outer  side  lobes  are  in  fact  lower  than  in  the  Villeneuve  pattern),  and  the 
excitation  distribution  shows  no  edge  brightening.  The  excitation  distribution  also  has  a  smaller 
dynamic  range  ratio  and  smaller  than  the  Villeneuve  distribution,  and  is  even 

marginally  more  efficient  (Table  1). 


4.  Conclusions 

The  method  described  in  this  letter  allows  edge  brightening  to  be  eliminated  from 

% 

Villeneuve  excitation  distributions  with  no  increase  in  SLL  or  loss  of  efficiency  while 
minimizing  dynamic  range  ratio  and  It  can  also  be  applied  with  advantage  to  the 

synthesis  difference  patterns,  shaped  beams  and  sum  patterns  with  inner  side  lobes  that  are 
depressed  relative  to  outer  side  lobes. 
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LEGENDS  FOR  THE  TABLE  AND  FIGURES 

Table  1.  Dynamic  range  ratio,  adjacent  excitation  ratio  and  efficiency  of  the  patterns  of  Fig.l. 

Fig.l.  Power  plots  of  a  Villeneuve  pattern  (SLL  =  -20  dB,  =  6;  dashed  curve)  and  the 
pattern  achieved  in  this  work  taking  the  Villeneuve  pattern  as  its  starting  point  (solid 
curve). 

Fi<^.2.  Excitation  distributions  of  the  patterns  of  Fig.l  (dashed  curve,  Villeneuve  pattern;  solid 

curve,  this  work). 
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.K  t  The  authors  report  the  application  of  the  simulated  annealing  optimization  technique  to  the 
^  "" 

edge  brightening. 


1.  Introduccidn 

En  la  literatura  previa,  Taylor  [1]  introdujo  una 
tecnica,  aplicable  a  distribuciones  lineales,  que 
permitia  generar  diagramas  de  tipo  suma  de  alta 
eficiencia  con  un  nivel  de  lobulos  laterales  deseado. 
Esta  tecnica  ha  sido  muy  utilizada  en  la  detemunacion 
de  las  excitaciones  de  arrays  lineales,  mediante  el 
muestreo  de  la  distribucion  de  abertura  contmua  en 
aquellas  posiciones  en  las  que  se  encuentran 
localizados  los  elementos  del  array.  El  mayor 
problema  que  presentan  estas  distribuciones  de  Taylor 
es  que  sus  amplitudes  de  excitacion,  en  muchos  casos, 
varian  muy  rapidamente  en  los  extremes  de  la 
abertura,  efecto  que  se  conoce  como  ‘  edge 
brigthening”  y  que  puede  hacer_  mviable  la 
implenientacion  fisica  de  la  antena.  Recienteraente,  se 
ha  aplicado  la  tecnica  de  “simulated  annealing”  para 
conseguir  una  topografia  6ptima  de  lobulos  laterales 
de  modo  que  se  eliminase  el  “edge  brightening 
manteniendo  una  alta  eficiencia  en  la  abertura  [2]. 

Para  arrays  discretes,  Villeneuve  introdujo  una 
tecnica  que  permitia  sintetizar  directamente  diagramas 
muy  similares  a  los  de  Taylor  [3].  Aunque  las 
distribuciones  de  amplitud  resultantes  son  ^  muy 
eficientes,  siguen  presentando  el  efecto  de  “edge 
brightening”.  Posteriormente,  McNamara  [4J 
introdujo  una  gcneralizacion  de  esta  tecnica  que, 
introduciendo  un  nuevo  parametro  relacionado  con  la 
velocidad  de  caida  de  los  lobulos  exteriores,  permitia 
reducir  el  “edge  brightening”  a  costa  de  disminuir  la 
eficiencia  de  la  abertura.  En  un  trabajo  posterior,  se 
demostro  que  tambien  era  posible  aliviar  cste  efecto 
mediante  el  relleno  de  los  ceros  de  los  diagramas  de 
radiacion  [5].  Sin  embargo,  esta  tecnica  conlleva  el 
uso  de  una  distribucion  de  abertura  compleja  y  lleva 
asociada  una  perdida  en  la  eficiencia  de  la  abertura. 

Este  trabajo,  que  es  una  extension  de  [2]  a  arrays 
discretos,  permite  sintetizar  diagramas  de  tipo  suma 
mediante  distribuciones  de  abertura  reales  que  no 
presentan  “edge  brightening”,  manteniendo  una  alta 
eficiencia.  La  tecnica  perturba  ligeramente  las  raices 
de  las  distribuciones  de  Villeneuve  de  forma  que  se 
minimice  una  funcion  de  costo  que  tenga  en  cuenta  el 
nivel  de  lobulos  laterales,  la  variabilidad  de  las 


excitaciones  y  todos  los  demas  parametros  de  interes 
eneldiseho. 


2.  Dcscripcion  del  metodo 

El  proceso  parte  de  una  distribucion  inicial  de 
raices  dadas  por  (como  por  ejemplo,  las 
correspondientes  a  una  distribucion  de  Villeneuve  [j]) 
y  las  va  perturbando  ligeramente  de  forma  iterativa. 
Dichas  perturbaciones  5b ^  se  calculan  mediante  la 
tecnica  de  “simulated  annealing”  a  traves  de  la 
minimizacion  de  una  funcion  de  costo  dada  por  la 
siguiente  expresion; 

C(66j,66j . 

/ 

en  donde/j,  que  penaliza  el  incremento  de  los  lobulos 
laterales  por  encima  del  nivel  de  diseno,  se  define 
como 

{SLL,, -  SLL,/  si  SLL,^  ^  SLL,^  ^2) 
0  en  otro  caso 


en  donde  y  SLL,,,  es  el  nivel  obtenido  y  deseado 
del  lobulo  i-dsimo  respectivamente; 
rango  dindmico  de  la  distribucion  de  abertura,  V,  una 
medida  de  la  suavidad  de  dicha  distribucion,  controla 
el  “edge  brightening”  {V=max{Rj} ,  en  donde  R.  es  la 
razon  entre  el  maximo  j  y  el  menor  de  los  minimos 
adyacentes  de  la  abertura);  77  es  la  eficiencia  de  la 
abertura,  definida  como  la  razon  entre  la  directividad 
pico  del  diagrama  obtenido  y  la  directividad  pico  de 
un  array  uniforme  con  el  mismo  niimero  de 
elementos.  Finalmente  c,  son  los  factores  de  peso  que 
pueden  ser  ajustados  en  funcion  de  las 
especificaciones  del  diseno. 

S.Resultados 

Para  ilustrar  el  metodo  se  ha  utilizado  un  array 
lineal  de  30  elementos  isotropicos  equiespaciados  272. 
Se  ha  partido  de  una  distribucion  de  Villeneuve  con 
un  nivel  de  lobulos  laterales  de  -20  dB  y  un  indice  de 
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transicion  n  =6  (se  ha  elegido  este  valor  puesto  que 
optimiza  la  eficiencia  de  la  distribucion  para  el  SLL  y 
niimero  de  elementos  considerado).  El  diagrama 
resultante  se  muestra  en  la  fig.  1,  mientras  que  en  la 
fig.  2  se  muestra  la  distribucion  de  abertura 
correspondiente,  la  cual  presenta  “edge  brightening”. 
El  diagrama  obtenido  tras  la  optimizacion  y  su 
distribucion  de  abertura  se  representan  en  trazo 
continue  en  dichas  figuras.  Tal  y  como  se  observa,  se 
ha  encontrado  la  topografia  optima  de  lobulos 
laterales  (manteniendo  el  nivel  maximo  de  diseno  de 
-20  dB)  que  permite  eliminar  el  “edge  brightening”  de 
la  distribucion  de  abertura.  En  la  tabla  I  se  muestran 
los  parametros  de  ambos  diagramas.  Se  observa  que  la 
distribucion  de  abertura  obtenida  presenta  un  menor 
range  dinamico  asi  como  una  menor  variabilidad  de 
las  excitaciones  entre  elementos  colindantes  que  la 
inicial,  manteniendo  una  elevada  eficiencia. 


n 

i 

j  Villeneuve 

1.75 

1.13 

0.966 

1 

j  Optimizado 

1.55 

1.09 

_ - 

0.967 

Tabla  I.  Comparacion  del  range  dinamico,  variabilidad  de  las 
c.xcitacioncs  y  eficiencia  del  diagrama  inicial  y  el  obtenido  Iras  la 
optimizacion. 


4.  Conclusiones 

Se  ha  introducido  una  tecnica  que  permite 
eliminar  el  efecto  de  “edge  brightening”  en  las 
distribuciones  de  abertura  correspondientes  a  los 
diagramas  de  radiacion  de  tipo  Villeneuve.  La 
solucion  obtenida  tambien  optimiza  el  rango  dinamico 


Fig.  1.  Diagrama  de  radiacion  correspondiente  a  la  distribucion  dc 
Villeneuve  (SLL— 20,  n  =6)  y  el  obtenido  tras  la  optimizacion. 


y  la  variabilidad  de  las  excitaciones,  ademas  de 
mantener  una  alta  eficiencia.  El  metodo  tambien 
permite  eliminar  el  “edge  brightening”  en  patrones 
con  lobulos  deprimidos,  siendo  tambien  aplicable  a 
diagramas  diferencia  y  de  haz  perfilado. 
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metaUic  sides.  The  element  is  fed  by  means  of  coa.xial  probes  at 
the  centre  of  the  lower  edge  of  one  or  more  of  the  vertical  sides.  A 
single  feed  is  sufficient  for  linear  polarisation  and  feeds  excited  m 
phase  quadrature  on  two  adjacent  sides  may  be  used  for  circular 
polarisation.  The  circularly  polarised  (CP)  element  in  the  Figure 
has  four-quadrant  symmetry;  however,  other  CP  implementations 
which  have  only  two  vertical  faces  have  been  tested,  and  linearly 
polarised  elements  with  a  single  vertical  face  have  also  been  evalu¬ 
ated.  In  these  other  cases  the  element  pattern  is  less  symmetncal, 
with  the  beam  peak  occurring  off  broadside.  This  beam  squint  has 
been  found  to  be  useful  in  some  array  applications  [1]. 


a  “  [gw 


Fig.  1  Radiator  geometry 

a  Top  view 
b  Side  view 


Fig.  2  Measured  and  predicted  RHCP  radiation  patterris  at  1.65  GH- 


- measured 

- predicted 


Fig.  3  Measured  return  loss 


Experimental  and  analytical  results.'  A  circularly  polarised  version 
of  the  new  radiator  has  been  analysed  using  the  WIPL  software 
package  [3],  which  implements  a  method  of  moments  solution  to 
the  integral  equation  problem.  In  this  case,  the  element  has  two 
vertical  faces  with  one  feed  probe  exciting  each  face.  The  probes 
are  fed  with  equal  current  magnitudes  and  relative  phases  of  90°. 
The  element,  which  has  width  w  =  3.8cm  and  height  h  =  2.54cm, 
has  been  positioned  at  the  centre  of  a  22cm  square  ground  plane. 
The  predicted  RHCP  radiation  pattern  at  1.65  GHz  and  the  corre¬ 
sponding  measured  pattern  are  given  in  Fig.  2.  The  protot^qje  ele¬ 
ment  was  excited,  by  a  rmcrostrip  circuit  incorporating  a  Wilkinson 
hybrid,  followed  by  transmission  lines  differing  in  phase  length  by 
90°  and  microstrip  transformers  matching  to  the  2400  element 


resonant  impedance.  The  return  loss,  which  is  plotted  in  Fig.  3, 
was  measured  without  the  hybrid,  at  the  input  to  one  of  the 
matching  transformers.  The  hybrid  absorbs  reflections  from  the 
element  Input  ports,  and  consequently  the  return  loss  at  the  hybrid 
input  port  is  not  a  good  measure  of  the  element  performance. 

Conclusion:  The  compact  radiator  described  here  has  been  used  in 
a  number  of  successful  antenna  development  programs  [1.  2]  and 
has  been  found  to  be  a  useful  alternative  to  microstrlp  patches  and 
other  low  profile  elements.  The  radiator  has  a  footprint  ot  only 

O. 04  square  wavelengths  while  achieving  a  bandvricta  of  >  17% 
with  6dB  return  loss'and  a  -6dB  CP  beamwidth  of  >  120°. 
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Design  of  efficient,  easily  feed-matched 
array  antennas  by  joint  optimisation  of 
excitations  and  element  geometry;  pencil 
beam  example 

J,A.  Rodriguez,  F.  Ares  and  G.  Francescheiti 

By  using  simulated  annealing  for  simultaneous  optimisation  of  the 
geometries  and  excitations  of  the  array  elements,  it  is  often 
possible  to  design  array  antennas  that  generate  specified  radiation 
power  patterns,  have  elements  with  real  active  impedances,  and 
are  easily  feed-matched.  The  example  of  a  twelve^lement  pencfl 
beam  antenna  fed  by  coaxial  cables  of  practical  dimensions  with 
no  additional  circuitry  is  described. 


Introduction:  The  efficiency,  and  indeed  the  realisability,  of  a  given 
array  antenna  distribution  of  excitations  depends  on  the  mutual 
coupling  among  the  array  elements.  Most  available  methods  of 
array  antenna  power  pattern  synthesis  make  no  attempt  to  take 
mutual  coupling  explicitly  into  account  during  the  calculation  of 
an  excitation  distribution  attaining  the  desired  power  pattern. 
Coupling  is  minimised  by  reducing  the  dynamic  range  of  the 
excitations  or  the  ratio  between  the  e.xcitations  of  adjacent  ele¬ 
ments  [1  -  4].  However,  to  design  an  efficiently  fed  antenna  that 
also  radiates  the  desired  pattern  efficiently,  it  is  not  sufficient  to 
ensure  that  the  excitation  distribution  has  small  dynamic  range; 
even  excitation  distributions  with  quite  small  dynamic  ranges  can 
be  very  costly  or  impossible  to  implement  efficiently  because  of 
coupling  problems  [5,  9],  This  is  because  the  active  impedances  of 
the  array  elements  and  the  impedances  presented  to  the  main  feed 
line  by  the  branch  lines  to  the  array  elements  depend  on  both  the 
excitation  distribution  and  the  geometry  of  the  elements  and/or 
array.  In  general,  both  excitations  and  geometry  must  be  varied  to 
find  a  combination  which  achieves  the  desired  radiation  pattern 
and  antenna  properties. 

In  the  past,  it  has  been  common  to  follow  a  three-stage  proce¬ 
dure.  Pint,  under  the  assumption  that  all  the  radiating  elements 
are  identical,  and  with  fixed  geometry,  an  excitation  distribution 
attaining  the  desired  radiation  pattern  is  obtained.  Next,  with  the 
excitation  distribution  fixed,  the  element  geometry  is  calculated  to 
ensure  efficient  radiation  (this  results  in  pattern  degradation,  but 
the  degradation  is  often  negligible  except  for  extreme  angles). 
Finally"  feed  matching  is  ensured  by  circulating  feed  network 
parameters. 
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In  this  Letter  we  describe  how  an  efficient,  easily  feed-matched 
array  producing  a  specified  radiation  power  pattern  can  be  synthe¬ 
sised  by  optimising  element  geometry  and  in-phase  excitations 
simultaneously.  We  take  as  our  example  a  linear  array  of  twelve 
centre-fed  dipoles  that  generates  a  pencil  beam  and  is  fed  by  coax¬ 
ial  cables  with  no  additional  circuitry;  the  method  is  readily  gener- 
alisable  to  planar  arrays  and  to  arrays  with  other  types  of 
radiating  elements  (e.g.  microstrip  dipoles  and  waveguide  slots,  the 
design  of  which  raises  the  same  difficulty  [6  -  9]). 


Theory:  Consider  a  linear  array  of  N  centre-fed  dipoles  located  a 
distance  h  from  a  ground  plane  and  connected  to  the  main  coaxial 
feed  line  by  coaxial  branch  lines.  For  practical  values  of  /i  {/i  <  TJ 
4,  where  X  is  the  radiation  wavelength),  the  image  principle  is 
applicable  even  for  a  finite  ground  plane,  which  rnay  Aus  be 
replaced,  for  the  purposes  of  analysis,  by  a  parallel,  identical  vir¬ 
tual  array  excited  180°  out  of  phase  with  respect  to  the  real  array. 
The  active  impedance  of  the  ith  dipole,  Zf,  is  therefore  given  by 

(1) 

j=i 

where  /„  is  the  current  in  the  ath  dipole  and  Zy  is  the  mutual 
impedance  between  dipoles  i  and  j  minus  the  mutual  impedance 
between  dipole  i  and  the  virtual  image  of  dipole  j  (_Zij  is  the  self¬ 
impedance  of  dipole  i  minus  its  mutual  impedance  with  its  image) 
[5].  Z,j  depends  on  the  geometry  of  the  dipoles  and  their  mutual 
geometric  relations. 

If  the  array  element  excitations  are  all  in  phase,  then  for  effi¬ 
cient  radiation  the  active  impedances  Z,"  should  aU  be  real  (and 
are  hereafter  denoted  by  Rf).  This  additional  constraint  requires 
the  simultaneous  optimisation  of  the  currents  I,  and  the  array  and/ 
or  dipole  geometry.  This  can  be  done  by  using  simulated  anneal¬ 
ing  [10]  to  minimise  a  cost  function  including  terms  for  deviation 
from  the  desired  pattern  and  for  the  ima^nary  parts  of  the  active 
impedances.  In  addition,  easy  feed  matching  can  be  ensured  if  the 
cost  function  also  includes  appropriate  terms  for  feed  network 
parameters.  Note  that  the  variation  of  element  geometry  means^ 
that  the  radiation  pattern  cannot  be  factorised  as  the  product  of 
an  array  factor  and  a  single  element  factor,  but  must  be  calculated 
from  the  individual  patterns  of  each  dipole. 


Example:  Consider  the  synthesis  of  a  pencil  beam  pattern  by  a  lin¬ 
ear  array  of  twelve  centre-fed  dipoles  of  radius  0.0048  X  laid  out  in 
the  x-y  plane  with  a  distance  'KIT.  from  each  other  with  their  cen¬ 
tres  on  the  y  axis.  The  array  centre  is  at  the  origin,  and  a  ground 
plane  lies  at  z  =  —h  —  — X/4.  In  the  y-z  plane  ((j)  —  7t/2),  the  radiation 
pattern  is  to  have  a  maximum  sidelobe  level  of-18dB  (this  max¬ 
imises  the  peak  directivity  of  the  twelve-element  Dolph-Chebyshev 
solution  we  assume  as  our  starting  point).  The  dipoles  are  to  be 
fed  in  phase  from  a  single  10£2  main  transmission  line  via  quarter- 
wave  coaxial  cables  with  no  extra  circuitiy.  We  assume  that  the 
radiation  patterns  of  the  dipoles  are  given  in  terms  of  their  lengths 


2/,  by 


Fi{9,(p)  —  sin{khcos6) 


cos(fcli  sin  9  cos  o)  —  cos{kli) 
y  1  —  sin"  9  cos^  p 


(2) 


where  k  -  2rJX,  and  that  their  self  and  mutual  impedances  are  as 
civen  in  [5]  (this  implies  that  dipole  lengths  can  be  varied  only  in 
the  interval  [1.3/A:,  1.7/A:]  in  which  the  assumed  expression  for  self¬ 
impedance  is  valid). 

To  verify  pattern  specification  when  optimising  element  currents 
and/or  geometry,  the  cost  function  to  be  minimised  can  include 
the  terms  c.ZXA^,)^  and  cJD,,  where  &S,  is  the  difference  between 
the  calculated  level  of  the  side  lobe  s  and  the  maximum  tolerated 
level,  if  this  difference  is  positive,  and  AS,  =  0  otherwise,  and  D,  is 
the  peak  directivity  of  the  pattern  (c,  and  c,  are  appropriate 
weights).  To  ensure  that  the  array  elements  have  real  active  imped¬ 
ances,  the  cost  function  can  include  the  term  CiLiiXff ,  where  A)"  's 
the  imaginary  part  of  Z"  (again,  is  an  appropriate  weight). 
Finally, '^the  feed  specification  means  that  the  feed  network  is 
determined  by  the  characteristic  resistances  /?,*'■  of  the  quarter- 
wave  branch  lines.  We  have  =  (Z“-Ri)’'’.  where  Z  are  the 
resistances  presented  to  the  mam  feed  line  by  the  branch  lines  and 
are  calculated  by  solving  the  equations  [5] 

IUIRj=:{ljR])lilfRt)  (3) 
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together  with  the  teed  matching  condition 

^(1/i?,)  =  l/R^^  (4) 

i 

To  meet  the  requirement  that  the  7?,®^  have  realistic  values,  e.g. 
<  140D,  the  cost  function  can  include  the  term  CjZ^g,-,  where  g,  = 
p_BL  _140,  if  this  is  positive,  and  g,  =  0  otherwise,  with  c,  an 
appropriate  weight.  If  1,  are  all  given  the  fixed  value  A/4,  the  radia¬ 
tion  pattern  specification  is  satisfied  by  the  Dolph-Chebyshev 
solution  [5]  listed  in  Table  1.  This  solution  is  inefficient  because  of 
the  coupling  indicated  by  the  complex  active  impedances  of  the 
dipoles  (also  listed  in  Table  1). 

Table  1:  Excitations,  lengths  and  active  impedances  of  elements  of 
12-dipole  array  when  dipole  length  is  kept  uniform  (=  A/ 
2)  and  excitations  are  optimised  to  generate  Dolph- 
Chebyshev  pattern  with  -18dB  side  lobes 


/ 

I, 

2/A 

r  0 
“  t 

1;  12 

1.000 

0.500 

n 

89.7  55.4; 

2;  11 

0.665 

0.500 

110.1-4.1; 

3;  10 

0.826 

0.500 

91.6  +  28.1; 

4;  9 

0.964 

0.500- 

101.2  20.8; 

5;  8 

1.064 

0.500 

97.0  -^23.2; 

6;  7 

1.117 

0.500 

98.6  22.5; 

Peak  directivity  =  57.9 


Table  2:  Excitations,  lengths  and  active  impedances  of  elements  of 
12-dipole  array  after  perturbation  of  dipole  lengths  of  Table  1  to 
afford  real  active  impedances  and  best  possible  set  of  branch  line 
characteristic  impedances,  keeping  excitations  fixed 


/ 

1 

2/A 

•ya 

i 

Rf’- 

£1 

a 

1;  12 

1.000 

0.461 

70.9 

96.6 

2;  11 

0.665 

0.498 

107.0 

145.2 

3;  10 

0.826 

0.479 

81.7 

117.0 

4;  9 

0.964 

0.482 

90.4 

100.3 

1.064 

0.481 

86.9 

90.8 

mm 

1.117 

0.481 

88.2 

86.4 

Peak  directivity  =  57.5,  computation  time  =  1.2  min 
Computation  time  refers  to  calculations  on  a  Pentium  II  running 


at  266MF1Z 


Table  3:  As  for  Table  2,  but  with  simultaneous  optimisation  of 
dipole  lengths  and  excitations 


i 

h 

2/A 

^  f 

RJ’ 

£1 

n 

1;  12 

1.403 

0.466 

74.7 

119.1 

2;  11 

1.426 

0.485 

95.4 

117.2 

3;  10 

1.456 

0.479 

82.7 

114.8 

4;  9 

1.444 

0.490 

95.7 

115.8 

5;  8 

2.092 

0.475 

83.1 

79.9 

6;  7 

1.875 

0.484 

90.1 

89.1 

Peak  directivity  =  57.0,  computation  time  =  2.1  min 


.4  solution  with  the  same  current  set  that  radiates  efficiently  and 
maintaim  the  directivity  can  be  found  by  perturbing  /.  with  the 
currents  fixed.  The  cost  function  C  =  c,Z,(AS,)^  -r  cJD^  -t- 
is  used,  to  move  to  zero  the  imaginary  part  of  the  active  imped¬ 
ances  while  penalising  directivity  loss.  However,  if  simultaneous 
satisfaction  of  the  feed  specification  is  attempted  by  using  the  cost 
function  C  =  c,S,(A5,)^  +  cJD^  +  c^{Xff  -t-  +c,Z^,  with  g,  =  Rf’- 
-140,  no  solution  is  attainable.  To  obtain  the  solution  given  in 
Table  2,  it  is  necessary  to  raise  the  limiting  values  of  Rf’-  to  145  Q 
(i.e.  g,  =  Rf’-  -145),  which  provides  characteristic  impedances 
associated  with  coaxial  cables  which  are  difficult  to  manufacture. 

In  contrast,  if  the  lengths  21,  and  currents  I,  are  simultaneously 
\-aried  by  using  the  same  cost  function  (with  g,  =  R,^’-  -140),  then 
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the  desired  power  pattern,  real  active  impedances  and  more  rea¬ 
sonable  Ri^'-  values  can  all  be  achieved  together  (Table  3). 
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Two-dimensional  beam-scanning  phase- 
shifterless  technique  using  linear  active 
leaky-wave  antenna  array 

Cheng-Chi  Hu,  C.F.  Jou  and  Jin-Jei  Wu 

.A  novel  two-dimensional  electronic  beam  scanning  technique 
using  a  linear  leaky-wave  antenna  array  with  coupled  oscillators  is 
introduced,  eliminating  the  need  for  phase  shifters.  The  measured 
H-plane  main  beam  can  be  continuously  scanned  from  70  to  40° 
as  the  frequency  varies  from  7.9  to  9.05  GHz.  By  detuning  the  free 
running  frequencies  of  the  end  elements,  the  measured  E-plane 
main  beam  can  be  continuously  scanned  from  -22  to  -26°. 

Introduction:  The  comple.xity  usually  associated  with  the  two- 
dimensional  (2D)  scanning  array  offers  a  special  challenge  for  the 
planar  active  phase  array  design.  In  1990,  Oliner  [1]  proposed  a 
2D  scanning  array  using  a  one-dimensional  (ID)  phased  array  of 
leaky-wave  line-source  antennas.  A  pencil  beam  can  scan  in  both 
elevation  and  azimuth  planes,  however,  phase  shifters  are  required 
in  this  design.  These  phase  shifters  usually  require  complicated 
control  circuitry  which  it  may  also  be  difficult  to  achieve  in  the 
HmltpH  space.  In  1994,  Liao  and  York  [2]  proposed  a  new  phase- 
shifterless  ID  beam-scanning  technique  using  a  patch  antenna 
array  with  coupled  oscillators.  By  controlling  the  free  running  fre¬ 
quencies  of  the  end  elements  of  the  array,  the  main  beam  can  scan 
in  the  azimuth  plane.  In  this  Letter,  we  extend  the  work  to  encom¬ 


pass  the  phase  control  technique  [2]  and  leaky-wave  antenna  char¬ 
acteristic  [1,  3,  4],  leading  to  a  novel  method  for  phase-shifterless 
2D  electronic  beam  scanning. 

microstrip 


tea/ll 

Fig.  1  Configuration  of  active  microstrip  leaky-wave  antenna  array 
11’  =  12mm,  L  -  100mm,  rf  =  1  X, 


Fig.  2  Geometry  and  coordinate  system  for  microstrip  leaky-wave 
antenna 


Design  and  measurement  results:  Fig.  1  shows  the  microstrip  reali¬ 
sation  of  a  two-element  phaise-shifterless  active  leaky-wave 
antenna  array  structure.  The  varactor-tuned  oscillator  is  based  on 
that  we  previously  designed  in  [4].  Coupling  circuits  are  designed 
to  provide  in-phase  coupling,  which  ensures  a  stable  in-phase 
mode  of  operation.  The  elements  within  the  array  are  coupled  to 
one  another  using  one  wavelength  long  transmission  lines.  To 
excite  the  first  higher  order  mode,  the  microstrip  leaky-wave 
antenna  is  fed  asymmetrically  [5],  The  circuit  is  designed  and  fab¬ 
ricated  on  an  RT/Duroid  substrate  with  a  dielectric  constant  of 
2.2  and  a  thickness  of  20mils.  An  NEC  NE42484  low  noise 
HEMT  is  used,  and  the  drain  is  biased  at  2.0V  with  a  drain  cur¬ 
rent  of  10mA.  The  GaAs  varactor  (M/A-COM  MA46410)  is  used 
as  a  tuning  varactor,  which  has  a  capacitance  ratio  of  10:1  and  a 
capacitance  of  0.5pF  at  4V. 

Elevation  plane  scanning:  To  understand  the  radiation  properties 
of  such  a  microstrip  leaky- wave  antenna  [1],  we  obtained  its  com¬ 
plex  propagation  constants  |3  -  ja  of  the  first  higher  microstrip 
mode  in  its  leakv'  range,  where  p  is  the  phase  constant,  and  a  is 
the  attenuation  constant.  Such  complex  propagation  constants 
represent  a  forward  leaky-wave  radiating  into  the  space  at  an 
angle  6„  =  cor'  (p/X'o),  where  0„  is  the  angle  of  the  Ixam  ma.xi- 
mum  measured  from  the  z-a.xis,  and  A'o  is  the  free-space  wave 
number.  In  addition,  the  scanning  angle  0,„  can  be  varied  with  fre¬ 
quency. 

For  a  tuning  voltage  of  1.0  to  10 V,  our  active  leaky-wave 
antenna  array  exhibits  a  tuning  bandwidth  of  7.9  to  9.05GHz  cor¬ 
responding  to  a  measured  main  beam  position  from  70  to  40°.  Fig. 
3  shows  the  experimental  results  of  the  radiating  patterns  scaned 
in  the  elevation  plane.  The  maximum  effective  radiated  power 
(ERP)  of  this  active  antenna  array  is  ~  20dBm  ±  2dBm  through¬ 
out  the  frequency  tuning  range.  The  difference  m  the  power  level 
of  the  main  beam  is  caused  mainly  by  the  varied  impedance  of  the 
microsirip  leaky-wave  antenna  as  function  of  frequency. 
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Abstract 

By  using  simulated  annealing  for  simultaneous  optimization  of  the  geometries  and  excitations  of  the  array 
elements,  it  is  often  possible  to  design  array  antennas  that  generate  specified  radiation  power  patterns,  have  elements 
with  real  active  impedances,  and  are  easily  feed-matched.  The  example  of  a  twelve-element  pencil  beam  antenna  fed 
by  coaxial  cables  of  practical  dimensions  with  no  additional  circuitry  is  described. 


1.  Introduction 

The  efficiency,  and  indeed  the  realizability,  of  a  given  array  antenna  distribution  of  excitations  depends  on  the 
mutual  coupling  among  the  array  elements.  Most  available  methods  of  array  antenna  power  pattern  synthesis  make  no 
attempt  to  take  mutual  coupling  explicitly  into  account  during  calculation  of  an  excitation  distribution  attaininz  the 
desired  power  pattern.  Coupling  is  minimized  by  reducing  the  dynamic  range  of  the  excitations  or  the  ratio  between  the 
excitations  of  adjacent  elements  [1-4].  However,  to  design  an  efficiently  fed  antenna  that  also  radiates  the  desired  pattern 
efficiently,  it  is  not  sufEcient  to  ensure  that  the  excitation  distribution  has  small  dynamic  rangej  even  excitation 
■distributions  with  quite  small  dynamic  ranges  can  be  very  costly  or  impossible  to  implement  efficiently  because  of 
coupling  problems  [5,9].  This  is  because  the  active  impedances  of  the  array  elements  and  the  impedances  presented  to 
the  main  feed  line  by  the  branch  lines  to  the  array  elements  depend  on  both  the  excitation  distribution  and  the  geometry 
of  the  elements  and/or  array.  In  general,  both  excitations  and  geometry  must  be  varied  in  order  to  find  a  combination 
attaining  the  desired  radiation  pattern  and  antenna  properties. 

In  the  past,  it  has  been  common  to  follow  a  three-stage  procedure.  First,  under  the  assumption  that  all  the  radiating 
elements  are  identical,  and  with  fixed  geometry,  an  excitation  distribution  attaining  the  desired  radiation  pattern  is 
obtained.  Next,  with  the  excitation  distribution  fixed,  the  element  geometry  is  calculated  to  ensure  efficient  radiation 
(this  results  in  partem  degradation,  but  the  degradation  is  often  negligible  except  for  extreme  angles).  Finally,  feed 
matching  is  ensured  by  calculating  feed  network  parameters. 

In  this  commumcation  we  describe  how  an  efficient,  easily  feed-matched  array  producing  a  specified  radiation 
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power  pattern  can  be  synthesized  by  optimizing  element  geometry  and  in-phase  excitations  simultaneously.  We  take 
as  our  example  a  linear  array  of  hvelve  centre-fed  dipoles  that  generates  a  pencil  beam  and  is  fed  by  coaxial  cables  with 
no  additional  circuitry,  the  method  is  readily  generalizable  to  planar  arrays  and  to  arrays  with  other  types  oi  radiatmg 
elements  (e.g.,  microstrip  dipoles  and  waveguide  slots,  whose  design  meets  the  same  difficulty  [6-9]). 


2.  Theory 

Consider  a  linear  array  of  centre-fed  dipoles  located  a  distance  h  from  a  ground  plane  and  connected  to  the  mam 
coaxial  feed  line  by  coaxial  branch  lines.  For  practical  values  ofh  (h  ^  M,  where  X  is  the  radiation  wavelength),  the 
image  principle  is  applicable  even  for  a  finite  ground  plane,  which  may  thus  be  replaced,  for  the  purposes  of  analysis, 
by  a  parallel,  identical  virmal  array  excited  180“  out  of  phase  with  respect  to  the  real  array.  TTie  active  impedance  of 

the  z-th  dipole,  Z,.“ ,  is  therefore  given  by 

where  4  is  the  current  in  the  zr-th  dipole  and  Z-  is  the  mutual  impedance  between  dipoles  z  andy  minus  the  mutual 
impedance  between  dipole  z  and  the  virtual  image  of  dipoley  (Z,.  is  the  self-impedance  of  dipole  z  minus  its  mutual 
impedance  with  its  image)  [5].  The  Z^  depend  on  the  geometry  of  the  dipoles  and  their  mutual  geometnc  relations. 

If  the  array  element  excitations  are  all  in  phase,  then  for  efficient  radiation  the  active  impedances  Z/  should  all 
be  real  (and  are  hereafter  denoted  by  i?/).  This  additional  constraint  requires  the  simultaneous  optimization  of  the 
currents  1,  and  the  array  and/or  dipole  geometry.  This  can  be  done  by  using  simulated  annealmg  [10]  to  minimize  a  cost 
function  including  terms  for  deviation  from  the  desired  pattern  and  for  the  imaginary  parts  of  the  active  impedances. 
In  addition,  easy  feed  matching  can  be  ensured  if  the  cost  function  also  includes  appropriate  terms  for  feed  network 
parameters.  Note  that  the  variation  of  element  geometry  means  that  the  radiation  pattern  cannot  be  factonzed  as  the 
product  of  an  array  factor  and  a  single  element  factor,  but  must  be  calculated  from  the  individual  patterns  of  each  dipole. 


3.  Example 

Consider  the  synthesis  of  a  pencil  beam  pattern  by  a  linear  array  of  twelve  centre-fed  dipoles  of  radius  0.0048/i 
laid  out  in  the  x-y  plane  with  a  distance  All  one  from  another  with  their  centres  on  the  y  axis.  The  array  centre  is  at  the 
origin,  and  a  ground  plane  lies  ztz  =  -h  =  -XJ4.  In  the  y-z  plane  ((])  =  n/2),  the  radiation  pattern  is  to  have  a  maximum 
sidelobe  level  of -18  dB  (this  maximizes  the  peak  directivity  of  the  twelve-element  Dolph-Chebyshev  solunon  we 
assume  as  our  starting  point).  The  dipoles  are  to  be  fed  in  phase  from  a  single  10  Q  main  transmission  line  via  quarter- 
wave  coaxial  cables  and  no  extra  circuitry.  We  assume  that  the  radiation  patterns  of  the  dipoles  are  given  m  terms  of 

their  lengths  2/,-  by 


-375  - 


F.(6,(}))=sm(fc/i  COS0)- 


cos  (W.  sin  6  cos  (j))  -  cos  (kl.) 
V 1  -sin^0  cos'4) 


(2) 


where  k  =  2ti/A;  and  that  their  self  and  mutual  impedances  are  as  given  in  [5]  (this  implies  that  dipole  ieneths  can  be 
varied  only  in  the  interval  [1.3/^',  1.7/k]  in  which  the  assumed  expression  for  self-impedance  is  vaJid). 

To  verify  pattern  specification  when  optimizing  element  currents  and/or  geometry,  the  cost  runction  to  be 
minimized  can  include  the  terms  (A5^)^  and  c^/D^  ,  where  AS^  is  the  difference  between  the  caictiiated  level  of 

X 

the  side  lobe  s  and  the  maximum  tolerated  level,  if  this  difference  is  positive  and  AS^=Q  otheruise,  and  D.  is  the  peak 

directivity  of  the  pattern  (c,  and  c,  are  appropriate  weights).  To  ensure  that  the  array  elements  have  real  active 

impedances,  the  cost  function  can  include  the  term  ,  where  X,"  is  the  imaginary  pan  of  z/  (agaicu  c-  is 

/ 

an  appropriate  weight).  Finally,  the  feed  specification  means  that  the  feed  network  is  determined  by  the  characteristic 
resistances  of  the  quarter- wave  branch  lines.  We  have  /?,“=(/?//?.)*' ,  where  /?,.  are  the  resistances  presented  to  the 

main  feed  line  by  the  branch  lines  and  are  calculated  by  solving  the  equations  [5] 

R,/Rj=(ijRjy(i-R;)  (3) 

together  with  the  feed  matching  condition 

(4) 

To  meet  the  requirement  that  the  j?"  have  realistic  values,  say  <  140  Q,  the  cost  function  can  include  the 
term  J]  g,,  where  g,. =/?.“- 140,  if  this  is  positive,  andg.  =  0  otherwise,  with  an  appropriate  weight  If  the  I  are  all 

I 

given  the  fixed  value  A/4,  the  radiation  pattern  specification  is  satisfied  by  the  Dolph-Chebyshev  solution  [5]  listed  m 
Table  1.  This  solution  is  inefficient  because  of  the  coupling  indicated  by  the  complex  active  impedances  of  the  dinoles 
(also  listed  in  Table  1). 

A  solution  with  the  same  current  set  that  radiates  efficiently  and  maintains  the  directivity  can  be  found  by 
perturbing  the  /,  with  the  currents  fixed.  The  cost  function  C=Cj^  is  used,  to  move  to  zero 
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the  imaginary  part  of  the  active  impedances  while  penalizing  directivity  loss.  However,  if  simultaneous  satisfaction  of 
the  feed  specification  is  attempted  by  using  the  cost  function  C=c^{Asy*cJD  with 

grRi  -140,  no  solution  is  attainable.  To  obtain  the  solution  given  in  Table  2,  it  is  necessary  to  raise  the  limiting  values 
of  the  R’’"  to  145  Q  (i.e.,  g,- =/?,“- 145),  which  provides  characteristic  impedances  associated  with  coaxial  cables 
difficult  to  manufacture. 

By  contrast,  if  the  lengths  2/,.  and  currents  /.•  are  simultaneously  varied  by  using  the  same  cost  function  (with 
g,=/?“-140),  then  the  desired  power  pattern,  real  active  impedances  and  more  reasonable  i?"  values  can  all  be 
achieved  together  (Table  3). 


i 

/,• 

2l;/X 

z°(Q) 

1;12 

1.000 

0.500 

89.7+55.4; 

2;11 

0.665 

0.500 

110.1-4.1; 

3;10 

0.826 

0.500 

91.6+28.1; 

4;9 

0.964 

0.500 

101.2+20.8; 

5;8 

1.064 

0.500 

97.0+23.2; 

1  6;7 

1.117  1  0.500 

98.6+22.5; 

Peak  directivity:  57.9 

laOtC  l.  CACJUlUUtia,  oina  u.w.irN^  ....f, -  -- 

elements  of  a  12-dipole  array  when  the  dipole  length  is  kept 
uniform  (=AJ2)  and  the  excitations  are  optimized  to  generate  a 
Dolph-Chebyshev  pattern  with  -18  dB  side  lobes.  Peak  directivity 
and  branch  line  characteristic  impedances  are  also  shown. 
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/,• 

21/X 

z;(Q) 

Rr(Q)  ! 

1;12 

1.000 

0.461 

70.9 

96.6  i 

2;11 

0.665 

0.498 

107.0 

145.2  1 

3;10 

0.826 

0.479 

81.7 

117.0  i 

4;9 

0.964 

0.482 

90.4 

100.3  j 

5;8 

1.064 

0.481 

86.9 

90.8  i 

6;7 

1.117 

0.481  88.2 

86.4  1 

Peak  directivity:  57.5  Computation  time:  1.2  min.  1 


Table  2.  Excitations,  lengths  and  active  impedances  of  the  elements 
of  a  12-dipole  array  after  perturbation  of  the  dipole  lengths  of 
Table  1  to  afford  real  active  impedances  and  the  best  possible  set  of 
branch  line  characteristic  impedances,  keeping  excitations  fixed. 
Computation  time  refers  to  calculations  on  a  Pentium  II  running  at 
266  Mhz. 
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/. 

21, a 

Rr(^) 

1;12 

1.403 

0.466 

lA.l 

119.1 

2;11 

1.426 

0.485 

95.4 

117.2 

3;10 

1.456 

0.479 

82.7  i  114.8 

4;9 

1.444 

0.490 

95.7  !  115.8 

5;8 

2.092 

0.475 

83.1 

79.9 

6;7 

1.875 

0.484 

90.1 

89.1 

Peak  directivity;  57.0  Computation  time:  2.1  min 

Table  3.  As  for  Table  2,  but  with  simultaneous  optimization  of 
dipole  lengths  and  excitations. 
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Abstract-A  new  technique  for  the  synthesis  of  linear  array  antenna 
patterns  is  presented.  This  technique  allows  to  fix  nulls  in  prescribed 
directions,  to  simulate  the  presence  of  failed  elements  and  to  obtain 
power  patterns  when  both  situations  occur  at  the  same  time.  These 
patterns  axe  synthesized  by  finding  the  optimal  configuration  of  the 
array  factor  roots  using  the  simulated  annealing  technique.  Examples 
of  fixing  nulls  and/or  failed  elements  in  both  sum  and  flat-topped  beam 
patterns  are  presented. 

1.  INTRODUCTION 

Given  a  hnear  array  antenna,  there  e.xist  well-known  analytic  tech¬ 
niques  that  can  be  used  to  find  a  radiation  pattern  with  specific  re¬ 
quirements  about  beamwidth  and  sidelobe  level.  However,  if  some  el¬ 
ements  fail,  no  analytic  means  exist  to  obtain  an  aperture  distribution 
that  compensates  for  the  degradation  of  the  pattern.  The  defective 
elements  destroy  the  symmetry  and  cause  sharp  variations  in  the  field 
intensity  across  the  array  aperture,  increasing  the  side  lobe  level  of 
the  power  pattern.  In  the  literature,  there  exist  some  solutions  for 
this  problem:  in  the  case  of  receiving  antennas  the  degradation  can  be 
offset  by  transforming  the  received  pattern  [1],  another  approach  uses 
neural  networks  to  produce  a  new  array  transfer  function  and  offers 
advantages  in  direction  finding  in  a  single  signal  environment  [2].  In 
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the  case  of  transmitting  antennas,  the  currents  of  the  non-defective  el¬ 
ements  ran  often  be  adjusted  to  produce  a  pattern  with  minimal  loss  of 
quality  with  respect  of  the  original  pattern  [3—5].  This  compensation 
for  the  defective  elements  can  be  achieved  by  ntimerically  finding  the 
excitation  of  each  non-defective  element  that  optimizes  some  objective 
function. 

As  fax  as  we  know,  however,  pattern  recovery  techniques  have  not 
hitherto  been  developed  for  or  applied  to  the  problem  of  recovering 
patterns  in  which  nulls  have  been  fixed  in  certain  specific  directions. 
The  synthesis  of  array  patterns  to  obtain  nulls  in  prescribed  directions 
while  also  having  some  other  generally  desirable  pattern  features  (like 
low  sidelobe  level,  beamwidth,  etc.)  has  received  considerable  attention 
[6-9],  but  not  the  problem  of  maintaining  fiLxed  nulls  in  the  presence 
of  failed  elements. 

In  this  paper,  we  have  used  the  simulated  annealing  technique  [10] 
to  minimize  an  objective  function  in  order  to  obtain  a  desired  linear 
array  power  pattern  with  some  defective  elements  and/or  with  some 
fixed  nulls.  The  objective  function  used  takes  into  account  several 
specifications  such  as  side  lobe  level,  directivity,  dynamic  range  ratio 
(llmox/Iminl)  and  beamwidth. 

2-  THEORY 

The  array  factor  of  a  linear  array  of  iV  -f- 1  radiating  elements  laid  out 
at  equal  intervals  d  along  the  z  axis  with  its  center  at  the  origin  is; 

N  jV 

f{w)=yi  n 

n=0  n=l 

where  In  is  the  complex  excitation  of  the  n-th  element,  w  =  with 
Tp  =  kdcos{6),  where  k  is  the  wavenumber  and  9  is  the  angle  from 
endfire,  and  the  roots  Wn  may  be  written  in  the  form: 

Wn  =  (2) 

Fixing  nulls  in  certain  directions  means  fixing  a  subset  W'o  of  roots  of 
F{w)]  for  each  root  Wi  in  Wo,  ai  is  left  in  its  initial  value  (af)  and  bi 
depends  on  the  angular  position  of  the  corresponding  null.  The  failure 
of  element  n  means  that  \In\  =  0.  We  denote  the  set  of  defective 
elements  by  H. 
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Given  a  set  Wo  of  null- fixing  roots,  a  set  of  defective  elements  H 
and  a  desired  side  lobe  level  SLL^,  it  is  often  possible  to  calculate,  for 
the  non-defective  elements,  excitations  producing  a  pattern  with  the 
desired  side  lobe  level  and  nulls,  good  peak  directivity  D,  low  dynamic 
range  ratio  |Jmai/-lmtnl  and  also  to  control  the  beamwidth.  This  can 
be  achieved  by  using  the  simulated  annealing  technique  to  minimize 
the  following  objective  function: 

C  =Ci\SLL  -  SLL£  +  Cn\Ima./Imin\- 

+  ^-^C4Y^\Ik\''-T  cs\BW  -  (3) 

keH 

where  SLL  is  the  side  lobe  level,  and  BW  and  BWj,  axe  the  cur¬ 
rent  beamwidth  in  the  optimization  and  the  desired  beamwidth  re¬ 
spectively.  The  coefficients  Ci  are  weights  controlling  the  optimization 
process.  These  coefficients  must  be  selected  according  to  the  required 
specifilcations,  i.e.,  presence  of  defective  elements,  fixed  nulls  or  both 
situations  at  the  same  time.  Among  all  the  weights,  ci  is  one  of 
the  most  important  ones  because  the  final  pattern  must  have  a  desired 
side  lobe  level  that  is  usually  the  most  difficult  specification  to  achieve. 
However,  ail  the  weights  are  closely  coimected. 

In  the  following  sections  we  analyze  ail  the  cases  that  the  presented 
algorithm  can  simulate:  null  fixing  in  prescribed  directions,  element 
failure  simulation  and  the  case  in  which  both  situations  occur  together. 

2.1  Fixed  Nulls 

The  nulls  in  a  power  pattern  are  related  with  the  roots  of  the  array 
factor,  F{w)  since  their  angular  positions  axe  determined  by  the  bn 
factors.  Then,  it  is  possible  to  fix  nulls  in  a  desired  directions  (angular 
positions)  associated  to  specific  bi’s.  Therefore  in  the  optimization 
process  it  is  necessary  to  specify  previously  both  the  location  of  the 
fixed  nulls  (which  root  will  contain  each  fixed  root)  and  their  angular 
positions.  The  values  of  On  axe  left  in  their  original  values  to  avoid 
unwanted  null-filling.  The  mathematical  expression  for  the  roots  vari¬ 
ation  is: 

bn  =  bi  'iwn  e  Wo 

bn  =  6°  -f  5bn  '^VJn  ^  Wo  (4) 

On  =  Vne{l,2,...,iV} 

where  the  superscript  o  indicates  the  original  value  of  On  and  bn, 
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the  suDcrscript  /  the  pre-fixed  value  of  and  dbfi  the  perturbation 
introduced  to  6° . 

With  this  perturbation  of  the  roots,  the  minimization  of  the  ob¬ 
jective  function  (3)  is  performed  in  order  to  obtain  the  desired  power 
pattern.  In  the  general  objective  function  we  set  the  weight  C4  =  0 
since  the  simulation  of  failed  elements  is  not  needed  in  this  case. 

2.2  Failed  Elements 

To  simulate  failed  elements  no  restrictions  are  imposed  to  the  vari¬ 
ation  of  the  roots,  i.e.,  all  the  roots  are  moved  and  therefore  both  an 
and  bn  for  each  root  are  perturbed.  This  is  necessary  because,  in  gen¬ 
eral,  there  will  be  some  roots  off  the  Schelkunoff  circle  when  the  power 
pattern  has  failed  elements  [11].  Therefore,  the  variation  of  the  roots 
is  given  by: 

In  the  objective  function  (3)  the  weight  C4  must  be  chosen  large  enough 
to  ensure  small  values  of  the  excitations  associated  to  the  failed  sim¬ 
ulated  elements  (represented  in  the  H  set).  The  other  weights  are 
chosen  according  to  the  desired  specifications  of  side  lobe  levels,  dy¬ 
namic  range  ratio,  directivity  and  beamwidth. 

2.3  Fixed  Nulls  and  Failed  Elements 

With  the  presented  algorithm  it  is  possible  to  obtain  power  pat¬ 
terns  with  failed  elements  while  maintaining  fixed  nuUs.  To  search  the 
solution  eq.  (5)  is  used,  maintaining  invariable  the  set  of  roots  associ¬ 
ated  with  fixed  nulls.  The  minimization  of  the  objective  function  (3) 
is  performed  using  a  high  value  for  the  weight  C4  as  in  the  previous 
subsection. 

3.  RESULTS 

In  this  section  we  present  some  results  chosen  between  all  the  examples 
performed  to  test  the  method  proposed.  The  method  was  applied 
to  a  sum  pattern  and  to  a  fiat-topped  beam  pattern  and  the  tests 
include  fixing  of  nulls,  failed  elements  and  failed  elements  with  fixed 
nrdls.  At  the  end  of  each  power  pattern  analysis  we  present  a  table 
summaxizing  the  results  about  the  desired  parameters  to  optimize,  and 
the  computation  time  required  for  this  optimization.  This  time  has 
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been  measured  in  a  personal  computer  with  a  Pentium  11  processor 
mnning  at  266  MHz. 

3.1  Sum  Power  Pattern 

As  a  first  example  we  considered  a  40-element  (A/2)-spaced  lin¬ 
ear  array  producing  a  Dolph-Chebyshev  sum  pattern  with  a  side  lobe 
level  of  —25  dB  [12].  Other  parameters  of  interest  are  the  beamwidth 
at  -3  dB,  9.5°  (in  -w),  the  dynamic  range  ratio,  3.16,  and  the  peak 
directivity  36.58. 

3.1.1  Fixing  Nulls 

The  first  operation  performed  was  to  fix  nulls  in  certain  directions. 
The  selected  directions  were  ip  =  14°  and  ib  -  15°  which  are  nearby 
the  main  lobe.  In  Figure  1  it  is  shown  the  resulting  power  pattern 
after  the  optimization  process.  The  arrow  indicates  the  fixed  nulls. 
The  corresponding  excitations  are  shown  in  Table  1. 

3.1.2  Failed  Elements 

Starting  from  the  initial  power  pattern,  we  simulated  the  presence 
of  defective  elements  trying  to  decrease  as  much  as  possible  the  excita¬ 
tions  associated  with  the  “failed”  elements  while  maintaining  other  fac¬ 
tors  of  interest  in  the  design,  such  as  sidelobe  level  and  beamwidth,  the 
most  invariable  as  possible.  As  an  example  we  simulated  the  failure  of 
the  elements  6*^^  and  32“^^.  These  elements  had  normalized  amplitude 
excitations  of  0.51  and  0.66  respectively  in  the  initial  Dolph-Chebyshev 
power  pattern.  If  the  excitations  of  these  elements  are  directly  set 
to  zero  the  pattern  suffers  important  degradations,  as  it  is  shown  in 
dashed  line  in  Figure  2,  with  an  increase  of  the  sidelobe  level  of  about 
5  dB.  However,  if  we  start  from  the  initial  pattern  and  optimize  the 
objective  function,  we  obtain  the  power  pattern  shown  in  solid  line  in 
Figure  2.  In  Table  2  the  final  currents  (amplitude  and  phase)  of  this 
optimization  are  shown. 

Other  examples  were  performed  with  successful  results,  for  instance 
we  obtained  a  good  solution  for  the  problem  of  simulating  the  presence 
of  a  failure  in  the  20*^^  element  which  is  the  most  excited  one.  The 
final  sidelobe  level  obtained  was  —24.3  dB. 
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3.1.3  Fixed  Nulls  Together  with  Failed  Elements 

In  the  following  example  four  elements  were  considered  as  failed 
elements:  the  34'^^,  35*^^,  36'^^  and  the  37'*"  (their  initial  normal¬ 
ized  amphtude  excitations  were  0.57,  0.31,  0.46  and  0.41  respectively). 
Furthermore  we  fixed  two  nulls  in  il>  =  —83.6'’  and  tb  =  —82.6°.  The 
degraded  and  the  optimized  power  patterns  are  shown  in  Figure  3  and 
the  obtained  excitations  in  Table  3. 

The  information  about  the  optimized  parameters  for  each  exam¬ 
ple  performed,  sidelobe  level,  dynamic  range  ratio,  peak  directivity, 
beamwidth  and  computation  time  is  shown  in  Table  4. 

3.2  Flat-Topped  Beam  Power  Pattern 

In  this  example  we  used  a  40-element  (A/2)-3paced  linear  array 
producing  a  starting  flat-topped  beam  pattern  with  -30  dB  side  lobes 
and  four  filled  nulls  with  ±0.5  dB  ripple.  The  beamwidth  at  -3  dB  is 
52°  (in  ij)  and  ll-max/ Imvn\  is  3.35.  This  initial  pattern  was  obtained 
by  the  Orchard-Elhott  method  [13]. 

3.2.1  Fixing  Nulls 

The  first  example  performed  in  this  subsection  consists  in  fixing 
two  nulls  in  the  positions  -37°  and  -35°  {ip  values),  nearby  the 
main  lobe.  The  roots  that  contain  these  nulls  are  the  16'^  and  the 
17'^.  After  the  optimization  process  the  pattern  obtained  is  shown  in 
Figure  4  where  the  resulting  dip  is  indicated  with  an  arrow.  In  Table  5 
the  excitations  of  the  resulting  pattern  are  shown. 

3.2.2  Failed  Elements  Simulation 

In  this  case  we  simulated  the  presence  of  three  defective  elements 
in  the  positions  3,  4,  and  5  that  had  an  initial  normalized  amplitude 
excitations  of  0.40,  0.48  and  0.53  respectively.  The  degraded  power 
pattern  with  these  failed  elements  is  shown  in  dashed  line  in  Figure  5. 
It  can  be  seen  that  the  sidelobe  level  has  increased  up  to  -20  dB,  10 
dB  above  the  initial  sidelobe  level.  After  performing  the  optimization 
process  a  new  power  pattern  was  obtained,  plotted  with  solid  line  in 
Figure  5,  that  has  again  a  sidelobe  level  of  —30  dB.  The  corresponding 
excitations  (amplitude  and  phase)  are  shown  in  Table  6. 
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Figure  3.  Degraded  and  optimized  power  patterns  with  four  defective 
elements  (34,35,36,37)  and  two  fixed  nulls  in  ip  =  -83.6°  and  ip  = 
-82.6°. 

Table  3.  Excitations  for  a  sum  power  pattern  with  4  defective  ele¬ 
ments:  the  34‘^,  35*^*',  36"^^,  and  37'*^,  and  two  fixed  nulls  in  ip  = 
-83.6°  and  ip  =  -82.6°. 


Amplitude  |  Phase  (deg) 


mm 


Phase  (deg) 


0 


0.836 


roeiag 


0.969 


0 


0.927 


21.344 


21.723 


19.785 


0.995  ! 

i  -20.388  ! 

I  32  i 

1  0.320  1 

15.962 


-20.619 

34 

0.000 

-27.772 

35 

0.000 

0.000 

-24.631 

37 

0.000 

-25.816 

38 

0.250 

-25.430 

39 

0.250 

-24.892 

40 

0.506 

-25.658 


-35.536 


-32.826 


1.294 


-23.517 


1.270 


-16.468 


-25.224 
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Table  4.  Resulting  dynamic  range  ratio,  peak  directivity,  beamwidth 
and  computation  time  in  the  sum  power  pattern  examples. 


Failed 

elements 

Fixed  nulls 
position  (0) 

m 

Peak 

directivity 

Beamwidth 

at  —3  dB 

■ 

— 

14°,  15° 

3.58 

35.76 

1-24’’ 

6th_32nd 

— 

4.45 

33.43 

9.75° 

7’11” 

20th 

_ 

6.09 

33.47 

6’42” 

34"‘',35'^, 

36‘^,3T‘'' 

-83.6°, 

-82.6° 

4.00 

31.93 

10.75° 

7’46” 

3.2.3  Fixed  Nulls  and  Failed  Elements 

The  example  performed  consisted  on  j&xing  a  null  in  -0  =  -36.2° 
and  simulate  the  11'^  element  as  defective  (that  had  a  normalized 
excitation  amplitude  of  0.88).  The  degraded  power  pattern  with  the 
element  as  defective  and  the  recovered  pattern  with  the  presence 
of  the  fixed  null  are  shown  with  dashed  and  solid  lines  respectively  in 
Figure  6.  In  this  Figure  it  can  be  observed  that  the  failed  element 
causes  an  initial  degradation  in  the  power  pattern  in  such  a  way  that 
the  fixing  null,  in  spite  of  being  located  in  the  desired  direction,  it 
cannot  be  seen  in  the  degraded  pattern.  In  the  optimized  pattern, 
also  showed  in  Figure  6,  the  fixed  null  is  indicated  with  an  arrow.  The 
excitations  of  this  recovered  flat-topped  beam  power  pattern  are  showm 
in  Table  7. 

Finally,  the  results  about  the  sidelobe  level,  dynamic  range  ratio, 
beamwidth  and  computation  times  for  these  flat-topped  beam  exam¬ 
ples  axe  shown  in  Table  8.  ■ 

4.  CONCLUSIONS 

The  simulated  annealing  technique,  along  with  an  appropriate  way 
of  varying  the  roots  of  the  array  factor,  allows  to  obtain  power  pat¬ 
terns  with  fixed  nulls,  defective  elements  and  power  patterns  with  fixed 
nulls  in  presence  of  defective  elements.  Furthermore,  the  obtained 
patterns  have  suitable  properties  about  sidelobe  level,  beamwidth,  dy¬ 
namic  range  ratio  and  peak  directivity. 

It  can  be  observed  that  the  recovery  of  the  patterns  with  clustered 
failed  elements  is  better  than  if  they  are  distributed  more  uniformly 
along  the  array.  This  is  related  to  the  fact  that  the  pattern  distortion 
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- Optimized  pattern 

—  Degraded  pattem 


Figure  5.  Degraded  and  optimized  flat-topped  beam  pattern  with 
three  fcdled  elements  (3'^'*,  4*^^,  and  5‘^). 

Table  6.  Excitations  to  produce  a  desired  flat-topped  beam  pattern 
with  three  defective  elements:  the  3^^,  4'^^,  and  5'^. 


Phase  (deg)  n 


21 


26 


-76.041  I  27  I 


261.166 


236.431  1  32 


193.874  13510. 


Phase  (deg) 


128.794 


138.155 


itaMUKi 


181.265 


212.264 


227.222 


235.591 


240.794 


244.654 


255.071 


265.840 


-67.611 


-33.803 


-1.845 


2.357 


m 


0.135 


Null  fixin 


in.  and  failed  elements  in  linear  arrays 


-  Optimixecl  paKem  1 
Degraded  pattern 


T..  6  De<^aded  and  recovered  Eat-topped  beam  power  pattern 

U-  ?:Lt  defective  and  a  fixed  nnU  rn  =  -36  2  .  ^ 

Table  T.  Excitations  for  a  fiat-topped  beam  pattern  with  a  f^ure  m 
iaoie  I.  „  r  t  in  Tb  =  -36.2  . 


Amplitude 

fueTTZlL 

0.163 

0.202 

Olg 

0.195 

Phase  (deg)  T 

- r679t 

- ZTiOT^ 

5.15/ _ 

- -098^ 

- 0  179 

n 

_23 _ 

22 

23 

24 

25 

26 

Amplitude 

tut! 

0372^ 

0.763 

0333^ 

0.999" 

TOGO" 

0.195 

0.116 

O.lOO 

O.lOO _ 

O.lOO 
— oM 

M 

"  =030^ 

- =0368"'^ 

- — 37:99r" 

riT~ 

28 

29 

03^ 
1  0307" 

1  025"^ 

1  iusr: 
1  0322^ 

T  0447" 

; - ThTTsstT 

- 

1  30 
31 

1  32 

Phase  (deg) 

9  HI)  .662" 
240-361 

— '  99.9.9  i0~~ 

217.41Y“ 

- 216.058' 

- 213:^ 

- 219.719 

- - 225.816" 

I  9:-tH.87T 

1  958.488 

1 - 

1 - 

-1 - Z1S9.8TT 
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Table  8.  Resulting  dynamic  range  ratio,  peak  directivity,  beamwidth 
and  computation  time  in  the  flat-topped  beam  power  pattern  examples. 


Failed 

elements 

Fixed  nulls 
position  {ijj) 

H 

Beamwidth 

at  —3  dB 

Computation 

time 

— 

li^l 

52° 

2’33” 

3rd,  4th,  5th 

— 

7.39 

51° 

9’ir’ 

11th 

-36.2° 

10.03 

52° 

lO’lT” 

caused  if  all  the  defective  elements  are  close  together  consists  largely 
in  a  relatively  localized  increase  in  side  lobe  level,  whereas  the  increase 
is  general  if  there  are  defective  elements  throughout  the  array. 

Fixing  nulls  in  the  power  pattern  also  has  effects  on  the  convergence 
of  the  algorithm,  but  not  as  important  as  simulating  elements  failure. 
The  examples  presented  demonstrate  the  capability  of  the  algorithm 
of  fixing  nulls  very  close  to  the  main  lobe  or  the  shaped  region,  which 
is  a  considerable  advantage  for  some  applications. 

Finally,  with  slight  changes  in  the  algorithm  it  is  possible  to  obtain 
power  patterns  with  symmetric  fixed  nulls,  and  the  computation  time 
decreases  drastically.  The  key  factor  is  to  optimize  the  pattern  in  a 
symmetric  way,  i.e.,  only  taking  into  account  half  of  the  array  factor 
roots.  Nevertheless,  if  failed  elements  are  introduced,  the  results  are 
worse  than  in  the  presented  algorithm  because  the  simulated  aimeahng 
uses  fewer  variables  to  find  the  optimal  solution. 
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Abstract 

Given  the  failure  of  one  or  more  elements  of  an  array  antenna  designed  to 
produce  a  radiation  pattern  including  nulls  in  arbitrarily  fixed  directions,  a 
pattern  with  the  desired  side  lobe  level  and  nulls  can  often  be  recovered  by 
modifying  the  excitations  of  the  non-defective  elements.  The  modified 
excitations  can  be  calculated  by  simulated  annealing  using  an  appropriate  cost 
function. 

1.  INTRODUCTION 

The  radiation  pattern  produced  by  an  array  antenna  is  in  general  severely 
degraded  by  the  failure  of  even  a  very  small  number  of  the  array  elements.  In  the 
case  of  receiving  antennas,  the  degradation  can  be  offset  by  transforming  the 
received  pattern  [1].  In  the  case  of  transmitting  antennas,  the  excitations  of  the 
non-defective  elements  can  often  be  adjusted  to  produce  a  pattern  with  minimal 
loss  of  quality  with  respect  to  the  original  pattern  [2-4].  As  far  as  we  know, 
however,  pattern  recovery  techniques  have  not  hitherto  been  developed  for  or 
applied  to  the  problem  of  recovering  patterns  in  which  nulls  have  been  fixed  in 
certain  specific  directions.  The  problem  of  null-fixing  itself  has  received 
considerable  attention  [5-8],  but  not  the  problem  of  maintaining  fixed  nulls  in  the 
presence  of  element  failure.  In  this  paper  we  report  that,  for  linear  arrays,  fixed 
nulls  and  side  lobe  levels  can  both  be  recovered,  with  minimal  increase  in 
beamwidth  and  in  the  dynamic  range  ratio  of  the  element  excitations,  if  the 
excitations  of  non-defective  elements  are  recalculated  by  simulated  annealing 
using  an  appropriate  cost  function. 

2.  THEORY 

The  array  factor  of  a  linear  array  of  N+l  radiating  elements  laid  out  at  equal 
intervals  d  along  the  z  axis  with  its  centre  at  the  origin  is 

(1) 

n=0  «=1 

where  7„  is  the  complex  excitation  of  the  n-th  element,  H’  =  exp(/i//) 
(\l/=kd  cos{9),  where  k  is  the  wavenumber  and  6  the  angle  from  endfire),  and 
the  roots  w„  may  be  written  in  the  form  w„  =  exp(a„  +jbn).  Fixing  nulls  in  certain 


directions  means  fixing  a  subset  Wo  of  roots  of  F(w);  for  w,-  e  Wo,  a,-  =  0  and  bi 
depends  on  the  angular  position  of  the  corresponding  null.  The  failure  of  element 
n  means  that  /„  =  0;  we  denote  the  set  of  defective  elements  by  H. 

Given  a  set  Wo  of  null-fixing  roots,  a  set  of  defective  elements  H  and  a 
desired  side  lobe  level  SLLd,  it  is  often  possible  to  calculate,  for  the  non¬ 
defective  elements,  excitations  producing  a  pattern  with  the  desired  side  lobe 
level  and  nulls,  good  peak  directivity  D  and  low  dynamic  range  ratio  l-fmax/fininl. 
This  can  be  achieved  using  simulated  annealing  to  rninimize  the  cost  function 

C  =  Cj  |5LL  —  SLLj  I  +  C2  max  ^  ^  min  \  ~jZ  ■*"  ^4  X  |  (2) 

L>  ksH 

where  SLL  is  side  lobe  level  and  the  c,-  are  weights  controlling  the  optimization 
process  (C4  must  be  large  to  ensure  that  the  "optimized  excitations"  of  the 
defective  elements  are  effectively  zero);  the  optimization  is  carried  out  by 
introducing  slight  perturbations  in  the  values  of  a„  and  bn  for  non-fixed  roots. 


v  (deg) 

Fig.l.  Sum  power  pattern  with  two  fixed  nulls  at 
\|r  =  14°  and  Y  =  15°. 


3.  EXAMPLES 

As  a  first  example,  we  consider  a  40-element  (A/2)-spaced  linear  array 
required  to  produce  an  initial  sum  power  pattern  with  -25  dB  side  lobe  level.  In 
this  array  two  nulls  in  \j/  =  14°  and  \|;  =  15°  were  fixed.  The  final  obtained 
pattern  is  shown  in  Figure  1.  Other  examples  including  fixed  nuUs  together  with 
failed  elements  (to  be  shown  in  the  symposium)  were  performed  using  the  initial 
sum  power  pattern.  Table  1  summarises  the  obtained  results  about  the  sidelobe 
level,  dynamic  range  ratio,  beamwidth  and  computation  times,  measured  in  a 


Fig.  2.  Degraded  and  recovered  flat-topped  beam  power 
pattern  with  the  11***  element  defective  and  a  fixed 
null  in  \i;=-36.2‘’ 


Pentium  II  at  266Mhz,  for  this  and  other  sum  power  pattern  examples 
performed. 

Another  example  consisted  in  using  the  same  linear  array,  but  producing  a 
starting  flat-topped  beam  pattern  with  -30dB  side  lobes  and  four  filled  nulls  with 
±0.5  dB  ripple.  In  Figure  2  it  is  shown  the  degraded  and  the  recovered  power 
pattern  when  the  1 F'’  element  is  simulated  as  defective  and  a  null  is  fixed  m  t]/  = 
-36.2°.  In  Table  2,  there  are  shown  the  results  about  this  and  other  flat-topped 
beam  power  pattern  examples. 


Failed  elements" 

Fixed  nulls 
position  (\g) 

JJmax/Jimiil 

.  Peak , 
dir^tivity 

Beamwidth  at 
-3dB 

Computation 

time 

— 

14°,  15° 

3.58 

35.76 

9.50° 

r  24” 

4.45 

33.43 

9.75° 

7’  11” 

20* 

— 

6.09 

33.47 

10.00° 

6’ 42” 

34*,35*,36*.37‘’ 

-83.6°,  -82.6° 

4.00 

31.93 

10.75° 

7’ 46” 

Table  1.  Dynamic  range  ratio,  peak  directivity,  beamwidth  and  computation  time  in 


several  sum  power  pattern  examples. 


Failed  elements 

Fixed  nuUs 
position  (\j/)  . 

Ilmax/lTniDl 

/■  .  > 

Beamwidth  at 
-3dB " '  . 

Computation  time 

— 

-37°,-35° 

3.41 

52° 

2’  33” 

— 

7.39 

51° 

9’  17” 

11* 

-36.2° 

10.03 

52° 

10’  17” 

Table  2.  Dynamic  range  ratio,  peak  directivity,  beamwidth  and  computation  time  in 
some  flat-topped  beam  power  pattern  examples. 


4.  FINAL  REMARKS 

1)  As  was  done  in  synthesizing  the  pattern  of  Fig.l,  the  method  described  in 
Section  2  can  be  used  to  fix  nulls  in  the  patterns  of  non-defective  antennas. 
This  approach  to  nuU-fixing  can  fix  nulls  very  close  to  the  main  beam. 

2)  When  more  than  one  element  fails,  the  pattern  recovered  by  the  method 
described  in  this  communication  is  generally  better  if  the  defective  elements 
form  a  cluster  than  if  they  are  distributed  more  uniformly  along  the  array  (see 
Tables  1  and  2).  This  is  related  to  the  fact  that  the  pattern  distortion  caused 
if  all  the  defective  elements  are  close  together  consists  largely  in  a  relatively 
localized  increase  in  side  lobe  level,  whereas  the  increase  is  general  if  there 
are  defective  elements  throughout  the  array. 
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Abstract-Methods  previously  developed  for  designing  linear  array 
monopulse  antennas  with  relatively  simple  feed  network  requirements 
axe  extended  to  allow  the  suppression  of  jamming  signals  in  both  sum 
and  difference  modes. 

1.  INTRODUCTION 

High-performance  monopulse  array  antennas  must  generate  both  sum 
and  difference  patterns  with  low  side  lobes,  high  directivity  and  narrow 
beamwidth.  They  are  also  often  required  to  operate  in  situations  in 
which  intentional  or  unintentional  Jamming  signals  are  being  received 
from  certain  directions.  To  counteract  such  interferences,  both  the  sum 
and  difference  patterns  of  the  antenna  must  feature  nulls  or  quasi-nulls 
(highly  depressed  responses)  in  these  directions.  Although  this  can  cer¬ 
tainly  be  achieved  by  synthesizing  the  sum  and  difference  patterns  by 
means  of  excitation  distributions  that  are  mutually  independent  [1], 
this  approach  maximizes  the  complexity  and  cost  of  feed  networks. 
Haupt  [2]  showed  how  nulls  can  be  introduced  into  given  sum  and  dif¬ 
ference  patterns  with  minimal  modification  requirements  by  assigning 
each  element  a  phase  shift  or  complex  factor  modifying  both  its  sum 
and  difference  excitations,  but  assimed  that  for  high  pattern  quality 
the  original  sum  and  difference  patterns  would  have  to  be  generated 
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by  independent  excitation  distributions.  For  lineax  arrays  with  2N 
elements  it  is  possible  to  synthesize  sum  and  difference  patterns  with 
common  arbitrary  nulls  using  a  sum  distribution  that  is  periodic  with 
period  N  and  a  difference  distribution  that  is  the  same  as  the  sum 
distribution  except  for  a  phase  inversion  affecting  half  the  elements, 
but  the  resulting  pattern  quality  can  be  rather  poor  [3].  Lee  [4]  sought 
a  compromise  between  pattern  qualit^’”  and  feed  network  simplicity  by 
modifying  Chebyshev  and  Bayliss  distributions  so  that  the  outer  ele¬ 
ments  have  the  same  excitations  for  both  sum  and  difference  modes 
(except  for  the  phase  inversion  affecting  one  side  of  the  distribution) 
while  the  sum  and  difference  excitations  of  the  inner  elements  are  ob¬ 
tained  by  phase-shifting  the  original  Chebyshev  and  BayHss  excitations 
respectively;  however,  this  compromise  clearly  favours  pattern  quality 
more  than  feed  network  simplicity. 

In  previous  work  [5,  6]  we  developed  synthesis  methods  that,  with¬ 
out  tackling  the  jammer  neutralization  problem,  significantly  reduced 
the  feed  network  complexity  of  monopulse  array  antennas  by  forcing 
certain  excitation  parameters  to  be  shared  among  elements  or  between 
modes.  Specifically,  we  showed  that  by  using  simulated  annealing  [7] 
to  optimize  appropriate  objective  functions,  high-quality  patterns  can 
be  synthesized  in  which  either  the  excitation  ampHtude  of  each  element 
is  the  same  for  both  the  sum  and  difference  patterns  ( “phase-only  con¬ 
trol”)  or  the  difference  pattern  is  obtained  by  multiplying  the  sum 
excitations  of  all  the  elements  in  each  of  a  series  of  subaxrays  by  a 
single  subarray  weight  (Fig.  1).  In  this  paper  we  report  the  exten¬ 
sion  of  these  methods  to  allow  the  introduction  of  arbitrauy  quasi-nulls 
common  to  the  sum  and  difference  patterns. 


2.  METHODS 


2.1  Phase-Only  Control 

We  assume  that  the  phase-only  control  method  described  in  [5]  for  pla¬ 
nar  arrays  has  been  used  to  calculate  the  amplitudes  A; ,  sum  pattern 
phases  ai  and  difference  pattern  phases  Pi  of  compromise  excitation 
distributions  for  a  linear  array  of  2N  elements,  and  that  the  sum  and 
difference  distributions  are  respectively  symmetric  and  antisymmetric 
about  the  centre  of  the  array  so  as  to  generate  symmetric  sum  and  dif- 
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4  (0)4' (6) 

Figure  1.  Subarray  method:  configuration  of  the  element  and  subar¬ 
ray  weighting  devices  of  one  side  of  a  monopulse  finear  array  antenna. 


ference  patterns.  We  wish  to  introduce  quasi-nulls  in  given  directions 
6j  (measured  from  broadside).  In  keeping  with  the  "phase-only'’  phi¬ 
losophy,  we  do  so  by  altering  the  phases  ai  and  Pi  by  increments  5ai 
and  5Pi  respectively,  and  we  maintain  the  symmetry  or  antisymmetry 
of  the  distributions.  The  increments  Sai  and  SPi  are  calculated  by 
using  simulated  annealing  to  optimdze  the  cost  function 

C(Sa,  5p)  =  Cl  •  4-  [SLX^o  - 

+  C2  ■  E  {[r 4-  [r 
j 

(1) 

where  da  and  56  are  the  vectors  of  dOj  and  Spi ,  SLL  indicates 
a  side  lobe  level  and  F  {dj)  the  value  of  the  power  pattern  in  the 
direction  dj  ,  subscript  S  indicates  the  sum  distribution  or  pattern 
and  subscript  A  the  difference  distribution  or  pattern,  su’bscript  d 
indicates  the  desired  value  of  a  variable  or  parameter  and  subscript  o 
the  value  obtained  with  the  current  arguments  5a  and  5P ,  and  ci 
and  C2  are  weighting  factors  controlling  the  balance  between  main¬ 
tenance  of  desired  side  lobe  levels  and  suppression  of  interferences  in 
the  directions  dj  .  A  broad  quasi-null  extending  over  a  finite  angular 
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interval  [^staxti  ^end]  can  be  obtained  by  specifying  a  large  number  of 
point  quasi-nulls  9j  within  this  interval. 

2,2  Subaxraying 

In  the  subarray  method  for  linear  arrays  described  previously  [5,  6] 
the  sum  pattern  is  generated  by  a  Taylor  or  Dolph-Chebyshev  distri¬ 
bution  and  the  real  subarray  weights  are  optimized  to  make  the 
difference  pattern  as  good  as  possible  using  a  cost  function  involving 
only  difference  pattern  characteristics.  This  approach,  which,  leaves  the 
original  sum  pattern  unchanged,  is  not  possible  when  nulls  or  quasi¬ 
nulls  are  to  be  imposed  in  both  the  difference  and  sum  patterns.  In 
this  situation,  for  a  2iY-element  array  the  cost  function  must  also  op¬ 
timize  the  sum  excitations  Ai  (constrained  to  be  real  and  symmetric 
about  the  centre  of  the  array)  and  control  the  characteristics  of  both 
the  sum  and  difference  patterns.  A  suitable  cost  function  is 

C{A,  p)  =  Cl  •  [[SLL^o  -  SLL^d]-  d-  [5LLao  -  5LLAdl“} 

+  C2  ■  E 

0 

-f-  C3  •  {\A  1}  (2) 

where  A  =  (Ax ,  •  -  • ,  A/^Y,  g  =  (pi,  •  •  • ,  gqYiQ  is  the  number  of  sub¬ 
arrays  on  each  side  of  the  array  centre,  about  which  the  difference 
excitation  distribution  is  constrained  to  be  antisymmetric),  Amax  and 
Amin  are  respectively  the  greatest  and  least  of  the  Ai,  B^ax  and  Bmm 
are  respectively  the  greatest  and  least  of  the  Bi  =  ,  and  C3  is  a 

weight  analogous  to  ci  and  co  ■  The  optimization  method  is,  as  usual, 
simulated  annealing. 

3.  EXAMPLES 

In  this  section  we  illustrate  the  above  techniques  by  applying  them  to 
a  Linear  array  of  50  equispaced  elements  one  half  wavelength  apart.  In 
ail  cases,  the  desired  power  level  specified  for  quasi-nulls  was  —60  dB  , 
which  is  in  most  situations  sufficient  to  neutralize  jammers.  The  sum 
and  difference  patterns  plotted  in  Figs.  2-6  have  been  normalized  to 
their  respective  peak  directivities.  None  of  the  calculations  took  more 
than  5  min  on  a  personal  comouter  with  a  Pentium  II  running  at  266 
MHz. 
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Figure  2.  Sum  and  difference  patterns  of  a  5 0-eiement  antenna  ob¬ 
tained  by  the  phase-only  control  method  [5]  specifying  —20  dB  side 
lobes  but  no  fixed  quasi-nuUs. 


3.1  Phase-Only  Control 

Fig.  2  shows,  as  functions  of  ■0  =  vrsin^,  the  compromise  sum 
and  difference  patterns  that  were  obtained,  without  specification  of 
nulls  or  quasi-nulls,  by  the  phase-only  control  design  method  previously 
described  for  planar  arrays  [-5].  The  side  lobe  levels  of  the  sum  and 
difference  patterns  are  —20.8  and  —19.7  dB  respectively,  and  the 
dynamic  range  ratio  |.d.maxMinin|  is  1.29. 

Fig.  3  and  Fig.  4  show  the  results  of  using  the  method  of  Section  2.1 
to  modify  the  excitation  distributions  corresponding  to  the  patterns 
of  Fig.  2  so  as  to  introduce  quasi-nulls  at  0  =  40°,  80°  and  120° 
(Fig.  3)  or  throughout  the  interval  [50°,60°j  (Fig.  4).  The  side  lobe 
levels  of  the  sum  and  difference  patterns  are  now  —19.6  and  —17.0  dB 
respectively  in  the  former  case,  and  respectively  —17.7  and  —14.4  dB 
in  the  latter. 

3.2  Subarraying 


For  best  results,  application  of  the  subarraying  technique  to  the  50- 
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Figure  3.  Sum  and  difference  patterns  obtained  by  using  the  phase- 
only  method  to  introduce  quasi-nulls  sA  ip  =  40°,  80°  and  120°  in  the 
patterns  of  Fig.  2. 


Figure  4.  Sum  and  difference  patterns  obtained  by  using  the  phase- 
only  method  to  introduce  a  broad  quasi-null  over  the  interval  50°  < 
Ip  <  60°  in  the  patterns  of  Fig.  2. 
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Figure  5.  Sum  and  difference  patterns  obtained  using  the  subarray 
method  to  fix  quasi-nulls  ■w  =■  40°,  80°  and  120°  while  ensuring 
side  lobe  levels  of  —  30  dB  . 


Figure  6.  Sum  and  difference  patterns  obtained  using  the  subarray 
method  to  fix  a  broad  quasi-null  over  the  interval  50°  <  tp  <  60° 
while  ensuring  side  lobe  levels  of  —30  dB  . 
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element  array  reqmred  the  use  of  four  subarrays  in  each  half  of  the 
antenna:  one  comprising  two  elements  and  the  other  three  compris¬ 
ing  three  each  (the  remaining  fourteen  elements  on  each  side  receive 
the  same  excitation  for  both  sum  and  difference  patterns).  Fig.  5 
shows  the  results  obtained  when  quasi-nulls  were  specified  at  ^  = 
40°,  80°  and  120°  and  side  lobe  levels  of  -30  dB  were  required  in 
both  sum  and  difference  patterns.  The  specified  side  lobe  levels  have 
been  achieved  for  both  patterns,  with  dynamic  range  ratios  of  1.32  and 
1.93  for  the  sum  and  difference  patterns  respectively. 

Fig.  6  shows  the  results  of  subarraying  to  achieve  a  broad  quasi-null 
covering  the  interval  50°  <  0  <  60°  together  with  side  lobe  levels 
of  -30  dB  in  both  patterns.  Again,  both  these  goals  have  been  met, 
with  dynamic  range  ratios  of  1.56  and  1.87  for  the  sum  and  difference 
patterns  respectively. 

3,3  Addendum:  Phase  and  Amplitude  Perturbation 

In  view  of  the  relatively  large  rise  in  side  lobe  levels  when  a  broad 
quasi-null  was  introduced  using  phase-only  control  (Section  3.1),  we 
resorted  to  perturbation  of  both  phases  and  amplitudes,  adding  the 
term  C4  •  (i/  -  max{15Ail/.4i})^  to  the  cost  function  so  as  to  limit  the 
relative  magnitude  of  the  perturbations.  .Table  1  lists  the  side  looe 
levels  and  dynamic  range  ratios  achieved  with  various  specified  values 
of  u  (the  desired  quasi-nuU  was  achieved  in  ail  cases).  As  might  be 
expected,  as  v  increases  the  side  lobe  levels  improve  but  the  dynamic 
range  ratio  also  increases. 


u  =  max{l<5Ai|/Ai} 

SLLzo 

SLL^o 

lAmax/ Atninl 

10% 

-18.6  dB 

-15.4  dB 

1.30 

25% 

—19.5  dB 

-16.4  dB 

1.36 

50% 

-19.9  dB 

—  17.7  dB 

1.57 

Table  1.  Side  lobe  levels  and  dynamic  range  ratios  achieved  by  per¬ 
turbing  both  the  phases  and  amplitudes  of  the  distributions  corre¬ 
sponding  to  Fig.  2  so  as  to  introduce  a  broad  quasi-null  over  the  inter¬ 
val  50°  <-ip  <  60°  ,  for  various  values  of  u  in  the  cost  function. 
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4.  CONCLUSIONS 

Of  the  two  techniques  extended  in  this  paper  to  allow  specification  of 
common  quasi-nulls  in  sum  and  difference  patterns,  the  better  app- 
pears  to  be  the  subaxraying  technique,  at  least  as  regards  achievement 
of  side  lobe  level  specifications.  This  may  be  attributed  to  its  reduc¬ 
ing,  by  means  of  its  subarray  weights,  the  discontinuity  present  in  the 
unweighted  difference  excitation  distribution  at  the  centre  of  the  ar¬ 
ray.  As  regards  ease  of  implementation  for  antennas  that  must  be  able 
to  adapt  to  changing  jamming  environments,  the  phase-only  control 
method  requires  2N  fixed  amplitude  attenuators  calculated  for  the 
starting  patterns,  together  with  either  2N  fixed  phase-shifcers  and 
2iY  variable  phase-shifters,  or  4iV  variable  phase-shifters;  the  subar¬ 
ray  method  requires  2iY  +  2Q  variable  amplitude  attenuators.  Both 
methods  appear  to  conserve  pattern  quality  significantly  better  than 
do  previously  published  methods  using  feed  networks  of  comparable 
complexity.  Both  can  be  extended  to  planar  arrays. 
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ABSTRACT 

We  describe  the  use  of  the  simulated  annealing  optimization  method  to 
introduce  minimal  modifications  in  the  excitations  and/or  geometries  of 
planar  and  circular  arc  antenna  arrays  so  as  to  depress  the  radiation  pattern  in 
ajbitrary  directions  while  maintaining  side  lobe  levels  and  beamwidth. 


1.  INTRODUCTION 

Given  an  array  antenna,  null  fixing  consists  in  modification  of  its  parameters 
so  that  the  radiated  power  is  zero  (null),  or  at  least  highly  depressed  (quasi¬ 
null),  in  certain  specified  directions.  Most  null-fixing  methods  have  been 
developed  for  linear  arrays  (see,  for  example,  (Trastoy  and  Ares,  1997)).  (Vu, 
1992)  devised  two  methods  for  fixing  nulls  in  planar  arrays,  but  the  first, 
based  on  the  method  of  separable  distributions,  only  controls  pattern 
parameters  in  the  two  principal  ^2J-cuts,  while  the  second,  which  uses  a 
cancelling  beam  as  in  some  techniques  developed  for  linear  arrays,  has  only 
been  applied  to  uniform  excitation  distributions;  neither  of  these  two  methods 
exercises  any  control  over  the  dynamic  range  ratio  of  the  array  excitation 
distribution,  |4ax/4inl.  or  is  able  to  null  a  finite  angular  segment  of  the  pattern 
rather  than  just  a  point  or  series  of  points.  Nulls  have  also  been  fixed  in  the 
radiation  patterns  of  circular  arc  arrays,  both  by  perturbing  the  excitations  of 
the  airay  elements  (Prasad  and  Charan,  1984),  (Vescovo,  1996),  (Ares  et  al, 
1996)  and  by  perturbing  their  radial  positions  (Hejres  and  Richie,  1996).  The 
excitation  perturbation  method  of  (Ares  et  al,  1996)  has  the  advantage  over 
the  others  of  minimizing  dynamic  range  ratio  and  the  extent  of  perturbation 


at  the  same  time  as  it  fixes  nulls. 

In  this  article  we  extend  to  planar  and  circular  arc  arrays  the  null-fixing 
method  described  previously  for  linear  arrays  (Trastoy  and  Ares,  1997). 
Briefly,  a  simulated  annealing  technique  based  on  the  simplex  method  (Press 
et  al,  1992)  is  used  to  optimize  perturbations  of  parameters  characterizing 
the  array  (specifically,  the  amplitudes,  phases  and/or  geometric  positions  of 
its  elements)  so  as  to  introduce  desired  quasi-nulls  while  controlling  other 
radiation  pattern  characteristics  (side  lobe  levels,  beamwidth,  etc.)  and 
minimizing  the  dynamic  range  ratio  of  the  excitation  distribution. 


2.  DESCRIPTION  OF  THE  METHOD 

We  start  with  an  antenna  array  of  given  geometry,  composed  of  radiating 
elements  with  given  excitations  and  positions.  By  perturbing  these  array 
parameters  1  i.e.  the  excitations  and/or  positions  of  the  elements,  the  radiation 
pattern  can  be  changed  so  that  quasi-nulls  appear  in  certain  desired  directions. 
In  general,  this  can  be  achieved  by  an  infinite  number  of  perturbation 
combinations:  the  aim  is  to  identify  a  combination  of  perturbations  that  not 
only  fixes  the  required  quasi-nulls,  but  also  maintains  or  improves  other 
pattern  characteristics  (side  lobe  levels,  beamwidth,  etc.)  and  facilitates 
implementation  by  minimizing  both  the  dynamic  range  ratio  of  the  excitation 
distribution  and  the  absolute  values  of  the  perturbations.  This  can  be  achieved 
by  using  a  suitably  versatile  optimization  procedure  such  as  simulated 
annealing  (Press  et  al,  1992)  (see  Section  3)  to  minimize  a  cost  function 
containing  a  group  of  terms  fixing  pattern  parameters  (including  quasi-nulls), 
a  group  of  terms  minimizing  perturbations,  and  a  term  minimizing  the 
dynamic  range  ratio: 

C^i2a,(Dd,-Dof*b\I^/I^f*c'Z\SAf  (1) 

J  =  1  "  =  1  , 

where  Ddi  and  Do,-  denote  respectively  the  desired  value  of  radiation  pattern 
parameter  i  (power  level  in  the  direction  where  the  quasi-null  is  required,  side 
lobe  level,  beamwidth,  etc.)  and  the  value  afforded  by  the  current  array 
parameters;  q  indicates  the  number  of  these  parameters  that  are  used;  5A„  is 
the  current  perturbation  of  the  n-th  array  parameter;  and  the  a,-,  b  and  c  are 
weighting  factors  controlling  the  relative  priorities  of  the  various  terms  in  the 
cost  function. 


3.  THE  OPTIMIZATION  METHOD 


The  cost  function  defined  in  equation  1  is  a  function  of  a  large  number  of 
variables,  upon  which  each  of  its  terms  depends  in  a  complex  fashion.  Such 
functions  generally  have  a  large  number  of  local  minima,  and  it  is  not  usually 
possible  for  them  to  be  minimized  by  conventional  optimization  procedures 
based  on  their  local  topography;  these  methods  will  usually  not  find  the 
global  minimum,  but  the  local  minimum  downhill  from  the  starting  point, 
i.e.  the  minimum  of  that  local  basin  (of  the  hypersurface  defined  by  the 
function)  which  contains  the  point  used  as  the  starting  point  of  the 
optimization  procedure.  The  optimization  technique  called  simulated 
annealing  overcomes  this  difficulty  in  a  fashion  inspired  on  the 
thermodynamical  process  it  is  named  after.  Broadly,  the  system  (i.e.  the 
current  best  estimate  of  the  position  of  the  minimum)  is  first  allowed  to  roam 
randomly  from  one  local  basin  to  another,  with  a  tendency  -  but  not  an 
absolute  obligation  -  to  occupy  successively  lower  positions;  then  the 
"energy"  of  this  motion  is  gradually  reduced  (by  decreasing  a  "temperature" 
parameter)  so  that  the  system  eventually  settles  into  a  single  basin;  and  finally 
the  lowest  point  of  that  basin  is  found,  usually  by  a  conventional  optimization 
method.  Like  other  "evolutionary"  optimization  procedures,  simulated 
annealing  cannot  guarantee  to  find  the  global  minimum  of  a  function,  but  it 
will  usually  find  a  much  better  local  minimum  than  any  conventional  method. 
Of  the  various  forms  of  simulated  annealing  that  have  been  proposed,  the 
variant  used  to  obtain  the  results  reported  below  in  Section  4  was  the 
algorithm  published  by  (Press  et  al,  1992),  which  is  based  on  the 
conventional  deterministic  optimization  procedure  known  as  the  simplex 
method  (Press  et  al,  1992).  In  minimizing  a  function y(x)  of  n  variables,  the 
simplex  method  does  not  start  with  just  a  single  initial  estimate  of  the  point 
at  which  the  minimum  lies,  but  with  n+1  estimates  Xi,...,x„+j  which  jointly 
define  a  simplex,  the  n-dimensional  equivalent  of  a  plane  triangle.  At  each 
iteration,  the  procedure  first  identifies  the  x,-  for  which  the  value  of  the 
function  is  greatest,  x^^,,  say,  and  implicitly  adopts  the  working  hypothesis 
that  the  face  of  the  simplex  that  is  opposite  to  this  point  divides  high  values 
from  low.  It  accordingly  constructs  a  new  estimate  Xo  lying  on  the 
prolongation  of  the  line  from  x^^^  to  the  centre  of  the  opposite  face,  and 
proceeds  to  test  the  validity  of  this  estimate  and  to  act  in  consequence.  For 
example,  if /(Xq)  is  less  than  all  the  other /(x,.)  (suggesting  that  a  good  line  of 
descent  has  been  found),  a  new  point  Xqo  is  constructed  further  along  the  same 
line  as  before,  the  simplex  is  re-defined  (excluding  and  including  the 
better  of  Xg  and  Xgo),  and  the  next  iteration  is  begun.  Depending  on  the  results 


of  similar  comparisons,  the  simplex  thus  changes  in  position,  shape  and 
volume.  Eventually,  it  becomes  smaller  and  smaller  until  all  the  function  /has 
very  similar  values  at  all  its  vertices,  at  which  point  the  least  of  these  values 
is  taken  as  the  minimum. 

When  the  simplex  method  is  used  for  simulated  annealing,  the  procedure  is 
the  same  as  above  except  that  comparisons  between  points  are  biased  in 
favour  of  acceptance  of  a  new  candidate  for  inclusion  in  the  simplex.  For 
example,  the  criterion  for  constructing  Xqo  is  no  longer /(x,)  >/(Xo)  for  all  i, 
hutfix;)  +  Y;  >y(Xo)  for  all  i,  where  7,-  is  a  random  variable  with  a  G(2,7) 
distribution,  7 being  the  temperature  parameter.  When  Tis  large,  then  so,  on 
average,  are  the  7,-,  which  means  that  Xqo  is  almost  always  constructed  (and 
then  almost  always  accepted).  The  simplex  thus  grows  until  prevented  from 
growing  further  by  meeting  high  values  ofJ{x)  as  \x\  grows.  Lowering  T every 
few  tens  of  iterations  reduces  the  probability  of  accepting  x^  and  increases 
the  probability  of  the  simplex  becoming  smaller  and  settling  in  a  local  basin. 
At  T  =  0  the  7;  are  all  zero  and  the  procedure  is  exactly  the  standard  simplex 
method. 


4.  RESULTS 

4.1  Application  to  Planar  Arrays 

To  establish  notation,  we  recall  that  the  far  field  F(6,(p)  of  a  planar  array  of 
N  elements  lying  in  the  (x,y)  plane  is  given  by 

F(e.<p)  =  E  (2) 

n  =  l 

where  and  y„  are  the  position  coordinates  of  the  n-th  element,  7„  and  a„  are 
respectively  the  amplitude  and  phase  of  its  excitation,  /  is  the  wavenumber 
and  FE(6,(p)  is  the  element  pattern.  We  take  as  an  example  the  case  of  a  468- 
element  array  composed  of  52  3x3-element  subarrays  of  circular  patches 
arranged  on  an  8x8-subarray  grid  with  three  subarrays  omitted  at  each  comer 
so  as  to  form  an  approximately  circular  global  array  boundary.  The  power 
pattern  of  each  circular  patch  is  given  by  which  falls  to 

-22  dB  at  0  =  90°. 

Starting  from  a  uniform  excitation  distribution,  a  distribution  with  a  dynamic 
range  ratio  |4ax/^miiil  of  just  2.07  producing  a  power  pattern  with  a  side  lobe 
level  of  -26.0  dB  can  be  synthesized  by  perturbing  the  excitation  amplitudes 
and  phases  and  controlling  side  lobe  level  in  ^-cuts  taken  every  5°.  A  quasi- 


null  with  a  radius  of  1°  centred  on  {d,(p)  =  (20°, 40°)  can  now  be  introduced 
in  this  power  pattern  by  further  perturbing  the  excitation  amplitudes  and 
phases  while  using  a  cost  function  that  penalizes  not  only  side  lobe  level  and 
dynamic  range  but  also  departure  from  -60  dB  at  any  of  a  series  of  points 
within  the  desired  quasi-null  sector.  Fig.l,  in  which  u  =  sin(^)cos(^?)  and 
V  =  sin((9)sin((p),  shows  the  resulting  power  pattern,  which  has  a  side  lobe 
level  of  -23.3  dB  and  a  quasi-null  in  the  desired  direction  (arrowed,  at 
u  =  0.26,  V  =  0.23)  and  is  produced  by  an  excitation  distribution  with  a 
dynamic  range  ratio  of  2.38.  If  only  phases  are  perturbed  in  the  second  stage 
of  this  process,  a  distribution  with  a  dynamic  range  ratio  of  2.07  producing 
a  power  pattern  with  a  side  lobe  level  of  -20.7  dB  is  obtained  (Fig.2). 
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FIGURE  1 .  A  broad  (1°)  null  fixed  at 
(00.  (Po)  =  (20°,40°)  in  the  power 
pattern  of  a  planar  array  antenna  by 
perturbing  excitation  amplitudes  and 
phases. 
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HGURE  2.  A  broad  (1°)  null  fixed  at 
(00,  (po)  =  (20°,40°)  in  the  power 
pattern  of  a  planar  array  antenna  by 
perturbing  excitation  phases  only. 


4.2.  Application  to  Circular  Arc  Arrays 

The  far  field  of  an  array  of  N  radiating  elements  lying  on  equispaced  rays  in 
a  sector  (-(p,(p)  of  a  circle  is  given  in  the  plane  of  the  circle  by 

«<p)=E  4  ® 

n=l 

where  R„  and  (p^  are  the  polar  coordinates  of  the  n-ih  element  and  the  other 


notation  is  the  same  as  for  equation  2.  Fig.3  shows  the  geometry  when  R„-R 
for  all  n. 


FIGURE  3.  Geometry  of  a  circular  arc  array 


In  this  case,  the  element  pattern  cannot  be  factored  out  of  the  expression  for 
F((p),  leaving  an  array  factor,  because  each  element  faces  in  a  different 
direction.  As  an  example,  consider  the  case  N  =15  and  (p  =  60°,  with  the 
elements  initially  at  equal  intervals  of  XII  along  the  arc  (F„  =  5.73X  for  all  n) 
and 


FF(0,(p)  =|-sine  [1  +2max(cos(p,-i)] 


(4) 


(this  is  the  element  pattern  used  for  (V2)-spaced  elements  by  (Jiao  et  al, 
1993),  who  adapted  the  findings  of  (Herper  et  al,  1985)). 

Starting  from  the  excitation  distribution  found  by  excitation 
perturbation  in  (Ares  et  al,  1996),  which  has  a  dynamic  range  ratio  of  7.54 
and  produces  a  pencil  beam  with  a  side  lobe  level  of  -40  dB,  our  aim  is  to 
perturb  the  excitations  and/or  positions  of  the  elements  so  as  to  introduce 
quasi-nulls  at  ^  =  20°,  22°  and  24°,  thereby  creating  a  broad  quasi-null 
centred  on  ^  =  22°.  Table  1  and  Figs.  4-7  show  the  results  obtained  by 
varying  a)  excitation  amplitudes  and  phases  (Fig.4),  b)  excitation  phases  only 
(Fig.5),  c)  radial  and  angular  element  coordinates  (Fig.6),  and  d)  excitation 
phases  and  angular  coordinates  (Fig.7);  in  all  cases,  perturbations  were 
symmetric  with  respect  to  ^  =  0°  so  as  to  reduce  the  number  of  variables. 
Best  results  were  achieved  perturbing  phases  and  position  angles,  which 
introduced  the  desired  quasi-null  sector  at  the  expense  of  only  a  2.64  dB  rise 


in  side  lobe  level.  Note  that  the  position  angle  perturbations  were  small 
enough  for  equation  4  still  to  be  a  valid  approximation  for  the  element 
pattern. 


TABLE  1.  Results  of  Fixing  a  Broad  Null  at  (p=22“  in  the  Power  Pattern  of  a 
Circular  Arc  Array  by  Perturbing  Various  Combinations  of  Array  Parameters 


Paramenters  perturbed 

SLL 

A 

Amplitudes  and  phases  (Fig.4) 

-36.38  dB 

8.31 

0.5?. 

Phases  (Fig. 5) 

-34.32  dB 

7.54 

0.5?. 

Radial  and  angular  positions 
(Fig.  6) 

-36.84  dB 

7.54 

0.48X-0.51?. 

Phase  and  angular  position 
(Fig.7) 

-37.36  dB 

7.54 

0.471  -  0.521 

SLL=Side  lobe  level;  A=distance  between  adjacent  elements  along  the  arc. 


FIGURE  4.  A  broad  (4°)  null  fixed  at  (p=22°  in  the  power 
pattern  of  a  circular  arc  array  by  perturbing  excitation 
amplitudes  and  phases. 


(j)  (degrees) 


,  FIGURE  5.  A  broad  (4°)  null  fixed  at  9=22“  in  the  power 
pattern  of  a  circular  arc  array  by  perturbing  excitation 
phases  only. 


FIGURE  6.  A  broad  (4“)  null  fixed  at  (p=22°  in  the  power 
pattern  of  a  circular  arc  array  by  perturbing  radial  and 
angular  element  coordinates. 
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FIGURE  7.  A  broad  (4°)  null  fixed  at  9=22“  in  the  power 
pattern  of  a  circular  arc  array  by  perturbing  excitation 
phases  and  angular  element  coordinates. 

It  was  also  found  that  quite  good  results  could  be  obtained  by 
perturbing  the  parameters  of  a  relatively  small  number  of  array  elements. 
Starting  from  the  excitation  distribution  found  in  (Ares  et  al,  1996)  for  the 
equispaced  array  with  =  5.73/1  for  all  n,  perturbation  of  the  phases  of  the 
eight  elements  whose  phases  were  changed  most  in  synthesizing  the  pattern 
of  Fig.5  (four  on  each  side  of  the  array)  afforded  the  pattern  shown  in  Fig.8, 
which  still  has  a  side  lobe  level  of  -29.6  dB. 


FIGURE  8.  A  broad  (4°)  null  fixed  at  9=22“  in  the  power 
pattern  of  a  circular  arc  array  by  perturbing  the  excitation 
phases  of  just  8  of  its  25  elements. 


5.  CONCLUSIONS 


By  perturbing  the  excitations  and/or  positions  of  some  or  all  of  the 
radiating  elements,  the  method  described  here  is  able  to  fix  quasi-nulls  over 
a  finite  sector  of  the  power  patterns  of  planar  and  circular  arc  array  antennas 
while  controlling  other  power  pattern  characteristics  and  minimizing  the 
dynamic  range  ratio  of  the  excitation  distribution.  Though  not  shown  in  this 
article,  multiple  point  quasi-nulls  can  be  fixed  by  the  same  method.  Although 
the  perturbation  of  element  positions  is  of  course  more  difficult  to  implement 
physically  than  perturbation  of  excitations,  it  shares  with  phase-only  control 
the  advantage  of  allowing  no  increase  in  dynamic  range  ratio,  and  for  circular 
arcs  affords  better  results  than  perturbation  of  excitations  alone. 
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Abstract-  A  planar  array  design  method  of  fixing  nulls  is  presented.  This  method  allows  to  fix  narrow  an 
broad  nulls  by  varying  amplitudes  and  phases  or  only  phases  of  the  excitations  of  the  radiating  elemen  s. 
Besides  it  is  possible  to  control,  at  the  same  time,  the  sidelobe  level  and  the  dynamic  range. 


1.  Introduccion 

Aunque  las  tecnicas  de  fijado  de  nulos  en  arrays 
lineales  han  sido  tratadas  en  grait  niimero  de  articulos, 
no  sucede  lo  mismo  en  arrays  pianos.  Vu  [1]  presento 
dos  tecnicas  para  ello.  La  primera  de  ellas,  basada  en 
la  utilizacion  de  distribuciones  separables,  es  una 
tecnica  no  optima,  puesto  que  el  control  de  los 
parametros  del  diagrama  se  realize  solo  en  los  dos 
cortes  principales.  La  segunda,  que  es  similar  a 
tecnicas  propuestas  para  arrays  lineales,  esta  basada  en 
la  formacion  de  un  haz  de  cancelacion  y  presenta 
mejores  resultados  aunque  es  aplicada  solo  a 
distribuciones  de  excitacion  uniforme  (control  s61o  en 
fase).  Sin  embargo,  ninguna  de  ellas  permite  el  control 
del  rango  dinamico  ni  el  fijado  de  nulos 

anchos. 

La  tecnica  aqui  propuesta  consiste  en  la  extension  a 
arrays  pianos  del  metodo  descrito  en  [2]  para  arrays 
lineales.  Esta  extension  pemiite  la  sintesis  de 
diagramas  de  radiacion  de  arrays  pianos  mediante  el 
control  de  las  excitaciones  de  los  elementos  radiantes, 
bien  variando  conjuntamente  la  amplitud  y  la  fase  de 
dichas  excitaciones,  o  bien  modificando  unicamente 
la  fase.  Aunque  el  objetivo  primero  de  esta  tecnica  es 
el  fijado  de  nulos  (tanto  puntuales  como  en  un  cierto 
entomo  de  un  punto  deseado),  este  metodo  permite  la 
modificacion  de  otros  parametros  importantes  del 
diagrama  de  radiacidn,  como  pueden  ser  el  nivel  de 
lobulos  laterales  o  el  ancho  angular  del  haz  principal. 
Ademas,  tambien  permite  la  minimizaemn  del  rango 
dinamico  de  las  excitaciones,  lo  que  favorece  la 
realizabilidad  del  array  resultante. 


2.  Dcscrlpcion  del  metodo 

La  expresion  del  patron  de  radiacidn  en  campo  lejano, 
para  un  array  piano  de  N  elementos  identicos  situados 
en  un  piano  paralelo  al  piano  XY,  es  la  siguiente. 


N 


.  ,  y-pjr  jen6cos<t) 

r(e,<f))=E4"  ® 

n=l 


donde  I„  denota  la  amplitud  de  la  excitacion  del 


elemento  n-esimo,  a„  su  fase,  x„,  e  y„  denotan  su 
posicion  sobre  el  piano,  FE(0,  cj))  el  factor  elemento  y 
P=2k/X. 

Se  parte  de  un  diagrama  inicial  que  es  modificado  con 
el  fin  de  conseguir  otro  con  nulos  fijados  en  unas 
ciertas  posiciones  angulares  y  que  posea,  ademas, 
unas  caractcristicas  deseadas.  Para  ello  se  varian  las 
excitaciones  de  los  elementos  radiantes,  bien  sus 
amplitudes  y  sus  fases: 


a  =a_+6a 

n  n  n 

0  bien  solo  sus  fases: 

a  =a°  -t-6a 

n  n  n 

Es  importante  que  estas  variaciones  ^  en  las 
excitaciones  sean  minimas,  condicion  esta  que 
podemos  poner  ennuestra  tecnica,  ya  que  esto  perrmte 
una  mayor  simplicidad  de  la  red  de  alimentacion 
requerida  para  el  sistema  radiante.  Este  parametro, 
junto  con  otros  caracteristicos  del  diagrama  y/o  de  las 
excitaciones,  se  introduce  en  una  funcion  de  coste  que 
sera  minimizada  mediante  la  tecnica  de  optimizacion 
de  “Simulated  Annealing”[3].  Dicha  funcion  de  coste 
incluye,  entonces,  un  primer  termino  para  conseguir 
unas  ciertas  caracteristicas  en  el  diagrama  de  radiacion 
tales  como  la  posicion  angular  y  la  profundidad  de  los 
nulos.  el  ancho  de  haz  o  el  nivel  de  lobulos  laterales. 
Un  segundo  termino  permite  minimizar  el  rango 
dinamico  de  las  excitaciones,  y  un  tercer  sumando 
permite  conseguir  que  las  excitaciones  de  los 
elementos  radiantes  varien  lo  menos  posible  con 
respecto  a  las  del  diagrama  de  partida.  Es  decir. 

donde,  por  Dd.  y  Do.se  denotan  los  valores  deseado 
y  obtenido  de  un  cierto  parametro  i  del  diagrama  de 
radiacion,  p  indica  el  numero  de  estos  parametros  que 
queremos  controlar,  y  Ei^  los  valores  final  e 
inicial  de  las  amplitudes  y/o  fases  del  elemento 
n-esimo  y  aj,  b  y  c  son  los  factores  de  peso  que  nos 
permiten  dar  prioridad  a  unos  sumandos  sobre  los 
otros. 
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3.  Resultados 

Para  ilustrar  la  tecnica  que  se  presenta,  hemos 
escogido  un  array  piano  de  468  elementos  radiantes 
equiespaciados  A./2,  colocados  segiin  un  enrejado 
rectangular  de  8x8  subarrays  eliminando  tres 
subarrays  en  cada  esquina  para  aproximar  a  un 
conromo  circular.  Cada  subarray  consta  de  3x3 
parches  circulares  cuyo  factor  elemento,  dependiente 
de  0,  puede  aproximarse  por  la  siguiente  expresidn: 

Inicialmente  se  parte  de  un  diagrama  de  radiacion 
correspondiente  a  una  distribucion  de  excitaciones 
uniforme  y,  mediante  la  vaiiacion  de  las  amplitudes  y 
las  fases  de  las  excitaciones  de  cada  uno  de  los 
elementos  radiantes,  se  sinteti2a  im  diagrama  en  el  que 
se  fijaran  posteriormente  los  nulos  y  que  presenta  un 
nivel  de  lobulos  laterales  SLL=-26.0  dB  y  un  rango 
dinamico  minimizado  hasta  j=2.07.  Los  nulos 
pueden  ser  tanto  muy  localizados  en  una  posicion 
angular  deteminada,  como  nulos  mas  anchos  que  se 
consiguen  fijando  varios  de  ellos  en  un  entomo  de  la 
posicion  deseada.  Los  ejemplos  que  aqui  se  muestran 
son  de  este  ultimo  tipo.  En  la  figura  1  se  presenta  la 
proyeccion  en  el  piano  XY  (primer  cuadrante)  de  un 
diagrama  con  un  nulo  de  banda  ancha  que  ocupa  un 
entomo  de  radio  1°  alrededor  de  la  posicion 
(0o.<1>o)=(2O“,4O”),  es  decir,  en  u=0.26,  v=0.23  (donde 
u=sen9cos(j)  y  v=sen0sen(j))  y  que  se  ha  fijado 
mediante  la  variacion  de  amplitudes  y  fases  de  las 
excitaciones.  El  diagrama  obtenido  presenta  un  nivel 
de  lobulos  laterales  de  -23.3  dB  y  un  rango 
dinamico  de  2.38. 


Figura  1.  Nulo  centrado  en  (0o,(|)o)=(2O°,4O°) 
obtenido  mediante  la  variacion  de  amplitudes  y 
fases. 


En  la  figura  2  se  muestra  ese  mismo  diagrama 
obtenido  mediante  solamente  la  variacidn  de  las  fases 
de  las  excitaciones.  Este  diagrama  de  radiacion 


presenta  un  nivel  de  lobulos  laterales  de  -20.7  dB 
y,  obviamente,  su  rango  dinamico  es  el  mismo  que  en 
el  patron  de  partida;  2.07. 


O-O  0.2  0.4  0.6  0.8 


Figura  2.  Nulo  centrado  en  (6o,4>o)=(20‘’,40°) 
obtenido  mediante  la  variacion  de  fases. 


4.  Conclusiones 

Se  ha  presentado  una  tecnica  para  el  fijado  de  nulos, 
tanto  sencillos  como  multiples,  en  el  diagrama  de 
radiacion  de  arrays  pianos.  Ademas,  posibilita  el 
control  del  nivel  de  lobulos  laterales  en  todo  el 
diagrama  3D,  asi  como  la  minimizacion  del  rango 
dinamico  de  las  excitaciones  y  de  las  variaciones  de 
dichas  excitaciones  (en  amplitud  y  fase)  con  respecto 
a  las  correspondientes  al  diagrama  de  partida. 
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and  digitisers,  sampling  each  element  signal  at  four  samples/bit 
(1.08  MHz);  a  direction  finding  DSP  subsystem  which  pre-proc¬ 
esses  and  stores  the  digitised  signals;  a  survey  vehicle  containing  a 
standard  DCS-1800  mobile  station  (MS)  and  a  GPS  receiver;  and 
a  DCS-1800  base  station  (BSS)  to  set  up  and  maintain  ±e  call  and 
to  provide  a  Hming  reference  for  the  location  estimation  subsys¬ 
tem.  The  experiments  were  carried  out  by  driving  the  survey  vehi¬ 
cle  over  a  number  of  test  routes  with  a  standard  call  in  progress 
between  the  MS  in  the  vehicle  and  the  BSS  in  the  field  trial  sys¬ 
tem.  The  signals  received  at  the  antenna  array  were  processed  in 
real  time  to  estimate  H  and  R,  which  were  time  stamped  and 
stored  for  later  analysis. 


Fig.  1  Histogram  of  DOA  estimates  over  2s  period 

Raw  DOA  and  range  estimates  were  computed  with  the  tech¬ 
nique  described  under  Tolar  location  estimation’.  The  number  of 
sources  was  fixed  to  one.  The  mobile  trajectory  was  estimated 
with  the  method  based  on  the  histogram  maxima  with  time  inter¬ 
vals  of  2s  and  the  threshold  of  reliability  set  to  40%.  Avera^g 
was  then  performed  with  a  window  of  size  20s,  or  ten  estimates. 
An  example  of  the  DOA  histograms  obtained  from  analysis  of  the 
experimental  data  is  shown  in  Fig.  1.  This  example  was  obtained 
from  an  experiment  in  which  a  co-channel  interference  signal  was 
present  in  addition  to  the  target  mobile. 

Results:  Since  differential  GPS  equipment  was  not  available,  the 
GPS  measurements  were  subject  to  a  bias  that  varies  between 
experiments.  Knowing  the  exact  routes,  the  measured  OPS  posi¬ 
tions  were  overlaid  onto  an  Ordnance  Survey  map  and  the  cor¬ 
rected  GPS  trajectories  were  calculated. 


Table  1:  DOA  error  in  degrees,  range  error  in  metres  and  r.m.s. 
position  error  in  metres 


(Ref.)  scenario 

DOA  error 

Range  error 

Position 

P(err  <  Xm) 

Mean 

0^^ 

1^3 

RMSE 

X=125 

deg 

m 

m 

m 

% 

% 

% 

(039)  moving  MS 

-0.4 

-53 

98 

133 

24.5 

59.6 

98.4 

(045)  moving  MS 
stationary  interferer 

0.3 

0.6 

2 

68 

114 

36.4 

82.0 

97.9 

Last  three  columns  represent  probability  that  position  error  is  <  X, 
with  X  =  50,  125  and  250m,  respectively 


The  position  estimates  derived  from  the  mobile  Icx^tion  algo¬ 
rithm  are  subject  to  a  DOA  bias  (due  to  misalignment  of  the 
anteima  array  and  phase  and  amplitude  mismatch  between  the 
receiver  channels)  and  a  range  bias  (due  to  fixed  delays  in  the 
eight  receivers.)  These  bias  terms  were  estimated  using  data  from  a 
’calibration’  experiment  in  which  the  mobile  station  was  stationary 
at  a  location  with  line  of  sight  to  the  base  station.  The  resulting 
bias  corrections  were  applied  to  the  estimates  derived  from  other 
experimental  data  sets.  Performance  was  evaluated  with  the  mean 
error  and  standard  deviation  on  the  DOA  and  range  estimates, 
and  the  r.m.s  position  error.  The  accuracy  of  the  location  esti¬ 
mates  was  also  assessed  by  counting  the  number  of  position  esti¬ 
mates  within  circles  of  increasing  radius  centred  on  the  true 
positions.  Table  1  shows  the  results  for  two  scenarios.  In  the  first 
one,  the  mobile  station  is  in  a  moving  vehicle  travelling  at  an 


average  speed  of  30mph.  The  duration  of  the  test  is  ~  5min  and 
the  distance  from  the  base  station  is  -  5  km.  The  second  scenario  is 
the  same  except  that  a  stationary  co-channel  interference  signal  is 
introduced. 

In  the  first  scenario,  the  estimated  location  is  within  250m  of 
the  true  location  for  98%  of  the  time,  and  within  50m  for  24%  of 
the  time.  The  second  scenario  gives  a  larger  percentage  of  esti¬ 
mates  within  50m  and  comparable  results  for  250m. 

Conclusion:  This  Letter  has  presented  a  technique  for  estimating 
the  location  of  a  mobile  radio  transmitter  using  an  antenna  array 
at  a  single  receiving  point.  Experimental  results  from  a  DCS-1800 
system  have  been  presented  which  compare  favourably  with  other 
published  results.  Of  course,  position  accuracy  is  very  dependent 
on  the  propagation  environment  and  future  experiment^  work 
will  include  the  application  of  the  technique  to  a  wider  range  of 
propagation  environments  (including  micro-cell  environments),  as 
well  as  the  refinement  of  the  location  estimation  algorithms. 
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Pattern  synthesis  of  array  antennas  with 
additional  isolation  of  near  field  arbitrary 
objects 

L.  Landesa,  F.  Obelleiro,  J.L.  Rodriguez, 

J.A.  Rodriguez,  F.  Ares  and  A.G.  Pino 

A  new  technique  for  synthesising  array  antennas  in  the  presence 
of  arbitrary  obstacles  in  the  near  field  region  is  presented.  This 
technique  allows  a  prescribed  far  field  pattern  to  be  synthesised, 
while  trying  to  avoid  coupling  with  near  field  obstacles  through 
the  constraint  of  minimising  the  radiated  power  over  the  obstacle 
surface.  Some  results  have  been  obtained  showing  that  the 
obstacles  are  finally  isolated  from  the  array  anteima  radiation. 

Introduction:  In  recent  works  [1  -  3],  the  synthesis  of  array  anten¬ 
nas  with  the  presence  of  obstacles  in  the  near  field  environment 
has  been  investigated.  In  [1],  a  solution  based  on  a  pattern  synthe¬ 
sis  procedure  imposing  near  field  nulls  inside  the  obstacles  was 
presented  by  Steyskal;  this  straightforward  approach  provides 
good  results  for  small  obstacles,  although  it  may  be  inefficient 
when  dealing  with  large  objects.  A  new  solution,  presented  in  [2], 
tries  to  take  into  account  the  obstacles  as  parasitic  elements, 
including  them  in  a  pattern  synthesis  algorithm,  which  allows  us 
to  take  advantage  of  their  scattering,  even  for  large  obstacles. 
Directivity  optimisation  is  achieved  in  [3]  following  a  similar  cou¬ 
pling  procure  to  take  into  account  the  presence  of  the  obstacles. 
These  two  last  approaches  [2,  3]  need  to  calculate  accurately  the 
electromagnetic  coupling  of  the  obstacle  in  order  to  obtain  the 
array  excitations.  This  complicates  the  synthesis  procedure  and, 
besides,  requires  the  use  of  an  accurate  numeric^  technique  to 
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calculate  this  coupling  matrix,  mainly  in  complex  and  three- 
dimensional  problems. 

In  this  Letter,  we  present  a  new  technique  that  generalises  the 
approach  presented  in  [1]  for  large  objects.  This  technique  allows 
us  to  synthesise  a  prescribed  far  field  pattern  while  attempting  to 
avoid  the  coupling  with  near  field  obstacles  through  the  constrtiint 
of  minimising  the  radiated  power  over  the  obstacle  surface  by  a 
simple  and  straightforward  procedure.  This  technique  has  been 
shown  to  be  very  efficient,  providing  results  as  good  as  those 
obtained  with  the  previous  techniques  [2,  3]  in  which  the  scattering 
of  the  obstacle  was  used  to  improve  the  radiation  pattern. 


T 

Xi 

array  coefficients 


Fig.  1  Linear  array  of  isotropic  sources  in  presence  of  arbitrary  obstacle 


Pattern  synthesis  technique:  In  this  pattern  synthesis  problem  we 
attempt  to  obtain  those  array  complex  coefficients  x,  that  best 
approximate,  in  a  mean  least  square  sense,  a  desired  far  field  pat¬ 
tern  fj,  simultaneously  minimising  the  near  field  radiated  power  in 
the  boundary  of  the  obstacles,  B.  We  consider  a  two-dimensional 
problem  involving  a  linear  array  of  N  isotropic  sources  with  V2 
spacing,  as  shown  in  Fig.  1. 

The  boundary  of  the  obstacle  is  sampled  at  a  set  of  M  matching 
points  (ten  matching  points  per  wavelength  have  been  used)  and  a 
matrix,  C,  that  relates  the  free  space  electric  near  field  in  these 
points  with  the  array  coefficients  is  calculated.  Thus,  the  field  in 
these  matching  points  is  obtained  as  a  column  vector  of  M  coeffi¬ 
cients: 

e  =  Cx  (1) 

where  x  is  a  column  vector  of  N  elements  that  contains  the  excita¬ 
tion  coefficients  x,  for  the  array,  and  C  is  a  matrix  of  MxN ele¬ 
ments  of  the  form 


(2) 


where  r,y  notes  the  distance  between  the  zth  matching  point  and  the 
^  element  of  the  array.  The  power  magnitude  in  the  boundary  of 
the  obstacles  is  defined  as 

Pb  =  \\Cx\\l  (3) 


which  will  be  minimised  in  order  to  isolate  the  obstacles.  The  syn¬ 
thesis  problem  can  therefore  be  formulated  as 

mini  [  \fd(.u)  —  f{u)fdu\  constrained  to  ||Cxll2  <  a 


(4) 


where  u  =  sin0  and  /is  the  free-space  far  field  pattern  of  the  array 
whose  coefficients  are  given  by  x. 

Considering  the  Fourier  relation  between  the  aperture  distribu¬ 
tion  and  the  far  field  [1,  4],  and  using  Lagrange  multipliers,  the 
last  equation  can  be  written  as 

mm{||xd-xl|i+A2||Cx||^}  (5) 

where  the  parameter  X  is  related  to  a  in  a  nonlinear  sense,  and  the 
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colunm  vector  Xj  represents  the  free-space  array  coefficients  for 
the  prescribed  pattern  /  This  equation  can  be  reformulated  in  an 
overdetennined  system  as 

(ac)*=(o)  ® 

which  can  be  solved,  in  a  mean  least  square  sense,  by  using  the 
normal  equation  or  the  QR  decomposition. 


Fig.  2  Geometry  of  array  antenna  in  presence  of  PEC  cylinder 
a  =  9X,b  =  0.8X,  N^4l,d=  O.SX,  D  = 


Fig.  3  far  field  patterns 

-  -  -  •  origiaal  Taylor  pattern 
- blocked  Taylor  pattern 

.  radiation  pattern  synthesised  with  proposed  method  (free 

space) 

- radiation  pattern  obtained  with  proposed  method  (blocked) 


In  the  solution  of  eqn.  6,  the  parameter  X  (which  must  be  real 
and  positive)  weighs  up  the  importance  of  the  boundary  power 
{P^  restriction  in  the  synthesis  procedure.  Small  values  of  X  corre¬ 
spond  with  a  poor  isolation  of  the  obstacles;  in  this  case,  although 
the  pattern  pr^ction  may  be  good  (compared  with  the  prescribed 
pattern),  the  scattered  field  radiated  by  the  obstacle  highly 
degrades  this  prediction,  so  the  total  far  field  pattern  becomes 
poor.  Otherwise,  for  large  values  of  X,  the  problem  is  dominated 
by  the  boundary  power  restriction;  hence,  although  the  obstacle’s 
radiation  does  not  disturb  the  far  field  pattern,  this  pattern  predic¬ 
tion  may  be  far  from  the  prescribed  pattern.  It  is  clear,  therefore, 
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that  a  good  selection  for  X  is  required;  namely,  X  must  be  ade¬ 
quately  selected  for  balancing  a  good  isolation  of  the  obstacle  and 
a  pattern  prediction  close  to  the  prescribed  one. 

Results:  Eqn.  6  is  used  in  this  Section  in  order  to  approximate,  in 
a  mean  least  square  sense,  a  Taylor  pattern  of  SLL  =  40dB  and  n 
=  3,  with  a  linear  array  antenna  blocked  by  a  perfect  electric  con¬ 
ducting  (PEQ  cylinder  whose  geometry  and  dimensions  are  shown 
in  Fig.  2.  Fig.  3  shows  the  radiation  patterns  obtained  in  the  fol¬ 
lowing  cases:  (i)  prescribed  Taylor  pattern  without  obstacle,  (ii) 
prescribed  Taylor  pattern  blocked  by  the  obstacle,  (iii)  free^pace 
radiation  pattern  synthesised  with  the  proposed  method  {X  =  0.7) 
and  (iv)  radiation  pattern  obtained  with  the  proposed  method 
including  the  scattering  of  the  obstacle.  It  is  important  to  note  the 
good  agreement  between  (iii)  and  (iv)  which  demonstrates  that  the 
synthesised  array  does  isolate  the  obstacle;  otherwise,  (iv)  provides 
an  important  improvement  with  respect  to  the  blocked  array  (ii). 
In  (ii)  and  (iv),  the  scattering  from  the  obstacle  is  obtained  by 
determining  the  induced  currents  over  the  stirface  .S  via  a  method 
of  moments  formulation  [5]. 

Conclusions:  A  new  method  for  synthesising  array  antennas  in  the 
presence  of  arbitrary  obstacles  has  been  presented,  which  allows 
the  synthesis  of  a  prescribed  far  field  pattern  while  trying  to  iso¬ 
late  near  field  obstacles  through  a  minimisation  performed  over  a 
power  magnitude  in  their  boundary.  This  techmque  is  straightfor¬ 
ward  to  apply  and  does  not  require  the  use  of  numerical  electro¬ 
magnetic  methods  for  characterising  the  coupling  of  the  obstacle. 
It  has  been  shown  to  be  very  efficient,  providing  results  that  can 
be  as  good  as  those  obtained  with  previous  more  rigorous  tech¬ 
niques  [2].  Although  numerical  results  have  been  presented  for 
PEC  obstacles,  the  method  is  general  for  any  kind  of  material 
obstacles. 
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Small  slot-coupled  circularly-polarised 
microstrip  antenna  with  modified  cross-slot 
and  bent  tuning-stub 

Kin-Lu  Wong  and  Ming-Huang  Chen 

A  novel  design  of  small  slot-coupled  circularly-polarised  circular 
microstrip  antenna  with  a  modified  cross-slot  cut  in  the  patch  and 
a  bent  tuning-stub  aligned  along  the  patch  boundary  is  proposed 
and  experimentally  studied.  Results  show  that,  for  fixed  circular 
polarisation  (CP)'  operation,  the  antenna  proposed  can  have  an 
antenna  size  reduction  of -80%,  as  compared  to  a  regular-size  CP 
design.  The  significant  size  reduction  of  the  proposed  antenna  is 
due  to  the  novel  modified  cross-slot  cut  in  the  patch,  and  CP 
operation  is  obtained  using  a  bent  tuning-stub  incorporating  a 
properly  oriented  coupling-slot  in  the  ground  plane  of  the 
microstrip  feed  line. 

Introduction:  Owing  to  recent  requirements  in  the  miniaturisation 
of  personal  communications  systems,  the  design  of  small  or 
reduced-size  microstrip  antennas  has  become  an  important  sub¬ 
ject.  Several  reduced-size  CP  designs  have  also  been  reported  [1,  2] 
in  which  CP  operation  is  achieved  by  adjusting  the  embedded  slot 
in  the  patch  centre  or  the  inserted  slits  at  the  patch  edges  to  be 
slightly  asymmetric.  Such  reduced-size  CP  designs  usually  exhibit 
stringent  fabrication  tolerances  when  the  anterma  size  reduction  is 
>  50%.  This  is  because,  in  such  cases,  very  small  asymmetry  in  the 
embedded  slot  or  inserted  slits  is  required  to  achieve  CP  opera¬ 
tion.  To  ease  the  problem,  the  reduced-size  CP  design  with  a  sym¬ 
metric  cross-slot  embedded  in  the  patch  centre  has  been  shown  [3] 
in  which  CP  operation  is  achieved  by  using  a  straight  mning-stub 
loaded  at  the  patch  boundary,  and  the  required  tuning-stub  length 
for  CP  operation  at  a  large  antenna  size  reduction  can  be  as  large 
as  ~26%  of  the  patch’s  linear  dimension,  which  significantly  eases 
the  stringent  fabrication  tolerances  as  required  m  (1,  2].  In  this 
Letter,  we  demonstrate  another  promising  reduced-size  CP  design 


using  a  novel  symmetric  slot  embedded  in  the  patch,  and  the  CP 
operation  is  adjusted  by  using  a  bent  tuning-stub  aligned  along  the 
patch  boundary  to  make  the  proposed  antenna  more  compact  in 
total  size.  Excitation  of  the  proposed  antenna  is  through  a  prop¬ 
erly  oriented  coupling-slot  in  the  ground  plane  of  the  microstrip 
feed  line.  Details  of  the  proposed  antenna  are  described,  and 
experimental  results  of  the  CP  performance  are  shown. 


Fig.  1  Exploded  view  and  top  view  of  circularly-polarised  circular 
microstrip  antenna  with  modified  cross-slot  and  bent  tuning-stub 

a  Exploded  view 
b  Top  view 

Antenna  design  and  experimental  results:  Fig.  1  shows  the  geome¬ 
try  of  the  proposed  design.  The  circular  patch  has  a  diameter  of 
25  mm  (D)  and  is  printed  on  a  substrate  of  thickness  1.6mm  (/i,) 
and  relative  permittivity  4.4  (e,).  A  modified  cross-slot  of  equal 
slot  lengths  is  cut  in  the  patch.  The  slot  length  is  chosen  to  be 
greater  than  0.8  times  the  diameter  of  the  circular  patch  in  order 
to  result  in  a  large  meandering  of  the  excited  fundamental  patch 
surface  current  path  to  lower  the  resonant  frequency  of  the  micros 
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ABSTRACT 

A  new  method  for  the  array  pattern  synthesis  in  presence  of  near-zone  scatterers  in  the  two- 
dimensional  scalar  case  is  presented.  The  synthesis  is  based  on  the  minimization  of  an  objective 
function  using  the  simulated  annealing  technique.  To  take  into  account  the  effects  of  the 
scatterers  the  physical  optics  approximation  is  used.  Examples  for  a  41 -element  linear  array  in 
presence  of  perfectly  electric  conducting  cylinders  are  performed. 

INTRODUCTION 

The  pattern  perturbations  caused  by  a  scatterer  in  the  proximity  of  an  array  antenna  is  a  problem 
of  practical  interest,  since,  for  example,  masts  or  wing  tips  often  are  within  the  radiative  near¬ 
zone  of  a  ship  or  airborne  antenna.  In  the  literature  exist  some  methods  to  solve  this  problem  in 
the  two-dimensional  scalar  case.  For  instance,  in  [1]  an  adapted  pattern  that  avoids  illuminating 
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such  near-field  scatterers  on  transmit,  or  rejects  the  corresponding  scattered  waves  on  receive,  can 
result  in  a  significantly  improved  antenna  far-field  response.  However,  this  method  imposes  a 
more  severe  field  constraint  than  necessary  and  it  does  not  allow  direct  sidelobe  control. 
Recently,  a  generalization  of  this  method  for  large  objects  has  been  presented  [2],  nevertheless  it 
shows  the  same  limitations  as  [1].  Another  solution  in  the  two-dimensional  case,  presented  in  [3], 
tries  to  take  into  account  the  obstacles  as  parasitic  elements,  including  them  in  a  pattern  synthesis 
algorithm,  which  allows  to  take  advantage  of  their  scattering  even  for  large  obstacles.  This 
method  only  considers  scatterers  with  cross  sections  oblong  in  the  direction  of  the  propagation  of 
the  wave,  which  has  been  demonstrated  that  reduce  the  blockage  width  (width  of  an  equivalent 
ideal  shadow  which  produce  the  same  forward  scattered  field  as  the  actual  cylinder)  with  respect 
to  cylinders  with  circular  sections  [4].  Furthermore,  the  improvement  on  the  sidelobe  level 
obtained  in  comparison  to  the  blocked  pattern  is  limited  to  about  8  dB. 

In  this  letter  we  present  a  new  method  based  on  the  use  of  the  simulated  annealing  technique  [5, 
6]  along  with  the  physical  optics  approximation  [7]  to  minimise  iteratively,  an  objective  function, 
which  is  the  difference  between  the  desired  specifications  of  the  problem  and  the  obtained  ones  in 
each  iteration.  This  optimization  is  accomplished  changing  the  amplitude  and/or  the  phase  of 
each  element  excitation.  The  physical  optics  approximation  is  used  to  compute  the  influence  of 
the  scatterer  because,  with  this  approximation,  we  get  acceptable  first-order  results,  only  taking 
minutes  for  each  iteration  of  run  time  on  a  computer  instead  of  hours  that  would  be  necessary 
using  other  more  accurately  techniques  such  as  moment  methods. 
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METHOD 


Let  us  consider  a  linear  array  of  N  equispaced  isotropic  sources,  laid  out  along  the  axis  (that 
leads  to  a  scalar  problem  and  is  representative  of  a  large  aperture  antenna  with  weight  control  in 
one  plane)  and  a  cylinder  with  radius  r  located  in  front  of  it  in  the  (pos^pos^)  position,  as 
indicated  in  Figure  1.  The  resulting  radiation  pattern  taking  into  account  this  scatterer  is  given  by: 


/■(«>)  = 

/=! 


,  Y  7-  JKxcm^os(p+y^sm^) 

m=l 


(1) 


where  x.  is  the  position  of  the  i-th  element  of  the  array,  I.  and  Cir.  the  amplitude  and  phase 
respectively  of  the  i-th  element  excitation  and  is  the  complex  excitation  induced  in  the  point 

m  in  the  cylinder  surface.  In  order  to  simulate  the  scattering  it  was  taken  M  points  from  the 
cylinder  surface  and  the  induced  excitations  in  these  points  were  calculated  using  the  physical 
optics  approximation  by  means  of: 


\/ieU 


(2) 


where  U  is  the  set  of  array  elements  that  “can  see”  the  m  point  in  the  cylinder  surface.  The  term 

r  im  is  the  relative  vector  between  the  i-th  element  in  the  array  and  the  m  point  in  the  cylinder 
surface. 

Once  known  the  power  pattern,  the  objective  function  was  constmcted  through  the  square  of  the 
difference  between  the  desired  sidelobe  level  and  the  obtained  sidelobe  level.  In  the  case  of 
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amplitude  and  phase  control,  it  was  also  included  the  square  of  dynamic  range  ratio  and 

the  square  of  the  difference  between  the  initial  amplitude  and  the  current  amplitude  in  every 
iteration.  Besides,  all  the  terms  are  multiplied  by  weight  constants  adjusted  to  obtain  the  best 
results.  The  minimization  of  this  cost  function  is  accomplished  varying  I.  and  a.  in  the  case  of 
using  amplitude  and  phase  changes  or  a,  in  the  case  of  phase-only  control. 


RESULTS 

We  started  from  a  linear  array  of  41  isotropic  sources  with  d  =  TJl  interelement  spacing.  Using 

initial  excitations  for  this  array  to  produce  a  -40  dB  Taylor  pattern  («=10),  two  examples  were 
performed  using  a  cylinder  with  r  =  1.0  and  taking  M  =  200  points  in  the  cylinder  surface.  In 
the  first  example  the  cylinder  location  was  {OX,  lOA,).  Since  the  cylinder  is  located  in  the  centre  of 
the  linear  array,  symmetric  amplitudes  and  phases  for  the  excitations  in  the  array  were  used.  In 
Figure  2  they  are  shown  the  blocked  pattern  that  has  a  sidelobe  level  of  about  -12  dB,  the 
optimized  pattern  with  phase-only  control  with  a  sidelobe  level  of  -23.4  dB  and  the  optimized 
pattern  with  amplitude  and  the  phase  control,  which  has  a  sidelobe  level  of  -29.7  dB.  In  this  last 
case  the  calculated  dynamic  range  ratio  was  6.94  and,  obviously,  for  the  phase-only  control  case 
the  dynamic  range  ratio  was  the  same  as  the  original  Taylor  pattern  (9.31).  Therefore  depending 
on  the  method  performed,  an  improvement  of  1 1.4  dB  (phase-only  control)  or  17.7  dB  (amplitude 
and  phase  control)  is  obtained. 

In  the  second  example  the  cylinder  is  located  in  (5X,,  lOX,),  therefore  we  do  not  used  symmetric 
amplitudes  and  phases.  In  Figure  3,  there  are  shown  again  the  blocked  pattern,  the  optimized 
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pattern  with  phase-only  control  and  the  one  with  amplitude  and  phase  changes.  The  sidelobe 
levels  in  both  sides  left  and  right  for  each  pattern  are,  respectively:  -15.3  dB  and  -15.1  dB  for  the 
blocked  pattern,  -29.7  dB  and  -24.2  dB  for  the  pattern  obtained  with  phase-only  control 
optimization  and  -34  dB  and  -29  dB  for  the  optimized  pattern  using  amplitude  and  phase  control. 
The  dynamic  range  ratio  in  this  last  case  is  9.3. 


CONCLUSIONS 

We  have  presented  an  array  pattern  synthesis  method  to  take  into  account  the  presence  of  near¬ 
field  scatterers  using  the  simulated  annealing  technique  and  the  physical  optics  approximation  in 
order  to  minimize  an  objective  function.  This  minimization  was  performed  in  two  ways:  changing 
both  amplitude  and  phase  of  the  excitations  of  the  array  or  only  changing  the  excitation  phase. 
This  method  can  be  applicable  to  the  more  complicated  three-dimensional  problem  directly,  in 
which  is  necessary  to  consider  finite  cylinders  and  to  take  into  account  the  polarization.  In  fact, 

dipole  planar  arrays  for  a  ground-mapping  radar  (with  a  esc  {6)cos{d)  vertical  pattern  over  a 
specific  range  and  a  pencil  beam  horizontal  pattern)  in  presence  of  finite  cylinders  is  being 
studied. 

Finally,  it  is  possible  to  apply  the  method  to  scatterers  with  arbitrarily  sections  (not  only  circular). 
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Figure  captions: 


Figure  1 :  Geometry  used  for  the  two-dimensional  problem  of  a  linear  array  antenna  in  presence 
of  a  cylindrical  scatter. 

Figure  2:  Power  patterns  optimized  with  phase-only  control  and  using  amplitude  and  phase 
changes  in  case  of  locating  a  cylinder  in  (OX,  lOX,).  The  initial  blocked  pattern,  in 
dashed  line,  is  also  shown. 

Figure  3:  Power  patterns  optimized  with  phase-only  control  and  using  amplitude  and  phase 
changes  in  case  of  locating  a  cylinder  off-centre  ((5X,  lOX,)  position)  in  relation  to  the 
centre  of  the  array.  The  initial  blocked  pattern  is  shown  in  dashed  line. 
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Abstract: 


A  method  for  the  array  pattern  synthesis  in  the  presence  of  near-zone  scatterers  in  the  three- 
dimensional  vector  case  is  presented.  The  synthesis  is  based  on  the  minimization  of  an 
objective  function  using  the  simulated  annealing  technique.  In  order  to  take  into  account  the 
effects  of  the  scatterers,  the  physical  optics  approximation  is  used.  The  polarization  in  both 
the  incident  and  the  scattered  fields  is  also  taken  into  account.  Examples  for  a  10x1 0-element 
planar  array  of  dipoles  in  the  presence  of  a  perfectly  electric  conducting  cylinder  are 
presented. 
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1.  INTRODUCTION 


The  perturbations  caused  by  a  scatterer  in  the  proximity  of  an  array  antenna  are  being 
a  practical  problem  with  applications  in  ship  or  airborne  antennas.  To  solve  this  problem, 
some  methods  that  deal  with  the  two-dimensional  problem  have  been  developed  in  the 
literature.  In  [1]  an  adapted  pattern  which  avoids  illuminating  near-field  scatterers  on 
transmit,  or  rejects  the  corresponding  scattered  waves  on  receive,  is  proposed  by  introducing 
nulls  in  the  near-zone  of  the  pattern.  However  the  method  may  impose  a  more  severe  field 
constraint  than  necessary  and  does  not  allow  direct  sidelobe  control.  The  method  described  in 
[2]  has  the  same  limitations  as  [1]  but  with  a  generalization  for  large  objects.  The  solution 
presented  in  [3]  uses  a  pattern  synthesis  algorithm  including  parasitic  elements,  trying  to  take 
into  account  their  effects,  which  allows  to  take  advantage  of  the  scattering  even  for  large 
obstacles.  Nevertheless  the  method  only  considers  scatterers  with  oblong  cross  sections, 
which  has  been  demonstrated  that  reduce  the  blockage  width  with  respect  to  scatterers  with 
circular  sections  [4].  These  scatterers  with  circular  sections  were  also  analyzed  in  [5]  using 
the  simulated  annealing  technique  to  avoid  as  much  as  possible  the  presence  of  the  scatterer 
and  taking  a  direct  control  to  the  sidelobe  level  of  the  scattered  power  pattern.  The  problem 
was  also  analyzed  in  two  dimensions  and  the  effects  of  the  scatterer  were  simulated  using 
physical  optics  [6]. 

In  this  paper,  we  present  an  extension  of  [5]  to  three  dimensions,  and  taking  into 
account  the  polarization  to  emulate  a  realistic  situation,  using  again  the  simulated  annealing 
along  with  the  physical  optics  approximation.  This  method  minimizes,  iteratively,  an 
objective  function,  which  involves  parameters  of  interest  in  the  optimization:  sidelobe  levels 
and  current  dynamic  range  ratio.  The  optimization  is  performed  by  changing  the  amplitude 
and/or  the  phase  of  each  element  excitation.  By  using  the  physical  optics  approximation  to 
compute  the  influence  of  the  scatterer,  acceptable  first  order  results  are  obtained  in  only 
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minutes  of  computation  time  instead  of  hours  if  more  accurately  methods,  such  as  moment 
methods,  are  used. 

METHOD 

Let  us  consider  a  planar  array  of  Px  Q  =  N  dipoles  with  21  =  7J2  length  located  in  the 
XZ  plane,  separated  dx.  in  the  x-direction  and  d,  in  the  z-direction.  In  front  of  the  dipole  array, 
a  cylinder  of  radius  r  and  height  h  is  located,  which  acts  as  a  scatterer,  as  it  is  shown  in 
Figure  1. 

The  field  pattern  produced  by  the  array  of  dipoles  in  the  presence  of  the  cylinder  is 
calculated  in  two  steps: 

1.  First  of  all  we  need  to  calculate  the  induced  complex  currents  in  the  cylinder.  We  use 
physical  optics  approximation  over  the  entire  cylinder,  which  has  been  previously  divided 
into  a  certain  number  of  cells  with  area  Aa^,  which  depends  on  the  cylinder  size.  The 
expression  for  every  induced  excitation  over  an  m-th  ceU  is  given  by: 

(1) 

where  is  the  normal  vector  to  the  m-th  ceU  surface,  is  the  complex  magnetic  field 
produced  by  the  i-th  dipole  in  the  m-th  cell  and  V  is  the  set  of  cells  that  are  not  in  the 
“shadow”  of  the  array  radiation.  The  expression  used  for  the  magnetic  field  is  given  by  [7]: 

-  2(cos ]  (2) 

where  all  the  parameters  are  explained  in  Figure  2. 

2.  Once  we  have  calculated  every  induced  current  in  the  cylinder,  we  can  obtain  the  total 
field,  which  is  the  sum  of  the  array  contribution  and  the  cylinder  contribution.  In  order  to 
calculate  the  scattered  field  we  consider  a  local  system  in  every  cell  of  the  cylinder.  The 
contribution  of  every  cell  is  added  (in  a  global  coordinate  system)  to  calculate  the  contribution 
of  the  entire  cylinder. 
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In  Figure  3  it  is  shown  that  exist  two  currents  in  every  cell,  one  in  the  Zm  direction  and 
other  in  the  local  (pm  direction.  However,  due  to  the  physical  optics  approximation  only  a 
current  in  the  Zm  direction  will  exist.  In  any  case,  and  trying  to  simulate  the  real  situation,  the 
existence  of  a  current  in  the  Xm  direction  and  a  current  in  the  z^-direction  is  considered.  The 
presence  of  a  current  in  the  direction  is  supposed  to  be  negligible  because  the  curvature  of 
every  local  cell  is  small  if  we  take  M  large  enough.  To  calculate  the  value  of  each  current  (in 
the  Zm  and  Xm  directions)  we  suppose  that  the  total  current  forms  an  angle  (3  with  the  Zm  axis, 
therefore: 

=4.,.cosi3 

In  other  words  we  assume  that  there  are  two  elementary  electric  dipoles  in  every  cell,  one 


along  the  Zm  direction  and  another  along  the  Xm  direction. 


With  this  assumption,  the  final  vector  field  pattern  is  given  by; 


Fie,<!>)  =  FEig)'^!,e“ 


'kidpXi  sindcoso+dpZfCosS) 


ik{cix„  sine„  cosci„+c(v„  sine„  sino„+c6,cos9„) 


where  the  i-th  dipole  in  the  array  is  located  in  (dpxi,  0,  dpzi)  and  the  m-th  cell  in  the  {clxm, 
clym,  clzni)  position,  the  element  vector  factor  FE{d)  =  FE(d)6  is  the  one  associated  to  a  dipole 
of  21  =  X/2  length  along  the  ^  direction; 


f  TT 

cos  —cosd 
FE(e)  =  — 

sin  6 


and  FE{6.(b)  =  FE^{d^,<l}J6  -f-  FE^{d^,<pJ(p  ,  where: 


®  J  =  T  Se  J6  +  cos -  cos(2^„  +  « J  +  cos(20)] 


FE(e„,(PJ  =  g,ie^,PJ 
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In  the  above  equations,  since  the  dipoles  in  the  array  produces  a  polarization  field  in  the  6 
direction,  the  factors  ge  and  are  the  copolar  and  crosspolar  components  of  the  element 
factors  of  each  cell  in  the  cylinder;  the  angles  Qm  and  (pm  are  the  angles  related  to  the  local 
system  in  the  m-th  cell.  Finally,  the  global  coordinate  system  and  the  local  system  of  a  m-th 
cell  are  related  through  a  clockwise  rotation  about  the  z-axis  of  am.  Therefore,  we  first 
calculate  the  element  factor  in  each  local  system  an  then  we  move  it  to  the  global  system  for 
adding  all  cell  contributions  [8]. 

In  order  to  manage  a  more  suitable  expression  we  separate  the  field  pattern  in  its 
copolar  (d  direction)  and  crosspolar  i<p  direction)  components: 

Fie,(p)  =  FgO  +  Fj  (8) 

In  the  presence  of  the  cylinder  the  initial  power  pattern  is  blocked  and  a  rise  in  the 
sidelobe  level  is  produced.  Our  goal  is  to  obtain  a  suitable  set  of  currents  for  the  array  of 
dipoles  so  that  the  cylinder  contribution  in  the  sidelobe  level  of  the  copolar  coniponent  (Fs) 
was  minimized.  The  crosspolar  component  {F^)  is  also  minimized  as  much  as  possible. 

The  initial  currents  of  the  dipoles  were  perturbed,  through  the  simulated  annealing 
technique,  to  accomplish  the  requirements  about  the  sidelobe  level  and  about  the  current 
dynamic  range  ratio  (Ilma-x/IminO.  which  should  be  maintained  close  (or  reduced)  to  low  values. 
The  perturbation  can  be  performed  varying  the  amplitude  and  phase  or  varying  only  the  phase 
of  each  dipole  excitation.  The  choice  of  the  better  amplitudes  and/or  phases  is  made  according 
to  the  minimization  of  a  given  cost  function,  which  involves  the  parameters  that  need  to  be 
controlled  (i.e.  sidelobe  levels  and  dynamic  range  ratio). 

EXAMPLES 

In  this  section  we  present  some  examples  selected  from  all  the  examples  performed. 
We  chose  a  dipole  array  of  10  x  10  elements  separated  0.7A,  in  both  directions  (z  and  x).  The  [3 
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angle  used  was  P  =  10  deg.  In  all  of  them  the  initial  power  pattern  used  is  a  Chebyshev  power 
pattern  (without  cylinder)  with  a  sidelobe  level  of  -40  dB  and  with  a  current  dynamic  range 
ratio  of  31.1.  Two  different  cylinders  with  0.51  radius  were  used;  one  was  61  high  and  the 
other  41  high.  Each  cylinder  was  located  centered  in  front  of  the  array  of  dipoles  at  distances 
of  51  and  81.  The  cylinders  were  divided  in  the  minimum  number  of  cells  as  possible, 
maintaining  the  shape  of  the  total  power  pattern,  i.e.,  more  cells  in  the  cylinder  did  not  alter 
the  pattern  in  a  significant  way. 

Example  1:  Cylinder  height  h  =  6X,  distance  from  the  array  d  =  8X.  Number  of  cells  in  the 
cylinder:  25  y.  15  -  375. 

After  locating  the  cylinder,  the  highest  sidelobe  level  for  the  copolar  component  was 
-25.19  dB,  and  the  crosspolar  -39.23  dB,  which  occurs  in  the  (j)  =  0  cut. 

1  .a)  Amplitude  and  phase  control. 

Performing  the  optimization  with  amplitude  and  phase  control,  the  final  sidelobe  level 
was  -30.31  dB,  which  implies  an  improvement  of  5.12  dB.  The  crosspolar  remained  almost 
unaltered  and  its  final  value  was  -38.05  dB.  In  case  of  the  current  dynamic  range  ratio  it  was 
lowered  to  16.97.  In  Figure  4  the  blocked  and  the  recovered  power  patterns  (copolar  and 
crosspolar  components)  in  the  (b  =  0  cut,  which  was  always  the  most  perturbed,  are  shown. 

In  order  to  have  a  more  general  view  of  the  sidelobe  level  improvement,  in  Figures  5  and  6  a 
3D  representation  of  the  blocked  and  the  optimized  copolar  components  of  the  power  patterns 
are  shown.  The  Figure’s  floor  is  set  to  -32  dB  and  in  both  of  them  u  =  sin  0  cos  b  and 
w  =  cos  0,  since  the  array  of  dipoles  is  located  in  the  XZ  plane.  The  Figure’s  floor  and  u  and  v 
will  be  the  same  in  the  rest  of  the  paper. 

Lb)  Phase-only  control. 

In  this  case  only  a  perturbation  in  the  phases  of  the  dipole  excitations  was  introduced. 
The  optimized  pattern  following  this  procedure  had  a  sidelobe  level  in  the  copolar  component 
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of -29.61  dB  (4.37  dB  less  than  the  blocked  pattern)  and  the  crosspolar  one  has  a  level  of 
-37.49  dB.  The  current  dynamic  range  ratio,  obviously,  remained  unaltered  to  31.1. 

The  (!)  =  0  cut  of  the  blocked  power  pattern  corresponding  to  the  copolar  component  and  the 
recovered  power  patterns  (copolar  and  crosspolar  components)  are  shown  in  Figure  7. 
Figure  8  shows  a  3D  plot  of  the  copolar  component  of  the  recovered  power  pattern. 

Example  2:  Cylinder  height  h  =  4X  distance  from  the  array  d  =  5A.  Number  of  cells  in  the 
cylinder:  20  x  10  =  200. 

In  this  case,  the  initial  copolar  blocked  component  in  the  presence  of  the  cylinder  had  a 
sidelobe  level  of -26.92  dB  and  the  crosspolar  of  -40.60  dB. 

2. a)  Amplitude  and  phase  control. 

The  optimization  using  amplitude  and  phase  perturbations  yielded  a  final  power 
pattern  with  a  sidelobe  level  of  -30.90  dB  for  the  copolar  component  (4.02  below  the  blocked 
pattern).  The  crosspolar  component  remained  in  -40.05  dB  and  the  current  dynamic  range 
ratio  for  the  final  pattern  was  22.34. 

As  in  the  previous  example,  Figure  9  shows  the  (p  =  0  cut  of  the  copolar  component  of  the 
blocked  and  recovered  power  patterns,  and  the  crosspolar  component  of  the  recovered  power 
pattern.  A  3D  representation  of  the  copolar  component  of  the  blocked  and  recovered  power 
patterns  are  shown,  respectively,  in  Figures  10  and  11. 

2.b)  Phase-only  control. 

The  searching  for  the  best  solution  using  phase-only  control  gave  an  optimized  pattern 
with  a  sidelobe  level  of  -28.07  dB  (1.15  dB  less  than  the  blocked  patter)  for  the  copolar 
component  and  -40.03  for  the  crosspolar  component.  The  initial  copolar  component  of  the 
blocked  pattern  and  the  copolar  and  crosspolar  components  of  the  optimized  power  patterns 
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are  shown,  for  the  ({)  =  0  cut,  in  Figure  12.  Finally,  Figure  13  shows  a  3D  representation  of  the 
copolar  component  of  the  recovered  pattern. 

Example  3:  Cylinder  height  h  =  6X,  distance  from  the  array  d  =  5X  Number  of  cells  in  the 
cylinder:  25  x  15  =  375. 

In  this  example  we  tried  to  test  the  importance  of  the  cylinder  height  locating  the  same 
cylinder  as  in  the  first  example  at  a  distance  d  =  5X  (the  distance  used  in  the  second  example). 
The  results  are  shown  in  Table  1. 

Finally,  examples  1  and  2  were  performed  again  but  using  |3  =  0  deg.,  i.e.,  the  current 
in  every  cylinder  cell  had  z  direction,  thus  the  power  pattern  has  only  a  copolar  component. 
The  results  are  summarized  in  Table  2. 

CONCLUSIONS 

A  very  important  factor  in  the  optimization  is  the  scatterer’s  height  (with  reasonable 
diameters).  The  distance  is  also  important  (large  distances  imply  better  results)  but  with  two 
scatterers  with  different  heights  located  at  the  same  distance  from  the  array,  we  obtained 
better  results  for  the  higher  scatterer.  Examples  with  greater  diameters  were  performed 
obtaining  similar  results  as  the  presented  ones,  confirming  the  theory  of  the  height 
importance. 

With  regard  to  the  optimization  methods,  with  amplitude  and  phase  control  we  got  results,  in¬ 
most  cases,  only  a  few  dB’s  better  than  using  phase-only  control.  This  implies  that  only 
changing  the  phase  of  the  currents  in  the  dipole  array,  which  is  physically  easier  than  change 
the  amplitude  and  the  phase,  we  can  obtain  good  results  that  improve  the  scattered  power 
pattern. 
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In  relation  to  the  (3  angle,  with  small  values  the  optimization  is  satisfactory  and  similar  to  the 
case  with  [3  =  0  deg.  (only  a  copolar  component).  As  the  value  of  P  increases  the  results  are 
worse,  since  there  are  more  difference  between  the  copolar  and  crosspolar  components  of  the 
scattered  field. 
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Figure  Captions; 

Figure  1 .  Array  of  dipoles  in  presence  of  a  scatterer  cylinder. 

Figure  2.  Situation  of  every  dipole  in  the  array  and  parameters  used  in  equation  (2). 

Figure  3.  Current  distribution  in  an  m-th  cell  of  the  cylinder. 

Figure  4.  ({)  =  0  cut  of  the  optimized  (crosspolar  and  copolar  components)  and  initial 

(copolar  component)  power  patterns  using  amplitude  and  phase  control,  h  =  6X, 
d  =  8X. 

Figure  5.  Surface  plot  of  the  initial  copolar  component  of  the  power  pattern,  h  =  6X,  d  =  8X. 

Figure  6.  Copolar  component  of  the  power  pattern  after  performing  the  optimization  using 
amplitude  and  phase  control  h  =  6X  d  =  8X. 

Figure  7.  (j)  =  0  cut  of  the  optimized  (crosspolar  and  copolar  components)  and  initial 

(copolar  component)  power  patterns  using  phase-only  control  h  =  6X,  d  =  8X. 

Figure  8.  Surface  plot  of  the  recovered  copolar  component  of  power  pattern  using  phase- 
only  control,  h  =  6X,  d  =  8X. 
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Figure  9. 


Figure  10. 
Figure  11. 


Figure  12. 


Figure  13. 


(t)  =  0  cut  of  the  optimized  (crosspolar  and  copolar  components)  and  initial 
(copolar  component)  power  patterns  using  amplitude  and  phase  control  h  =  4X, 
d  =  5X. 

Surface  plot  of  the  initial  copolar  component  of  the  power  pattern,  h  =  4X,  d  =  5X. 
Copolar  component  of  the  power  pattern  after  performing  the  optimization  using 
amplitude  and  phase  control,  h  =  4X,  d  =  5X. 

({)  =  0  cut  of  the  optimized  (crosspolar  and  copolar  components)  and  initial 
(copolar  component)  power  patterns  using  phase-only  control,  h  =  4X,  d  =  5X. 
Surface  plot  of  the  recovered  copolar  component  of  power  pattern  using  phase- 
only  control,  h  =  4X,  d  =  5X. 


Table  captions: 

Table  1.  Results  obtained  for  a  cylinder  with  h  =  6X,  d  =  5X  and  =10  deg. 
Table  2.  Summary  of  results  using  jS  =  0  deg. 
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Figure  7 
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Figure  12. 


Copolar  component 

Crosspolar  component 

Optimization 

method 
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(dB) 
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(dB) 

Initial  SLL 

(dB) 

Final  SLL 

(dB) 

Ilmax/Iminl 

Amplitude 

and  phase 

control 

-25.35 

-30.45 

5.10 

-38.97 

-38.30 

15.90 

Phase-only 

control 

-25.35 

-28.97 

3.42 

-38.97 

-38.38 

31.10 

Table  1. 


Optimization 

Method 

h 

(X) 

d 

(^•) 

Initial 

sidelobe 

level 
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level 

Improvement 

(dB) 

IlmaxTIminl 

Amplitude  and 

Phase  Control 

6 

8 

-25.08 

-30.19 

5.11 

18.15 

Phase-Only 

Control 

6 

8 

-25.08 

-29.53 

4.45 

31.10 

Amplitude  and 

Phase  Control 

4 

5 

-26.96 

-30.96 

4.00 

31.67 

Phase-Only 

Control 

4 

5 

-26.96 

-28.39 

1.43 

31.10 

Table  2. 
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where  n^is  the  normal  vector  to  the  m-th  surface,  is  the  complex  magnetic  field  produced  by  the  i-th  dipole 
in  the  m-th  cell  and  U  is  the  set  of  cells  that  are  not  in  the  "shadow”  of  the  array  radiation. 

2.-  Once  we  have  calculated  every  induced  current  in  the  cylinder,  we  can  obtain  the  final  field  pattern,  which  is 
the  sum  of  the  array  contribution  and  the  cylinder  contribution: 

F(0  ({))=F£:(0)^  ^^;i:(rfpx5inecos<t)-4jzcose)^^  ^  ^y*(c/i^sinecos<t)-cty„sinesin(t)+dz„cos0) 

i-l  m-I 

where  the  element  vector  factor  FE(6)&  is  the  one  associated  to  a  dipole  of  21  =  A/2  length  along  the  z  direction 
and  F£(0^,(j)J=F£Q(0^,(|)J0+FE^(0^,(j)j4>.  where  FEq  and  FE^  are,  respectively,  the  copolar  and  crosspolar 
components  of  the  vector  element  factor  for  each  cell  in  the  cylinder.  The  i-th  dipole  in  the  array  is  located  in 
{dpxi,  dpyi,  dpzi)  position  and  the  m-th  cell  in  the  iclx„,  cly„,  clzj  position. 

Using  the  physical  optics,  the  induced  current  in  every  cell  has  only  a  z  component.  Therefore,  to  simulate 
the  actual  presence  of  the  copolar  and  crosspolar  components  of  the  scatterer  field  we  make  the  assumption  that 
in  every  cell  the  current  forms  an  angle  (P)  with  the  z„  axis  which  is  the  z  axis  of  the  local  coordinate  system  of 
the  m-th  cell.  Therefore,  there  exist  two  components  of  the  current,  one  in  the  Zm  direction  and  the  other  in  the 
x„  direction,  i.e,  the  field  associated  to  every  cell  is  supposed  to  be  the  field  produced  by  two  elementary  electric 
dipoles,  one  along  the  local  x„  axis  and  another  along  the  local  z„  axis.  We  do  not  include  a  dipole  in  the  local  y„ 
direction  because  the  curvature  of  each  cell  is  small  enough  to  consider  a  plannar  cell. 

To  manage  a  more  suitable  expression  of  the  total  field,  we  divide  it  in  its  copolar  and  crosspolar 

components,  so  that  the  total  field  is:  F(0,<j))=FgS+F^$.  In  the  presence  of  the  cylinder,  the  initial  power  pattern 
is  blocked  and  a  rise  in  the  sidelobe  level  is  produced.  Our  goal  is  to  obtain  a  suitable  set  of  currents  for  the  array 
of  dipoles  so  that  the  cylinder  contribution  in  the  sidelobe  level  of  the  final  pattern  copolar  component  (F^)  is 
minimised.  The  crosspolar  component  (F^)  is  also  minimised  as  imuch  as  possible. 

The  initial  excitation  currents  of  the  dipoles  are  perturbed,  through  the  simulated  annealing  technique, 
to  accomplish  the  requirements  about  the  sidelobe  level  and  even  about  the  current  dynamic  range  ratio  (Ilmax/ImmD. 
which  should  be  maintained  (or  reduced  to)  low  values.  This  perturbation  can  be  performed  by  varying  the 
amplitude  and  phase  or  varying  only  the  phase  of  each  dipole  excitation.  The  choice  of  better  amplitudes  and/or 
phases  is  made  according  to  the  minimisation  of  a  given  cost  function,  which  involves  the  parameters  that  need 
to  be  controlled  (i.e.  sidelobe  levels  and  current  dynamic  range  ratio). 

EXAMPLE: 

An  example  selected  between  all  examples  performed  is  presented  in  this  section.  We  chose  a  dipole  array 
of  10  X  10  elements  separated  0.7  X  in  both  directions  (z  and  x).  The  angle  |3  was  set  to  10°.  The  initial  power 
pattern  used  is  a  Chebyshev  power  pattern  (without  cylinder)  with  a  sidelobe  level  of  -40dB  and  a  current  dynamic 
range  ratio  of  3 1 . 1 .  As  scatterer,  we  use  a  cylinder  with  height  h  =  6^,  lotaced  at  a  distance  firom  the  array  d  =  82., 
and  the  number  of  cells  used  in  the  cylinder  was  25  x  15  =  375.  After  locating  the  cylinder,  the  highest  sidelobe 
level  for  the  copolar  component  is  -25.19  dB,  and  the  crosspolar  -39.23  dB,  which  occurs  in  the  (J)  =  0  cut. 

Performing  the  optimisation  with  amplitude  and  phase  control,  the  final  sidelobe  level  is  -30.3  IdB,  which 
implies  an  improvement  of  5. 12  dB.  The  crosspolar  component  remained  almost  unaltered  and  its  final  value  was 
-38.05  dB.  The  dynamic  range  ratio  was  lowered  to  16.97. 

In  Figure  2  the  blocked  and  the  recovered  power  patterns  (copolar  and  crosspolar)  in  the  4)  =  0  cut  (the  most 
perturbed)  are  shown.  In  Figures  3  and  4  we  show  a  3D  representation  of  the  initial  and  recovered  copolar 
components  of  the  power  patterns.  Using  phase-only  control  the  final  sidelobe  level  was  -29.61dB  (4.37  dB  less 
than  the  blocked  pattern)  for  the  copolar  component  and  the  crosspolar  component  had  a  level  of -37.39  dB. 


AMPLITUDE  AND  PHASE  CONTROL 


•  eicckwd  p«iMn 

R«eak*r*d  p«t*m 
-I CroMpdw 


— ^ - 


Theta  (deg) 


Figure  2.  (}>  =  0  cut. 


Figure  3.  Initial  Power  Pattern 


Figure  4.  Optimised  Power  Pattern 


CONCLUS  ^  Qf  all  an  important  factor  in  the  optimisation  is,  with  reasonable  diameters,  the  scatterer  s  height. 
The  distance  is  also  important  (larger  distances  imply  better  results)  but  with  two  scatterers  with  different  heights 
Sated  at  the  same  disfance  from  the  array,  we  obtain  better  results  for  the  higher  scatterer.  In  addition  several 
examples  using  scatterers  with  the  same  height  and  different  diameters  gave  sum  ar  results.  Besides,  with 
amplitude  and  phase  control  we  get  results  only  about  IdB  better  than  using  phase-only  control. 

ReLding  the  p  angle,  with  small  values  the  optimisation  is  satisfactory  and  similar  to  the  case  with  p  -  0,  which 
the  exiLnce  of  only  a  copolar  component.  As  the  value  of  P  increases  the  results  are  worse,  since  there 
are  more  difference  between  the  copolar  and  crosspolar  components  of  the  scattered  field.  Finally,  in  e 
symposium  we  will  also  shown  results  using  this  technique  but  minimising  the  field  in  every  cell,  i.e.,  minimising 

the  induced  currents. 
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Phase-only  synthesis  of  continuous  linear 
aperture  distribution  patterns  with 
asymmetric  side  lobes 

A.  Trastoy  and  F.  Ares 

A  method  for  synthesising  linear  aperture  distribution  patterns 
with  different  side  lobe  levels  (SLLs)  on  both  sides  of  the  mam 
beam  is  proposed.  This  phase-only  control  method  was  apphrf  to 
a  Taylor  distribution  as  well  as  to  a  uniform  Mphtude 
distribution  and  achieved  differences  in  the  SLL  on  both  sides  of 
the  main  beam  as  great  as  23dB.  improving  on  a  previous 
procedure  that  had  achieved  a  gap  of  little  more  than  lOdB. 

Introduction:  Radar  and  communication  applications  ^st  in 
which  asymmetric  side  lobes  provide  a  system  advantage,  ^ides, 
phase  control  is  less  complicated  and  far  cheaper  to  pro^de  than 
amplitude  control.  Focusing  on  these  ideas,  Frey  and  Elliott  [Ij 
froze  the  amplitude  distribution  and  the  shape  of  the  phase  distri¬ 
bution  of  a  modified  Taylor  pattern.  They  calculated  an  analytic 
function  that  fitted  the  shape  of  that  phase  distribution  and,  by 
controlling  the  magnitude  of  the  swing,  they  obtained  differentials 
as  great  as  10  dB  tetween  the  two  side  lobe  levels  (SLLs)  on  each 
side  of  the  main  beam.  In  this  Letter,  a  method  based  on  the 
numerical  calculation  of  the  odd  phase  distribution  nece^ary  to 
achieve  an  asymmetric  pattern  with  a  greater  difference  between 
the  levels  on  the  two  sides  of  the  main  beam,  through  vanation  of 
the  roots  in  the  modified  Taylor  distribution,  is  described  and 
applied  to  a  Taylor  distribution  as  well  as  to  the  important  case  of 
uniform  amplitude  distribution. 

Description  of  method:  Given  a  linear  Taylor  or  uniform  amplitude 
distribution,  we  freeze  the  amplitude  of  the  aperture  distnbuuon. 
We  then  slightly  modify  the  associated  roots  u„  of  a  modified  Tay¬ 
lor  pattern  developed  by  Elliott  [2]; 


5o(u)  = 


nn— 1 

n  (1  -  U/Un) 

sin TTtr  n=-(n£,-i) _ _ 

TCU  /I  /  \ 

n  (1  -  u/n) 

n=-(ni.-l) 


and  use  the  following  formula  to  calculate  their  associated  aper¬ 
ture  distribution; 


1 

KO  =  —.  E  ^o(m)e-^' 


Tn=  — (fii,— l) 

where  m  is  an  integer,  2a  is  the  length  of  the  aperture,  and  and 
rtf.  are  positive  integers  that  denote  the  transition  roots  on  the 
right  and  left  sides  of  the  main  beam,  respectively. 


-10  -8  -6  -4  -2  0  2  4 

u=(2a/X)  cos  0 

Fig.  1  Obtained  pattern  corresponding  to  -25dB,  h 
tude  distribution 


7  Taylor  ampli- 


The  frozen  amplitude  distribution  and  the  phase  distribution 
(odd  about  the  centre)  obtained  from  eqn.  2  produce  an  asymrnet- 
ric  pattern.  The  simulated  annealing  technique  [3]  is  used  to  find 
the  location  of  the  roots  whose  phase  distribution  provides  the 
highest  difference  between  SLLs  on  both  sides  of  the  mam  beam 


through  the  minimisation  of  an  objective  function  that  takes  into 
account  the  SLLs  on  both  sides  of  the  main  beam. 


-5  ^  -3  -2  -1 


Fio.  2  Phase  distribution  of  pattern  shown  in  Fig.  I 
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U  =  (2a/>.)  cos  0  2553 

Fig.  3  Obtained  pattern  corresponding  to  uniform  amplitude  distribution 


Fig.  4  Phase  distribution  of  pattern  shown  in  Fig.  3 

Results:  To  illustrate  the  method  described  above,  two  examples 
are  presented.  In  both  cases,  the  pattern  is  produced  by  an  aper¬ 
ture  of  length  2a  =  10  L  and  only  the  SLL  of  the  lower  side  of  the 
main  beam  was  included  in  the  objective  function,  with  the  other 
side  being  uncontrolled.  The  first  example  starts  from  a  linear 
Taylor  distribution  calculated  for  an  SLL  of  -25  dB  and  n  =  7. 
We  set  =  fit  =  7  in  eqns.  1  and  2  and  we  obtain  the  phase  dis¬ 
tribution  shown  in  Fig.  2  that  corresponds  to  the  pattern  shown  in 
Fig.  1  and  presents  a  difference  between  SLLs  of  19dB.  It  can  be 
seen  that,  by  using  the  actual  phase  distribution  instead  of  that 
obtained  by  the  fitness  function  described  in  [1],  an  improvement 
in  the  gap  of  ~9dB  is  obtained.  In  the  second  example  (Fig.  3), 
starting  from  a  uniform  amplitude  distribution,  we  obtain  a  gap  of 
23  dB  between  both  sides  of  the  main  beam  (this  time  an  improve¬ 
ment  of  13dB  over  the  results  reported  in  [1]  has  b^n  achieved). 
Its  phase  distribution  is  shown  in  Fig.  4.  Note  that  m  both  exam- 
pies  the  main  beam  peak  is  shifted  from  u  —  0.  This  effect  can  be 
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Mkjh)  -  y  ]n'‘hk  k  =  0, 1 . p  -  1 


(1) 


I  .  u  1 

gressive  phase  to  the  aperture  distribution. 

Conclusions:  A  method  for  synthesising  as3mmetric  patterns  with 
phase-only  control  has  been  presented.  Previous  results  obtained 
by  Frey  and  Elliott  [1]  have  been  improved  and  patterns  with  gaps 
of  -20dB  achieved.  The  computational  time  for  the  two  examples 
shown  above  was  -lOmin  using  a  Pentium  200MHz  processor.  It 
is  possible  to  e.xtend  this  method  in  order  to  fix  nulls  with  phase- 
only  control  and  also  to  lower  the  whole  side  lobe  level  of  both 
linear  and  circular  apertures  using,  in  this  case,  an  even  phase  shift 
about  the  centre  for  linear  apertures  and  a  ((^symmetric  phase  dis¬ 
tribution  for  circular  apertures,  reported  in  [4]  for  patterns  with 
null  filling. 
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Property  of  vanishing  moments  of 
orthogonal  M-band  compactly  supported 
interpolating  scaling  function 

Zhang  Jiangkang  and  Bao  Zheng 

It  is  shown  that  the  number  of  vanishing  moments  of  the  even 
number-order  orthogonal  Af-band  compactly  supported 
interpolating  scaling  function  increases  by  one.  Consequently,  the 
asymptotic  performance  of  the  Mallat  projection  is  the  same  as 
the  orthogonal  projection  for  the  lo-approximation  of  smooth 
functions. 

Introduction:  The  main  advantages  of  interpolating  wavelet  trans¬ 
forms  are: 

(i)  The  wavelet  series  transform  of  a  signal  in  the  multiresolution 
subspace  is  exactly  consistent  with  its  discrete  wavelet  transform 
[1]  when  the  synthesis  scaling  fimction  has  the  interpolating  prop¬ 
erty,  i.e.  (p(n)  =  8{n),  where  5(«)  is  the  Kronecker  sequence. 

(ii)  The  first  N  moments  of  the  iV-order  interpolationg  scaling 
function  are  vanishing  except  for  the  first  moment.  However,  in 
the  two-band  case,  orthogonal  wavelets  do  not  exist  which  have 
both  the  interpolating  property  and  compact  support  except  for 
the  Haar  wavelet. 

Therefore,  in  [2,  3],  we  constructed  high-order,  orthogonal,  M- 
band  and  compactly  supported  interpolating  scaling  functions.  In 
this  Letter  we  reveal  furthermore  that  the  orthogonality  brings  the 
interpolating  scaling  function  an  extra  vanishing  moment  prop¬ 
erty. 

Definition  of  moments:  The  p-order  moments  of  a  scaling  filter  and 
the  scaling  function  are  defined  as,  respectively, 
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f  t%-{t)dt 

J  — oo 


A:  =  0,1 . p-1  (2) 


It  is  well-known  that  if  (p(0  is  an  A-order  interpolating  scaling 
function,  then  it  has  the  Coiflet  property  [4]  or  vanishing  moment 
property,  i.e. 


.\h{h)  =  A4(p)  =  S{k)  k  =  0, 1 . .V  -  1  (3) 


In  general  cases.  Mfh)  ^  0  and  iV/dip)  it  0.  However,  in  the  case 
of  an  even  number-order,  orthogonal  and  jV/-band  interpolating 
scaling  function,  the  orthogonality  forces  the  number  of  vanishing 
moments  to  increase  by  one. 


Property:  Throughout  this  Letter,  we  assume  that  the  interpolat¬ 
ing  scaling  function  (p(/)  satisfies 


L 

Pit)  =  \/A7^/n.<p(Mt-A-) 

Jfc=0 

and  that  the  scaling  filter  h  satisfies 


H{=) 


1  +  Z  +  ---  +  Z 

M 


.\t-i 


Q(-) 


(4) 

(5) 


and 


H{1)  =  v/M  (6) 


Funhermore,  (p(f)  is  assumed  to  satisfy  the  orthogonal  condition 


k 


An  immediate  consequence  of  eqn.  7  is 


A/-1 

E 

r=0 


H 


=  M 


(7) 

(8) 


Now  we  claim  that  the  moments  of  (p(/)  have  the  following 
property: 


Theorem-.  If  the  interpolating  scaling  filter  h  and  the  corresponding 
interpolating  scaling  function  <p{t)  satisfy  eqns.  4-8,  then, 

Mkih)  =  6{k)  k  =  0, 1, ...,  N  -  1  (9a) 

Mkiv>)  =  S{k)  k  =  0, 1, ..., N  -  1  (96) 

In  particular,  if  iV  is  an  even  number,  then 

Mk {h)  =  0  k  =  N,  N  +  2,  -V  -l-  4, ...,  2.V  -  2  (lOo) 

A4((p)  =  0  A:  =  iV,A--l-2,A'-F4,...,’2i\'-2  (106) 


Remarks-.  By  using  the  theorem  and  Unser’s  methods  [5],  it  can  be 
shown  that  the  Mallat  projection  of  ft)  onto  the  multiresolution 
subspace  F/ip) 

jfe  '  ' 

has  the  same  asymptotic  performance  as  the  corresponding 
orthogonal  projection 

Pjf  =  ^{f{t),V>J.k)v>j,k{t) 

k 

i.e.  both  M/f)  and  P/are  /^-approximations  of  the  smooth  func¬ 
tion^/)  as  the  scale  approaches  zero.  In  this  sense,  the  orthogonal 
A/-band  compactly  supported  interpolating  wavelet  compares 
favourably  with  the  two-band  compactly,  supported  interpolating 
wavelet. 


Proof  of  theorem:  We  know  that  the  scaling  filter  and  its  poly¬ 
phases  satisfy  the  relationship 

M-l 

H{z)  =  Y,  (11) 

r=0 
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mented  using  lossless  orthogonal  beam-forming  networks  such  as  the 
Butler  matrix  [2].  The  major  objections  to  the  Woodward-Lawson 
method  have  been  its  inability  to  control  either  the  amplitude  of  the 
ripple  in  the  shaped  region  of  the  pattern,  or  the  heights  of  side  lobes 
in  the  unshaped  region  [2,  3].  Another  telling  objection  is  that  the  ex¬ 
citations  it  prescribes  often  have  very  high  dynamic  range  ratios,  and 
may  not  be  physically  realizable  in  the  presence  of  mutual  coupling 
[4-6].  These  weaknesses  are  ultimately  due  partly  to  the  fact  that 
the  Woodward-Lawson  method  constrains  all  the  beams  with  which  it 
synthesizes  the  desired  pattern  to  have  the  same  phase,  so  that  the  far 
field  pattern  is  real,  and  partly  to  its  zero  tolerance  of  deviation  from 
the  desired  field  at  the  points  at  which  it  samples  this  field.  When  the 
component  beams  can  differ  in  phase  (because  a  real  far  field  is  not 
necessary)  and/or  perturbation  of  sample  amplitudes  is  allowed,  the 
extra  degrees  of  freedom  that  thereby  become  available  should  in  prin¬ 
ciple  allow  calculation  of  excitation  distributions  that  are  more  easily 
realized  and/or  afford  better  side  lobe  levels  and  ripple  characteristics. 
In  this  paper  we  describe  how  this  may  be  achieved  using  simulated 
annealing  optimization  methods  [7].  We  know  of  no  previous  attempt 
to  improve  the  Woodward-Lawson  method  in  this  way. 

2.  DESCRIPTION  OF  THE  METHOD 

Briefly,  the  Woodward-Lawson  method  for  a  linear  array  of  N  antenna 
elements  spaced  a  distance  d  apart  consists  in  first  sampling  a  desired 
field  pattern  F{9)  in  2M-M  directions  9i  spaced  cos“^(A/JVd)  apart 
(where  9  is  measured  from  endfire  and  M  is  such  that  the  visible 
region  is  just  covered),  and  then  approximating  F{9)  by 

=  (1) 

i^—M 

where  Ui  =  {ixd/ \{cos{9)  —  cos(0i)) .  Each  of  the  component  quasi- 
sinc  beams  sm{N Ui) !  5m[ui)  is  realized  by  an  excitation  distribution 
with  uniform  amplitude  and  uniform  phase  differences  between  the  el¬ 
ements,  and  WL{9)  is  produced  by  the  weighted  sum  of  these  distri¬ 
butions,  with  weights  F{9i)jN .  WL{9)  reproduces  F{9)  exactly  at 
each  sample  point  because  all  the  component  quasi-sinc  beams  except 
the  i-th  are  zero  at  9i . 


Shaped  power  patterns  produced  by  equispaced  linear  arrays 
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Figure  1.  Flat-topped  power  patterns  obtained  by  optimizing  the  am¬ 
plitudes  of  a  linear  combination  of  sin(iVa;)/sin(x)  beams  (real  field) 
and  then  optimizing  their  phases  (complex  field). 

As  noted  above,  the  original  Woodward-Lawson  method  assumes 
that  all  the  F{di)  are  in  phase,  and  the  coefficient  of  the  i-th  quasi- 
sinc  beam  is  exactly  F{9i) .  In  this  work  we  perturbed  the  amplitude 

of  each  F{9i)  by  6i  and/or  multiplied  all  but  one  by  a  phase  factor 
■ 

(2) 

We  used  simulated  annealing  to  determine  the  set  of  4>i  and/or  5i 
minimizing  a  cost  function  that  penalized  both  the  deviation  of  F[6) 
from  desired  pattern  characteristics  (side  lobe  levels,  ripple  levels  and 
beamwidth)  and  the  dynamic  range  ratio  of  the  excitation  distribution. 

3.  EXAMPLES 

A)  Suppose  that  we  require  a  20-element  (A/2)-spaced  linear  array 
to  produce  a  fiat-topped  power  pattern  with  a  —3  dB  beamwidth 
of  cos"^(0.4)  and  side  lobe  levels  of  — 40dB  for  the  two  inner 
lobes  and  — 20dB  for  the  others.  Fig.  1  (“real  field”  pattern) 
shows  that  the  specified  side  lobe  levels  and  ripple  levels  of  about 
±0.01  dB  can  be  achieved  by  optimizing  just  the  amplitudes  of 
the  beam  coefficients;  the  optimal  amplitudes  are  listed  in  Table 
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1.  If  the  amplitudes  are  now  fixed  at  these  values  and  the  phases 
of  the  beam  coefficients  are  optimized  (the  optimal  phases  are 
also  listed  in  Table  1),  there  is  virtually  no  change  in  the  power 
pattern  (Fig.  1)  but  the  dynamic  range  ratio  of  the  excitation 
distribution  falls  from  17.09  to  6.83,  an  improvement  of  60%. 
Fig.  2  shows  the  phase  pattern  of  the  resulting  complex  field. 


Table  1.  Beam  coefficients  corresponding  to  the  flat-topped  power 
patterns  shown  in  Fig.  1  (the  coefficients  for  i  =  ±10  are  not  shown 
because  for  a  symmetric  pattern  these  beams  cancel  each  other). 


|B 

Amplitudes 

Phases  (deg) 

Real  field 

Complex  field 

0 

1.134 

0 

-70.124 

±1 

1.130 

0 

-43.706 

±2 

0.822 

0 

-5.389 

±3 

0.078 

0 

19.760 

±4 

0.011 

0 

20.807 

±5 

0.115 

180 

201.042 

±6 

0.099 

0 

29.516 

±7 

0.076 

180 

215.331 

±8 

0.051 

0 

40.600 

±9 

0.026 

180 

215.337 

B)  Using  the  same  20-element  (A/2)-spaced  linear  array  as  above, 
we  desire  to  synthesize  a  pattern  with  a  cosecant  squared  beam 
in  the  sector  cos^  £  [0.1  — 0.5],  four  — 30dB  side  lobes  to  the  left 
of  the  cosecant  beam  and  no  other  side  lobe  higher  than  —  20  dB. 
Fig.  3  (“real  field  pattern”)  shows  the  result  of  optimizing  just  the 
amplitudes  of  the  quasi-sinc  beam  coefficients:  maximum  ripple 
in  the  cosecant  region  is  ±0.6  dB.  The  optimal  amplitudes  are 
listed  in  Table  2.  As  before,  if  these  amplitudes  are  fixed  and 
the  phases  of  the  beam  coefficients  are  optimized  (the  optimal 
phases  are  also  listed  in  Table  2),  there  is  no  significant  change 
in  the  power  pattern  (Fig.  3)  but  the  dynamic  range  ratio  of  the 
excitation  distribution  falls  from  12.05  to  8.22,  an  improvement 
of  32%.  Fig.  4  shows  the  phase  pattern  of  the  resulting  complex 
field. 


Shaped  power  patterns  produced  by  equispaced  linecir  cirrays 
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Figure  2.  Phase  pattern  of  the  complex  field  of  Fig.  1. 


W=COS(0) 


Figure  3.  Cosecant  squared  patterns  obtained  by  optimizing  the  am¬ 
plitudes  of  a  linear  combination  of  sin(Na:) /  sin(x)  beams  (real  field) 
and  then  optimizing  their  phases  (complex  field). 
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Table  2.  Beam  coefficients  corresponding  to  the  cosecant  squared 
patterns  shown  in  Fig.  3. 


i 

Amphtudes 

Phases  (deg) 

Real  field 

Complex  field 

-10 

0.074 

180 

195.126 

-9 

0 

19.304 

-8 

0.022 

0 

20.278 

-7 

0.078 

180 

192.667 

-6 

0.114 

0 

14.165 

-5 

0.131 

180 

194.168 

-4 

0.038 

0 

0.411 

-3 

0.042 

180 

179.153 

-2 

0.040 

0 

-3.660 

-1 

0.002 

0 

28.173 

0 

0.927 

0 

-3.324 

1 

1.224 

0 

1.562 

2 

0.651 

0 

8.696 

3 

0.443 

0 

-1.586 

4 

0.326 

0 

12.995 

5 

0.261 

0 

11.997 

6 

0 

-13.026 

7 

0.123 

180 

166.678 

8 

0.036 

0 

-14.167 

9 

0.118 

0 

-8.275 

10 

0.057 

0 

13.531 

4.  CONCLUSIONS 

By  means  of  simulated  annealing,  the  amplitudes  and/or  phases  of  a 
hnear  combination  of  orthogonal  quasi-sinc  {sm{Nx)/ sin(x))  beams 
can  be  optimized  so  as  to  achieve  a  shaped  power  pattern  with  specified 
beamwidth  and  side  lobe  topography  and  low  ripple  while  minimizing 
the  dynamic  range  ratio  of  the  excitation  distribution.  This  proce¬ 
dure  overcomes  the  major  weaknesses  of  Woodward-Lawson  pattern 
synthesis  when  field  phase  is  irrelevant.  Side  lobe  level  and  ripple  de¬ 
pend  mainly  on  beam  coefficient  amphtudes,  dynamic  range  ratio  on 
phases.  If  the  beam  amplitudes  of  standard  Woodward-Lawson  pat¬ 
terns  are  fixed  and  only  their  phases  optimized,  the  dynamic  range 
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Figure  4.  Phase  pattern  of  the  complex  field  of  Fig.  3. 


ratio  can  improve  by  70%  or  more  (results  not  shown).  The  new  syn¬ 
thesis  method  may  be  extended  to  direct  broadcast  satellite  phased 
arrays  by  combining  it  with  the  sampling  function  method  developed 
by  Richie  and  Kritikos  [8]. 
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The  problem  and  the  solution: 

Multi-beam  antenna  arrays  have  applications  in  the  field  of  electronic  countermeasures,  in  satellite  communications,  and  in 
adaptive  nulling  [  1].  For  these  systems,  the  most  natural  method  for  synthesizing  multiple  beams  is  the  Woodward-Lawson  method; 
which  can  be  implemented  using  lossless  orthogonal  beam-forming  networks  such  as  the  Butler  matrix  [2].  The  major  objections 
to  the  Woodward-Lawson  method  have  been  its  inability  to  control  either  the  amplitude  of  the  ripple  in  the  shaped  region  of  the 
pattern,  or  the  heights  of  side  lobes  in  the  unshaped  region  [2,3] .  Another  telling  objections  is  that  the  excitations  it  prescribes  often 
have  very  high  dynamic  range  rations,  and  may  not  be  physically  realizable  in  the  presence  of  mutual  coupling  [4,5].  These 
weaknesses  ^e  ultimately  due  partly  to  the  fact  that  the  Woodward-Lawson  method  constrains  all  the  beams  with  which  it 
synthesizes  the  desired  pattern  to  have  the  same  phase,  so  that  the  far  field  pattern  is  real,  and  partly  to  its  zero  tolerance  of 
deviation  from  the  desired  field  at  the  points  at  which  it  samples  this  field.  These  major  defects  of  the  Woodward-Lawson  of 
shaped  pattern  synthesis  (lack  of  control  over  side  lobes  and  ripple,  and  high  dynamic  range  ratios)  can  be  overcome  by  perturbing 
the  field  samples  weighting  the  component  sin(Nx)/sin(x)  beams.  If  areal  field  is  not  necessary,  further  improvement  (especially 
as  regards  dynamic  range  ratio)  can  be  achieved  by  optimizing  the  phases  of  the  coefficients  of  the  component  beams. 

Example: 

Suppose  that  we  require  a  20-element  (X/2)-spaced  linear  array  to  produce  a  flat-topped  power  pattern  with  a  -3  dB  beamwidth 
of  cos-'(0.4)  and  side  lobe  levels  of  -40  dB  for  the  two  inner  lobes  and  -20  dB  for  the  others.  Fig.  1  (“real  field”  pattern)  shows 
that  the  specified  side  lobe  levels  and  ripple  levels  of  about  ±0.01  dB  can  be  achieved  by  optimizing  just  the  amplitudes  of  the 
beam  coefficients.  If  the  amplitudes  are  now  fixed  at  these  values  and  the  phases  of  the  beam  coefficients  are  optimized,  there  is 
virtually  no  change  in  the  power  pattern  but  the  dynamic  range  ratio  of  the  excitation  distribution  falls  from  17.09  to  6.83,  an 
improvement  of  60%.  Results  for  a  cosecant-squared  pattern  wil  be  shown  in  the  conference. 


Figure  1.  Flat-topped  power  patterns  obtained  by  optimizing  the 
•amplimdes  of  a  linear  combination  of  sin(Nx)/sin(x)  beams  (real  field) 
and  then  optimizing  their  phases  (complex  field). 
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Abstract 

The  widely  used  first-order  polynomial  representation  of  the  voltage  reflection  coefficient 
in  multisection  quarterwavelength  transformers  is  assumed.  The  roots  of  this  polynomial 
are  displaced  iteratively  to  obtain,  through  the  simulated  annealing  technique,  low 
variability  between  the  characteristic  impedances  of  adjacent  transformers.  Besides,  for 
each  problem  configuration,  the  bandwidth  is  maintained  close  to  its  maximum  possible 
value.  An  example  is  presented  obtaining  an  important  decrease  in  the  characteristic 
impedances  variability  of  the  transformers. 

I.  INTRODUCTION 

Quarter-wavelength  transformers  are  used  as  intermediate  matching  sections  between  a 
primarily  resistive  load  Rl  and  a  waveguide  structure  of  characteristic  impedance  Zq^Rl- 
They  are  also  used  to  connect  two  transmission  lines  with  different  characteristic 
impedances,  caused  by  a  change  in  dimensions  or  dielectric.  When  a  match  over  a  narrow 
frequency  band  is  sufficient,  a  single-section  quarter-wavelength  transformer  may  be 
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adequate.  However,  for  good  performance  over  an  extended  band  of  frequencies  it  may 
be  necessary  to  use  two  or  more  sections  in  tandem. 

Elliott  [1]  obtained  results  for  a  tandem  of  quarterwavelength  transformers  with  a  ripple¬ 
like  frequency  response  with  a  moderate  loss  in  bandwidth.  Orchard  procedure  [2]  was 
used  to  calculate  the  roots  of  the  polynomial  associated  to  the  reflection  coefficient. 
Several  solutions  were  possible  and  the  one  that  implies  less  variability  between  the 
characteristic  impedances  of  adjacent  transformers  was  chosen.  However,  using  this 
technique  we  have  not  the  guarantee  to  obtain,  for  a  desired  frequency  response,  the 
solution  with  less  variability  and  with  low  loss  in  bandwidth  in  an  optimal  sense. 

In  this  letter  we  propose  a  new  technique  to  obtain  a  desired  frequency  response  with 
low  variability  in  characteristic  impedances  of  consecutive  transformers.  This  implies 
that  the  dimensions  (or  the  dielectric)  of  adjacent  transformers  are  similar  and  the 
physical  implementation  is  easier  since  discontinuity  effects  are  minimised.  Furthermore, 
the  bandwidth  is  also  taken  into  account  optimising  its  value  according  to  input 
parameters.  This  procedure  consists  in  minimising  a  cost  function  through  the  simulated 
annealing  technique  [3]. 

n.  THEORY 

Given  N  sections  of  common  length  I  and  characteristic  impedances  Zn,  connecting  a 
mismatched  load  Rito  a  line  of  characteristic  impedance  Zo,  it  is  well  known  that,  to  first 
order,  the  voltage  reflection  coefficient  seen  at  the  end  of  the  main  line  (i.e.  at  the  input 
to  the  first  transformer  section)  is  given  by  [1]: 

n=0 
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where 


■'n+l 


^n+1  ^/i 


By  means  of  the  following  change  of  variable: 

\l/  =  -2pi 
w  =  e^'‘' 

Equation  (1)  can  be  transformed  in: 


r=ir.w" 

n=0 

which  can  be  factorised  to  give: 

n=l 


(2) 


(3) 

(4) 

(5) 

(6) 


The  synthesis  of  a  desirable  set  of  section  reflection  coefficients  Fn  can  be  obtained  by 
first  determining  a  desirable  set  of  roots  w„  and  the  study  of  the  frequency  behaviour  of 
the  input  reflection  coefficient  F  can  be  performed  by  using  (6).  By  means  of  (2)  is 
possible  to  calculate  the  characteristic  impedance  of  each  quarter-wavelength 
transformer. 

The  synthesis  procedure  consists  on  calculate  the  appropriate  reflection  coefficients  in 
order  to  obtain  a  desired  frequency  response  and  also  to  adjust  the  bandwidth  to  desired 
specifications.  The  searching  of  these  coefficients  will  be  accomplished  through  the  roots 
Wn.  In  (4)  was  shown  the  expression  of  these  roots,  but  in  general,  they  adopt  the 
following  form: 

^  (7) 
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where  the  factor  an  indicates  the  radial  displacement  of  the  nth  root  in  the  unit  circle 
representation  of  the  roots. 

With  the  assumption  of  negligible  losses,  all  section  characteristic  impedances  are  pure 
real  and  also  the  reflection  coefficients,  therefore  the  roots  must  occur  in  complex 
conjugate  pairs  or  be  in  the  real  axis.  In  this  case: 

Pn  n=l  «=1 


where  p  =1  F I  and  p  „=  I  F 


( o/  Y 

It  can  be  seen  from  the  above  expression  that,  except  for  level,  y  does  not  depend 

y  Pn 

on  the  a„  sign,  thus  there  are  2^^  sets  of  roots  w„  that  will  produce  the  same  frequency 


response. 

The  simulated  annealing  technique  [3]  allows  us  to  obtain  desirable  frequency  response 
while  maintaining  low  variability  between  the  characteristic  impedances  of  each 
transformer.  In  this  technique  there  exist  two  key  factors;  the  construction  of  a  cost 
function  that  the  algorithm  mmirnises  and  the  perturbations  introduced  to  the  variables  in 
this  cost  function,  which  are  the  roots  w„.  As  starting  solution  a  Chebyshev  distribution 
of  the  section  reflection  coefficients  was  chosen,  which  has  all  the  lobes  at  the  same 
height  and  yields  the  maximum  bandwidth  in  the  frequency  response  for  a  given  N.  The 
cost  function  used  was: 


C  =  \M„ +\Z,„-Z,f^k,  +\P,-PJtk, 


where  Mid  and  Mjo,  are  the  desired  and  the  obtained  maxima  in  the  frequency  response  of 
each  transformer,  the  second  term  is  the  maximum  difference  between  the  characteristic 
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impedance  of  two  consecutive  transformers,  which  we  define  as  variability.  Finally  Pq 
and  Pd  are  the  obtained  and  maximum  desired  percentage  of  loss  in  bandwidth 
respectively  of  the  multisection  transformer.  These  percentages  are  taken  in  relation  to 
the  bandwidth  corresponding  to  the  Chebyshev  distribution.  The  factors  k  are  weights 
used  to  give  more  importance  to  one  factor  than  to  another  in  order  to  search  for  the  best 
solution.  In  the  weights  selection  the  greater  factor  is  but  a  reasonable  relation  with  the 
others  factors  (k;  and  ks)  is  always  maintained. 

III.  EXAMPLE 

In  this  example  a  tandem  of  A  =  6  transformers  was  chosen,  the  load  at  the  end  of  the 
line  was  Rl  =  140  ohms  and  the  input  impedance  was  Zo=  60  ohms.  With  this  load 
configuration  the  selected  value  of  p^ax  was  0.1,  so  the  maxima  (Mm)  were  aU  set  to  0.1. 
The  bandwidth,  defined  as  the  frequency  range  over  which  the  magnitude  p  of  the 
reflection  coefficient  does  not  exceed  the  specified  level  Pmax,  was  controlled  by  the 
percentage  of  deviation  from  the  bandwidth  of  the  Chebyshev  pattern  (Pd  parameter  in 
(9)).  In  Table  I  there  are  shown,  for  several  values  of  Pd,  the  results  about  initial  and 
final  beamwidth  (BWi  and  BWf  respectively),  initial  and  final  variability  (W  and  Vf 
respectively)  and  the  relative  percentage  of  gain  or  loss  in  each  factor.  In  Fig.  1  it  is 
shown  the  frequency  response  when  Pd  is  chosen  to  2%.  It  can  be  seen  that  as  long  as 
the  desired  percentage  of  loss  increases,  the  variability  decreases. 


IV.  CONCLUSIONS 

The  variability  of  the  multisection  quarter-wavelength  transformers  characteristic 
impedances  was  improved  by  using  the  simulated  annealing  technique.  Furthermore  the 
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bandwidth  of  the  frequency  response  was  taken  into  account  to  acquire  more  practical 
and  realisable  solutions  for  the  construction  of  multisection  quarter-wavelength 
transformers.  Due  two  the  multiplicity  of  solutions  there  are  2*^  sets  of  roots  Wn  that  will 
give  the  same  optimal  response  function  p(pl).  Each  of  these  sets  gives  a  different 
solution  for  the  section  reflection  coefficients  and  characteristic  impedances.  In  addition, 
the  method  is  also  applicable  in  cases  where  aU  the  zones  in  the  frequency  range  have  not 
the  same  lobe  level,  i.e.,  different  pmax’s  in  different  zones  in  the  passband. 

At  present  the  extension  of  this  method  to  tapered  transmission  lines,  described  in  [4],  is 
under  study.  Finally  this  method  is  directly  applicable  to  the  problem  analysed  in  [5] 
about  multihole  directional  couplers. 
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Legends  for  Figure  and  Table. 

Figure  1.  Frequency  response  for  a  tandem  ofN=6  quarter-wavelength  transformers, 
pmax  =  0.1.  The  maximum  desired  percentage  of  loss  in  bandwidth  is  Pd  =  2%.  The  load 
used  was  Rl  =  140  and  the  input  impedance  Zo  =  60. 

Table  I.  Results  obtained  in  the  example,  N  =  6. 
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Abstract 

The  widely  used  first-order  polynomial  representation 
of  the  voltage  reflection  coefficient  in  multisection 
quanerwavelength  transformers  is  assumed.  The  roots 
of  this  polynomial  are  displaced  iteratively  to  obtain, 
through  the  simulated  annealing  technique,  low 
variability  between  the  characteristic  impedances  of 
adjacent  transformers.  Besides,  for  each  problem 
configuration,  the  bandwidth  is  maintained  close  to  its 
maximum  possible  value.  Examples  for  several 
numbers  of  transformers  and  maximum  reflection 
coefficients  are  presented  obtaining  an  important 
decrease  in  the  characteristic  impedances  variability  of 
the  transformers. 


the  solution  with  less  variability  and  with  low  loss  in 
bandwidth  in  an  optimal  sense. 

In  this  paper  we  propose  a  new  technique  to  obtain  a 
desired  frequency  response  with  low  variability  in 
characteristic  impedances  of  consecutive 
transformers.  This  implies  that  the  dimensions  (or 
the  dielectric)  of  adjacent  transformers  are  similar 
and  the  physical  implementation  is  easier  since 
discontinuity  effects  are  minimised.  Furthermore,  the 
bandwidth  is  also  taken  into  account  optimising  its 
value  according  to  input  parameters.  This  procedure 
consists  in  minimising  a  cost  function  through  the 
simulated  annealing  technique  [3]. 

n.  THEORY 


I.  INTRODUCTION 

Quarter- wavelength  transformers  are  used  as 
intermediate  matching  sections  between  a  primarily 
resistive  load  Ri  and  a  waveguide  structure  of 
characteristic  impedance  Zq  ^  Rl-  They  are  also  used  to 
connect  two  transmission  lines  with  different 
characteristic  impedances,  caused  by  a  change  in 
dimensions  or  dielectric.  When  a  match  over  a  narrow 
frequency  band  is  sufficient,  a  single-section  quarter- 
wavelength  transformer  may  be  adequate.  However, 
for  good  performance  over  an  extended  band  of 
frequencies  it  may  be  necessary  to  use  two  or  more 
sections  in  tandem. 

Elliott  [1]  obtained  results  for  a  tandem  of 
quarterwavelength  transformers  with  a  ripple-like 
frequency  response  with  a  moderate  loss  in  bandwidth. 
Orchard  procedure  [2]  was  used  to  calculate  the  roots 
of  the  polynomial  associated  to  the  reflection 
coefficient.  Several  solutions  were  possible  and  the 
one  that  implies  less  variability  between  the 
characteristic  impedances  of  adjacent  transformers  was 
chosen.  However,  using  this  technique  we  have  not  the 
guarantee  to  obtain,  for  a  desired  frequency  response, 


Given  N  sections  of  common  length  I  and 
characteristic  impedances  Zn,  connecting  a 
mismatched  load  R^  to  a  line  of  characteristic 
impedance  Zq,  it  is  well  known  that,  to  first  order, 
the  voltage  reflection  coefficient  seen  at  the  end  of 
the  main  line  (i.e.  at  the  input  to  the  first  transformer 
section)  is  given  by  [1]: 


where 


n=:0 


r  = 


(1) 

(2) 


By  means  of  the  following  change  of  variable; 


If/ =  -2p  I  (3) 

w=e^'‘'  (4) 


Equation  (1)  can  be  transformed  in: 
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r=|;r,>v-  (5) 

n=0 

which  can  be  factorised  to  give: 

r  =  r^f[iw-w„)  (6) 

rt=i 

The  synthesis  of  a  desirable  set  of  section  reflection 
coefficients  r„  can  be  obtained  by  first  determining  a 
desirable  set  of  roots  w„  and  the  study  of  the  frequency 
behaviour  of  the  input  reflection  coefficient  F  can  be 
performed  by  using  (6).  By  means  of  (2)  is  possible  to 
calculate  the  characteristic  impedance  of  each  quarter- 
wavelength  transformer. 

The  synthesis  procedure  consists  on  calculate  the 
appropriate  reflection  coefficients  in  order  to  obtain  a 
desired  frequency  response  and  also  to  adjust  the 
bandwidth  to  desired  specifications.  The  searching  of 
these  coefficients  will  be  accomplished  through  the 
roots  w„.  In  (4)  was  shown  the  expression  of  these 
roots,  but  in  general,  they  adopt  the  following  form: 


to  the  variables  in  this  cost  function,  which  are  the 
roots  w„.  As  starting  solution  a  Chebyshev 
distribution  of  the  section  reflection  coefficients  was 
chosen,  which  has  all  the  lobes  at  the  same  height 
and  yields  the  maximum  bandwidth  in  the  frequency 
response  for  a  given  N.  The  cost  function  used  was: 

(9) 

where  Mjd  and  Mio,  are  the  desired  and  the  obtained 
maxima  in  the  frequency  response  of  each 
transformer,  the  Second  term  is  the  maximum 
difference  between  the  characteristic  impedance  of 
two  consecutive  transformers,  which  we  define  as 
variability.  Finally  Po  and  Pj  are  the  obtained  and 
maximum  desired  percentage  of  loss  in  bandwidth 
respectively  of  the  multisection  transformer.  These 
percentages  are  taken  in  relation  to  the  bandwidth 
corresponding  to  the  Chebyshev  distribution.  The 
factors  kj  are  weights  used  to  give  more  importance 
to  one  factor  than  to  another  in  order  to  search  for 
the  best  solution.  In  the  weights  selection  the  greater 
factor  is  kz  but  a  reasonable  relation  with  the  others 
factors  (k/  and  ks)  is  always  maintained. 


(7) 


where  the  factor  indicates  the  radial  displacement  of 
the  nth  root  in  the  unit  circle  representation  of  the 
roots. 

With  the  assumption  of  negligible  losses,  all  section 
characteristic  impedances  are  pure  real  and  also  the 
reflection  coefficients,  therefore  the  roots  Wn  must 
occur  in  complex  conjugate  pairs  or  be  in  the  real  axis. 
In  this  case: 


'  P'' 


N 


n=l 


(8) 


where  p  =1  F I  and  p  „  =  I  F  „1. 

It  can  be  seen  from  the  above  expression  that,  except 


for  level. 


Pn 


does  not  depend  on  the  a„  sign. 


thus  there  are  2^  sets  of  roots  that  will  produce  the 
same  frequency  response. 

The  simulated  annealing  technique  [3]  allows  us  to 
obtain  desirable  frequency  response  while  maintaining 
low  variability  between  the  characteristic  impedances 
of  each  transformer.  In  this  technique  there  exist  two 
key  factors:  the  construction  of  a  cost  function  that  the 
algorithm  minimises  and  the  perturbations  introduced 


m.  EXAMPLES 


A.  Example  1 

In  this  example  a  tandem  of  =  6  transformers  was 
chosen,  the  load  at  the  end  of  the  line  was  Rl  =  140 
ohms  and  the  input  impedance  was  7^=  60  ohms. 
With  this  load  configuration  the  selected  value  of 
Pmax  was  0. 1,  so  the  maxima  (My)  were  all  set  to  0. 1 . 
The  bandwidth,  defined  as  the  frequency  range  over 
which  the  magnitude  p  of  the  reflection  coefficient 
does  not  exceed  the  specified  level  pn,„,  was 
controlled  by  the  percentage  of  deviation  from  the 
bandwidth  of  the  Chebyshev  pattern  (Pj  parameter  in 
(9)).  In  Table  I  there  are  shown,  for  several  values 
of  Pd,  the  results  about  initial  and  final  beamwidth 
(BW'i  and  BWf  respectively),  initial  and  final 
variability  (Vj  and  Vf  respectively)  and  the  relative 
percentage  of  gain  or  loss  in  each  factor. 


Prt 

BWi 

BWf 

% 

Vi 

Vf 

% 

2.0 

4.92 

4.84 

-1.62 

18.51 

12.55 

-32.18 

1.5 

4.92 

4.88 

-0.81 

18.51 

13.70 

-26.01 

0.5 

4.92 

4.92 

0.00 

18.51 

15.46 

-16.49 

Table  1.  Results  obtained  in  the  example,  N  =  6. 

In  Fig.  1  it  is  shown  the  frequency  response  when  Pj 
is  chosen  to  2%.  It  can  be  seen  that  as  long  as  the 
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desired  percentage  of  loss  increases,  the  variability 
decreases. 


Figure  1.  Frequency  response  for  a  tandem  of  N  =  6 
quarter-wavelength  transformers,  =  0.1.  The  maximum 
desired  percentage  of  loss  in  bandwidth  is  Pj  =  2%.  The 
load  used  was  7?^  =  140  and  the  input  impedance  =  60. 

B.  Example  2. 

In  this  case  the  parameters  of  the  problem  were:  N  =7, 
Rl  =  100  ohm,  Zo  =  60  ohm  and  p,„ax  =0.05.  The 
results  for  three  values  of  Pj  are  shown  in  Table  11. 
The  initial  and  final  frequency  response  are  shown  in 
Figure  2. 


Pd 

BW, 

BWf 

% 

Vi 

Vf 

% 

2.0 

4.96 

4.88 

-1.62 

6.98 

5.34 

-23.42 

1.5 

4.96 

4.92 

-0.81 

6.98 

5.33 

-23.52 

0.5 

4.96 

4.96 

0.00 

6.98 

5.49 

-21.38 

Table  11.  Results  obtained  in  the  example,  N=  7. 


This  example  gives  similar  results  for  each  Pj,  all  over 
the  20%  of  improvement  in  the  variability.  The  most 
important  factor  in  this  case  is  the  greater  difference 
between  \Rl-ZJARi+Z„\,  which  is  the  reflection 
coefficient  at  =  0,  and  the  desired  pn,ax-  Other 
examples  performed  showed  the  same  behaviour,  the 
smaller  is  this  difference,  the  better  is  the  improvement 
in  the  variability. 


Figure  2.  Frequency  response  for  a  tandem  of  N  =  1 
quarter-wavelength  transformers,  pn^j,  =  0.05.  The 
maximum  desired  percentage  of  loss  in  bandwidth  is  Pj  = 
2%.  The  load  used  was  /?£,  =  100  and  the  input  impedance 
Zo  =  60. 

IV.  CONCLUSIONS 

The  variability  of  the  multisection  quarter- 
wavelength  transformers  characteristic  impedances 
was  improved  by  using  the  simulated  annealing 
technique.  Furthermore  the  bandwidth  of  the 
frequency  response  was  taken  into  account  to  acquire 
more  practical  and  realisable  solutions  for  the 
construction  of  multisection  quarter-wavelength 
transformers.  Due  two  the  multiplicity  of  solutions 
there  are  2^^  sets  of  roots  w„  that  will  give  the  same 
optimal  response  function  p(Pl).  Each  of  these  sets 
gives  a  different  solution  for  the  section  reflection 
coefficients  and  characteristic  impedances.  In 
addition,  the  method  is  also  applicable  in  cases 
where  all  the  zones  in  the  frequency  range  have  not 
the  same  lobe  level,  i.e.,  different  Pnux’s  in  different 
zones  in  the  passband. 

At  present  the  extension  of  this  method  to  tapered 
transmission  lines,  described  in  [4],  is  under  study. 
Finally  this  method  is  directly  applicable  to  the 
problem  analysed  in  [5]  about  multihole  directional 
couplers. 
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